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I. INTRODUCTION 

Flexible polymerized 2-dimensional films (tethered or polymerized membranes) [1] have 
very interesting statistical properties (for a review see [2-4]). In these objects there is a 
competition between entropy which favors crumpled or folded configurations as for polymers, 
steric interactions (self-avoidance) which tend to swell the membranes, and bending rigidity 
which favors flat configurations. Internal disorder, inhomogeneities and anisotropy may also 
play an important role, that we shall not discuss here (see the chapters 10-12 in [3] and [4] 
for a recent review of these effects). 

If one does not take into account self-avoidance, theoretical arguments (mean-field and 
renormalizat ion-group calculations) and numerical simulations show two phases: (1) A high- 
temperature/low-rigidity crumpled phase, where the membrane is crumpled with infinite 
Hausdorff dimension, and where bending rigidity is irrelevant; (2) a low-temperature (flat) 
phase with large effective rigidity and Hausdorff dimension two [5, 6]. 

In the high-temperature (crumpled) phase, steric interactions (self-avoidance) are phys- 
ically relevant, and will swell the membrane, as for polymers. Two scenarios are possible: 
Either the membrane is fiat with Hausdorff dimension two (as for high bending rigidity), or 
it is crumpled swollen with Hausdorff dimension larger than two. For large imbedding di- 
mension d a Gaussian variational approximation can be argued to become exact, predicting 
a Haussdorff dimension dn = d/2. It is non-trivial whether this swollen phase exists down 
to d = 3. Numerical simulations indicate that only a fiat phase exists in d = 3; for details 
see the discussion in [3, 4]. However, simulations are for rather small systems. It therefore 
remains important to have a solid theoretical understanding. 

The theory in question is a generalization of the Edwards model for polymers [10, 11]. It 
was proposed in [7-9], and is a model of self-avoiding manifolds (or membranes), hereafter 
denoted the SAM model. It is amenable to a treatment by perturbation theory (in the 
coupling constant of steric interactions) and to a perturbative renormalization group analysis 
which leads to a Wilson-Fisher like e-expansion for estimating the scaling exponents and the 
critical properties of the swollen phase. 

This SAM model is quite interesting at the theoretical level for several reasons: 

1. It can only be defined as a non-local field-theory over the internal 2-dimensional space 
of the manifold, with infinite-ranged multi-local interactions. Therefore the applicabil- 
ity of renormalization theory and of renormalization group techniques is a non-trivial 
issue. A proof of perturbative renormalizability to all orders was finally given in 
[12, 13]. 

2. The model is in fact defined through a double dimensional continuation, where both 
the dimension of space d and the internal dimension D of the manifold are analytically 
continued to non-integer values. The physical case of two-dimensional membranes is 
always in the strong-coupling regime where the engineering dimension of the coupling, 
e = 2D — ^^^d is e = 4 for any space dimension d. 

3. The analytical study of this model at the non-perturbative level is still in its infancy, 
since it is a technically quite difficult problem. A first step was made by the two 
present authors for the large orders of perturbation theory in [14]. It is this issue of 
the large-order asymptotics of the SAM model that we treat in this paper. 

For quantum mechanics [15-17] and for local quantum field theories [18] (such as the Laudau- 
Ginzburg- Wilson 0^ theories) the large-order asymptotics of perturbation theory are known 
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to be controlled by (in general complex) finite-action solutions of the classical equation of 
motion called "instantons" . More precisely the large-order asymptotics are described by 
semi-classical approximations around these instantons. We refer for instance to [19] for a 
review of this "instanton calculus" . 

In [14] we have shown that similarly for the SAM model there exists an instanton which 
controls its large-order asymptotics. This instanton is a scalar field configuration in the 
external ci-dimensional space, which extremizes a highly non-local effective-action functional, 
and which cannot be computed exactly. We also showed that remarkable simplifications 
occur in the large-d limit, which suggests that a systematic 1/d expansion can be constructed 
to study the instanton, but also that already the first 1/d correction to the large-d limit is 
plagued with infrared (IR) divergences whose origin was unclear. In [14] we only studied the 
instanton at the classical level, i.e. the (non-local) equation of motion and the properties of 
its solution, the instanton. 

In this article we present the full semi-classical analysis of the instanton for the SAM 
model, derive its connection with the large-order asymptotics, and study the UV divergences 
and renormalization necessary for the instanton. For this purpose, many new calculational 
techniques had to be developed, hence the length of the paper and its technical character. 
More precisely, the main new results are: 

1. We first show in much more details than in [14] how the instanton emerges from the 
functional integral which defines the continuum SAM model. In particular we treat 
properly and carefully the zero-modes for the instanton, how the contour of functional 
integration has to be deformed in the complex saddle-point method, as well as various 
normalization problems for the functional integration. This is done in Sect. III-A,B. 

2. Using this, we obtain the contribution of the fluctuations around the instanton in the 
semi-classical approximation as the determinant of a non-local kernel operator in ci- 
dimensional space, and derive the normalization factor for the large-order asymptotics 
(Sect. III-C,D,E). 

3. We analyze completely the UV divergences of this determinant, and show that in the 
renormalized theory these UV divergences for the instanton determinant factor are 
canceled by the one-loop perturbative counterterm of the renormalized theory, making 
the final asymptotics UV finite. This is an important check of the consistency of the 
SAM model, since the original proof of renormalizability is only valid in perturbation 
theory. (In a field theoretic language it is not a background-independent proof). 
This is done in Sect. IV. Our argument is based on the extension of the perturbative 
renormalizability argument to the general case of ensembles of interacting manifolds 
in an external background potential. 

4. In [14] the instanton equation was solved within a variational approximation. In 
sect. V we study how this approximation can be applied to the exphcit calculation of 
the instanton determinant factor. We first show that a direct variational calculation 
gives a result which is too naive, and does not take properly into account the UV 
fluctuations. We then propose a systematic framework to construct an expansion 
around the variational approximation, developing ideas that we proposed in [14]. We 
then show that this framework gives the leading term for the instanton determinant 
factor in the large-ci limit. 
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5. We are thus able to construct a systematic 1/d expansion for the instanton calculus, 
and show that this expansion is well defined as long as the SAM is super renormaliz- 
able, i.e. e > (no UV divergences in perturbation theory, apart from vacuum energy 
terms). The leading and first subleading terms are computed explicitly for the deter- 
minant factor and the normalization factor of the zero mode of the instanton. These 
calculations involve a new non-trivial diagrammatic expansion. This is done in Sect. VI 

6. Finally we study the 1/d expansion for the renormalized theory at e = 0. We show 
that, at variance with the instanton calculus for local field theories, some subtle issues 
arise for the SAM model. Indeed, we show that already at leading order in the limits 
d oo and e ^ do not commute, and that some care is needed in order to obtain 
the instanton determinant factor for the renormalized theory at large d. We then show 
that the subleading terms of the \/d expansion are plagued with IR divergences at 
e = 0, generalizing the results of [14]. We analyze completely these IR divergences 
at the first subleading order, and show that many compensations occur, leaving a 
single IR-singular term associated with a single eigenmode for the fiuctuations around 
the instanton, namely the unstable eigenmode generated by global dilation for the 
instanton. This analysis of the renormalized theory is done in Sect. VII. 

To summarize, we have performed a non-trivial check of the consistency of the model, in 
particular of its renormalization, in a non-perturbative regime, and we have developed the 
tools to compute the large-order asymptotics of the SAM model. 

Appendices contain more technical computations and details about the normalizations. 
In particular in appendix C we explicitly check that in the special case of the self-avoiding 
polymer (D = 1) we recover the large-order asymptotics of the Edwards model obtained by 
field theoretical methods (using instanton calculus and the well known equivalence between 
the Edwards model and the 0(n = 0) 0^ Landau Ginzburg model [20, 21]). This provides a 
check of the consistency of the SAM model. 

II. THE MODEL 

A. The non-interacting manifold 

First we define the model for the Gaussian non-interacting manifold (free or phantom 
manifold). Of course this model reduces to a massless free field, but we reconsider it closely 
in order to fix properly the normalization for the measure and for the definition of the 
observables, and for the treatment of the zero modes. 

1. The model and its action 

We consider a manifold M. with a finite size, as a closed /^-dimensional manifold M., with 
a fixed internal (or intrinsic) Riemannian structure, given by a metric tensor g = g^y{x). 
X = {x^\ fi = 1, . . . ,D) describes (a system of) local coordinates on Ai. From now on the 
internal volume of Ai, Vol(7W) and its internal size L are defined as 





(2.1) 
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with g = det[g]. The manifold is embedded in external (or bulk) rf- dimensional Euclidean 
space M'^. This embedding is described by the field r = {r"; a = 1, . . . ,d} 

M^R'^ , X ^ r(x) (2.2) 

We shall use dimensional regularization in this paper by considering that the internal di- 
mension D of the manifold is < -D < 2 non-integer. See the reference paper [13] for a more 
precise discussion of how we can define a finite membrane within dimensional regularization. 
In practice we can restrict ourselves to the case of a square D-dimensional torus of size L, 
Td = (M^/(L ■ Z)^), with fiat metric g^^ = S^u- 

We first consider the free non-interacting manifold (phantom membrane). The manifold 
may fiuctuate freely in external c?- dimensional space. Its free energy is given by the Gaussian 
local elastic term 5*0 , which is the integral of the square of the gradient of the field r 

So[r] = jj^^^'f = j d^'x^^^g^^d.r-d.r . (2.3) 

This is nothing but the Euclidean action for a free massless field (with d components) living 
on M.. The manifold may (and does) freely intersect itself, as does a random Brownian walk 
in d < 4 space dimensions. 



2. The partition function 

The partition function for the free manifold is thus given by the functional integral 

Zo = j P[r]e-^°[^1 (2.4) 

where P[r] is the standard functional measure for the free massless field r (see Appendix A 
for details and the normalization used in this paper). 

There is an infinite factor in Zq (the volume of bulk space Vol(M'^)) coming from the 
translational zero mode of the manifold. This can be isolated by choosing a specific point xq 
on the manifold and a specific point rg in bulk space, and by defining the partition function 
Zq for a marked manifold 

Zo = Zo(ro) = y"l?[r](5'^(r(xo)-ro)e-^»[^] (2.5) 

Zo{ro) is infra-red (IR) finite and does not depend on the choice of ro or of xq. We have 
formally 

= J d'ro Zo(ro) = VoliR") Zq . (2.6) 

The partition function Zq is found to be related to the determinant of the Laplacian 
operator over M. through 

r 1 -'^/2 

Zo(ro) = det' [-A] ■ 27r/Vol(7W) , (2.7) 

where det' [—A] is the product of the non-zero eigenvalues of (minus) the Laplacian operator 
A = g~^^^d^g^'^dy on A4.. Vol(AI) = J d^xy^ is the internal volume of the manifold. This 
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last term comes from the proper treatment of the translational zero mode of the Laplacian 
(see Appendix A). 

The determinant det' [—A] is ultra-violet (UV) divergent, and is defined through a zeta- 
function regularization (for a manifold A4 with non-integer dimension D 7^ 1 or 2 this is 
equivalent to dimensional regularization) 

log(det'[-A]) = -C'(0) , C(s) = tr'((-A)-^) . (2.8) 

The zeta-function ({s) is defined by analytic continuation from Re(s) large, tr' means the 
trace over the space orthogonal to the kernel of A. C{s), scales with the size of Ai as 

C(.) = [Vol(A^)]^^/^ C(^) , (2.9) 

where the "normalized zeta-function" ({s) depends on the shape of the manifold but not on 
its size (scale invariance). In the absence of a conformal anomaly, as this is the case for the 
generic case of D non-integer we have the exact identity 

C(0) = -1 . (2.10) 

(This factor comes from the contribution of the subtracted zero mode in the determinant). 
Hence the partition function reads 

Zo(ro) = [Vo\{M)r'-^ 

The last term is a "form factor" depending on the shape of Ai. 

For 2-dimensional manifolds {D = 2), the conformal anomaly gives an additional scale 
factor of the form \LfiQ\^'^^^, where L = Vo\{A4y^^ is the size oi /A, Hq the regularization 
mass scale, required to define properly the measure in the functional integral (see Appendix 
A), and x the Euler characteristics of the membrane. We shall not discuss this any further, 
since this is not relevant for the problem treated here, where we consider manifolds with 
D <2. 

B. The interacting self-avoiding manifold 

1. The action 

The steric self-avoiding interaction is introduced by adding a 2-body repulsive contact 
interaction term of the form 

/ / SM^) - r(y) = / d^xv^ / d^v^5'(Kx) - r(y)) 

JxJy JM JM 

(where 6'^{r) is the Dirac distribution in the external space ) to the action, which is now 

6 > is the self- avoidance coupling constant. This is similar to what is done in the Edwards 
model for polymers. 



.C'(o) 
"2^ 



(2.11) 
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2. The partition functions 

The partition function for the self-avoiding manifold is 

Z(ro,6)= y"p[r](5'^(r(xo)-ro)e"^[^'^l ; Z{b) = Jd\Z{,o,b)- (2-13) 

These partition functions are defined in perturbation theory, within a dimensional regular- 
ization scheme, i.e. by analytic continuation in the internal dimension D. 

If the internal coordinate x has engineering dimension 1, the external coordinate r has 
engineering dimension uq given by (i.e. [r] ~ [x]*^") 

z/o = i2-D)/2 (2.14) 

and the coupling constant b has engineering dimension — e (i.e. [b] ~ [x]^*^) with 

e = 2D -{2- D)d/2 (2.15) 

As usual in polymer and membrane problems, we shall consider mainly the normalized 
partition function 3(^), defined by the ratio of the interacting partition function for the 
interacting manifold M., divided by the partition function for the same manifold M., but 
free. 

m = Z{b)/ZQ = Z(ro,6)/Zo(ro) . (2.16) 

Let 

L = (yo\{M)Y^^ ; Vo\{M) = I (2.17) 

JM 

be the internal size of the manifold. The normalized partition function has a perturbative 
series expansion in powers of 6, of the form 

oo 

3(6) = 1 + 5^3fc [bUf (2.18) 
fc=i 

where the coefficients 3fc depend only on the shape of the manifold, on its internal dimension 
and on the external dimension d. These coefficients are given by the expectation value in 
the massless free theory defined by the free action 5*0 of the bi-local operators corresponding 
to the interaction term 

■■■II ( n ^'(^(^') - ^(^^))>o ' (2-19) 

xi,yi J Jxk,yk i=l k 



k\ 2^ 

with (. . . )o the expectation value w.r.t. 5o[r], see (2.3) 



3. Observables and correlation functions in external space 

We shall be mainly interested in correlations functions which correspond to observables 
which are global for the manifold (i.e. do not depend on the internal position of specific 
points on the manifold), but which may be local in external space (i.e. do depend on the 
position of specific points in the external space). These observables are the simplest ones. 
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In particular for the case D = 1 (polymers) these observables have a direct interpretation 
in terms of correlation functions of local operators in the corresponding local field theory in 
external (i- dimensional space. 

The observables involve the manifold density p(r). We define the manifold density at the 
point ri, p(ri), as the functional of the field r(x) 

p(ri;r) = [ 5\r{x) - r,) . (2.20) 

Jx 

The A^-point density correlator 7^'^^)(ri ■ ■ ■ r^v) for the interacting manifold is defined as 

N 

7^W(rl...r^;6) = J V[r] H pir,; r) e-'^'^'K (2.21) 



i=l 



Obviously the one-point density correlator is related to the partition function (for a one-point 
marked manifold) by 

7^«(ro;&) = Vo\{M)Z{ro;b) . (2.22) 

Ratios of density correlators define expectation values of densities. For instance, the e.v. 
(expectation value) of a product of density operators p(rj) for a manifold constrained to 
be attached to a point Tq is the ratio 

(p(ri)---p(r;v))ro;. = 7^(^+^)(ro,rl■■■r;v;6)/7^«(ro;&) • (2.23) 

As for the partition functions, it is more convenient to normalize the density correlators 
with respect to the partition function for the free manifold. We thus define the normalization 
for the normalized density correlators, by 

^(^^ ih---rN;b) = 7^(^) ir,---r^;b)^Vo\iM)Zo. (2.24) 

In particular the normalized 1-point correlator coincides with the normalized partition func- 
tion 

9^('Hri;&) = m (2.25) 

and is independent of ri. 

These observables have a perturbative series expansion in the coupling constant b. In 
particular they scale with the size L of the manifold as 

7^(^) {r,---rN;b,L) = . n^^) (n L-"" , ■ ■ ■ , L-"" ; bL') (2.26) 

and 

SK(^)(ri---r^;6,L) = L^^-^^^^"^) ■ ^n^^^Ci ■ ■ ■ , r^v &i^') (2.27) 

4. Global quantities and gyration radius moments 

We define the moments of order k for the gyration radius (in short the k-th gyration 
moment) Rgyr of the manifold by 

Jxi Jxo. 
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The standard gyration radius is Rgyr = y Rgyl- The expectation value T^gyl of the k-th 
gyration moment R^l for the interacting manifold is thus (for k > 0) 

<i = (R£1) = ^^^1^ 1 1 Ki - ^^1' ^^'H^i' b) . (2.29) 

C. UV divergences and perturbative renormalization 

Using dimensional regularization, the perturbative expansion for the partition function 
and the observables is known to be UV finite for 

2D 

< D <2 and D <e i.e. d < - — — . (2.30) 

As long as we deal with finite-size manifolds (L < oo), perturbation theory is free from 
infra-red divergences (which occur for infinite manifolds since perturbation theory is made 
around the free- manifold theory, which is a free massless scalar field in D < 2 dimensions). 
The perturbative expansion suffers from short-distance (UV) divergences when 

D < € . (2.31) 

These UV divergences come from the short- distance behavior of the expectation values which 
appear as integrands in the integrals, when the distance between several points Xj and yj 
goes to zero. Using dimensional regularization these divergences appear as poles in e {d 
being fixed), or equivalently as lines of singularity in the {d, D) plane. 

As proved in [13], these UV divergences can be studied within a multi-local operator 
product expansion (MOPE) which generalizes Wilson's OPE of local field theory. As a 
consequence, these UV divergences are proportional to the insertion of multi-local operators, 
and are amenable to renormalization theory. 

The MOPE formalism and dimensional analysis show that for < e < there is a finite 
number of divergences, with poles at 

e = D/n , n G N+ . (2.32) 

These divergences are proportional to insertions of the identity operator 1 (with dimension 
0). The model is super-renormalizable for e > and these divergences are subtracted by 
adding to the action a local counterterm proportional to the volume of the manifold (i.e. to 
the integral of the identity operator 1). 

AS(r) oc /" 1 = Vo\{M) . (2.33) 

J X 

These divergences and the corresponding counterterm are constant terms, independent of the 
configuration of the manifold, i.e. of the field r, and they disappear in the observables given 
by ratios of correlators such as the e.v. (p(i'i) • ■ ■ p(i'Af))ro and the normalized correlators 

For e = the model has an infinite number of divergences. These divergences are pro- 
portional to the insertion of the two operators present in the original action S. This means 
that the theory is renormalizable, and that it can be made UV-finite by adding to the action 
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counterterms of the same form than those of the original action. In other words, one can 
construct in perturbation theory a renormahzed action 

SAr,b„^) = ^ / (Vr)= + '-^^^ f f A'M - r(y)) , (2.34) 

t/ X o y. o y 

where 6r is the dimensionless renormahzed couphng constant, Z{h^) and Zf,{bj.) the wave- 
function and the couphng-constant renormahzation factors, and fi is the renormahzation 
momentum scale, while r(x) is now the renormahzed field. This renormahzed action is such 
that the renormahzed correlation functions 

N 

7ef)(ri---r^;6r,/i) = V[r] l[pir,; r) e-'^^'''-'^^ (2.35) 

^ i=i 

have a perturbative series expansion in bj- which is UV finite for e > and stays finite for 
e = 0. For a finite manifold with size L the renormahzed perturbation theory is still IR 
finite. 

From the standard arguments of renormahzation group (RG) theory the renormahzed 
theory describes the universal large-distance scaling behavior of self-avoiding manifolds. One 
can write RG equations which tell how the observables scale with the size of the manifold 
for e > 0. When expressed in terms of the renormahzed observables and the renormahzed 
coupling, these RG equations have a regular limit (at least in perturbation theory) as e — 
0+. As a consequence one can construct an e-expansion a la Wilson-Fisher for the scaling 
exponents. 



D. Effective non-local model in external space 

As shown in [14], to study the large-order behavior of the SAM model as well as its 
large-d behavior, it is necessary to reformulate the model as an effective non-local model 
for an auxiliary composite field V{r) in the external ti-dimensional space. We recall this 
reformulation. 



1. Auxiliary fields and effective action 

First we recall the auxiliary field p(r) (local manifold density) defined in (2.20), 



P(r) 



5'^(r(x) - r) 



(2.36) 



and its conjugate field ^(r), which is the Lagrange multiplier for the above constraint, such 
that 

' 5"(r(x) 



V[V]V[p] exp (^J d'^rV{r 



p(r) - 5'^(r(x) - r) 



M 



(2.37) 



p is a real field, while V is imaginary, and has to be integrated from —ioo to +ioo in the 
functional integral. Equivalently the functional measures for p and V are formally 



l[dp{r 



l[dp{r 



V[V] 



n^^. (2.38) 



2i7r 
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We now insert (2.37) in the functional integral. Since the interaction term can be written 



as 



5^(r(x) - r(y)) 



X 



P(r) 



(2.39) 



the integral over the field p is Gaussian and can be performed explicitly. We obtain for the 
partition function 



Z{h) 



V[r]V[V] exp 



-(Vxr)2 + y(r) 



1 

2& Jr 



V{r) 



Note that the functional measure 'D\V] over V[r\ is now normalized so that 



/PMexp(l/nr) 



(2.40) 



(2.41) 



and depends explicitly on the coupling constant h. 

This functional integral describes a free (not self-interacting) manifold fiuctuating in a 
random annealed potential V{r). This is a simple generalization of the reformulation of the 
SAW problem into a random walk in a random annealed potential. 

Now we integrate over the field r(x) and define the effective free energy Fa4[1^] for the 
non- interacting (phantom) manifold M. in the external potential ^(r) by 



(2.42) 



exp(-FM[y]) = Jv[r]exp - j^(^(y^rf + V{r) 
We are left with the effective action for the field V , Sm\^]i which is given by 

Sm[V] = Fm[V] - Y^Virf 



(2.43) 



and is a non-local functional of the potential V. The partition function is now given by a 
functional integral over the potential V alone 



Z{b) 



V[V] exp 



Fm[V] + 



1 

2b J, 



V{r) 



V[V] exp[-SM[V]] 



(2.44) 



2. Correlation functions for global observables 



The same transformation can be used to compute the density correlators 7^^^^ and the 
corresponding correlation functions as expectation values of observables with the effective 
action S[V]. Indeed inserting a density operator p(r) in the original functional integral (2.13) 
over r(x) amounts to insert a functional derivative with respect to the conjugate field V{r) 
in the functional integral (2.44) over V{r). 



so that 
7^(^)(rl,■■■ ,rAr;fe) 



P(r) 



V[V] exp[-FM[V]] 



6 



6V{r) 



6 



exp 



2b I 



(2.45) 



(2.46) 
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For instance the 2-point correlator is 



n^'\h,r,-b) = J V[V] 



= ^Zib){V{r,)V{r,)) + ^5^1-^2) (2.47) 

Similarly, for the moments of order k of the gyration radius (defined by (refdef GyrMom) ) 
we get (for k > 0) 

"^'S. = '^y^. jjj'^-'^t {nn)V(r,)) , (2.48) 
where ( ) denotes the average over V with the effective action ^xfV] given by Eq. (2.44). 

III. LARGE ORDERS OF PERTURBATION THEORY AND INSTANTON CALCULUS 
A. Instanton and large orders in Quantum (Field) Theory 

1. Instanton semi-classics 

To fix our normalizations let us first recall the basics of instanton calculus in quantum 
mechanics and quantum field theory. We consider a model defined by the functional integral 
over a field 0(r) with a classical action S[(p] and a (dimensionless) coupling constant g. The 
partition function is 

Z = j e-i^[^l . (3.1) 

The functional measure ^^[0] over is defined from the so-called DeWitt metric G on classical 
field configuration (super) space. We choose it to be local and translationally invariant, so it 
must be of the form 

II ■ II 2 is the L2 norm. This metric depends explicitly on an (arbitrary) normalization mass 
scale /io = /Uo/y^. The corresponding measure over field (super)space is (formally) = 



J|^d0(r)\/det G. It is such that 



/ 



p[<^]e-2f Jr<^^ = 1 . (3.3) 



(Note that the factor of g has been introduced for convenience, to have the same functional 
dependence on g for the measure in (3.3) and the Boltzmann factor in (3.1).) 

We assume that there is a classical vacuum (field configuration) 0o which minimizes the 
action S", which is constant (0o('') = 0o) £^nd which is unique (no zero modes around the 
classical vacuum). In the semi-classical approximation the contribution of 0o fo fhe partition 
function is simply 



classical vacuum 



e4^I^o] (Det[y'[0o]^a)"'^' > (3-4) 
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where S" is the Hessian operator, with kernel 

^"['^Kri,r2) = tJ^Wt-,- (3-5) 
6(j){ri)6(f){r2) 

Now we assume that there are also instanton configurations which contribute to the func- 
tional integral. These instantons are non-constant field configurations (pmsti"", z"") which 
are classical solutions of the field equations, and thus local extrema of the action S, i.e. 
'S"[0inst] = 0, with a finite action S'inst = 5'[0inst]- In general, the set of instantons with 
action S'inst is a finite-dimensional subspace, called the instanton moduli space. We denote 
z = {z°',a = l,m} a (local system of) collective coordinates on the m-dimensional mod- 
uli space of the instantons with action S'inst- The collective coordinate z must include the 
position of the instanton rinst {d moduli), its size if the action S is scale invariant, and in 
addition the internal degrees of freedom of the instanton if needed. 

The contribution of the instanton to the functional integral is also given by a semi-classical 
formula. We must separate the integration over the instanton moduli space Aiz from the 
integration over the field fluctuations transverse to the moduli space A4z, since the Hessian 
S"'[0inst] has now 5inst zero-modes da4>mst = The moduli space integration is then 

done explicitly. For that purpose, we must consider the restriction of the metric G to the 
instanton moduli space Ai^- The corresponding metric tensor hat in the coordinate system 
z is defined by 

||d(/.instir = dz'^dz'K.iz) ; d0inst = ^d^^ (3.6) 
where d^inst is an instanton fiuctuation. Hence the metric on moduli space is 



hah(z) 



d(t>: 



inst 



dz" 



d(pmst\ /io /" ^0inst(r, 2)<90inst(r,2) 



dz'' J 2ng J, dz'' dz^ 



(3.7) 



and the corresponding measure is d/i(z) = d^'z ^ydet{h). The contribution of the fiuctuations 
orthogonal to the moduli space Aiz is evaluated by the saddle-point method. The final result 
for the contribution of the instanton to the partition function is 



Mz 



— 1 /2 

d^Za Vdet{hab{z)) e-i^['^'-'l (def [5"[0inst]//i^])^ , (3. 



where det' [S"[0inst]] is the product of the non-zero eigenvalues of S"[0inst]- The det{h) gives 
a power of the coupling constant g'"^^"^, where m is the number of instanton zero-modes. 

Similarly, let us now consider the expectation value for an observable 0[(f)], for instance 
a product of fields O = 0(i'i) ■ ■ ■ (f){fn)- The expectation value is given by 

(O) = ^ 1 1?[0]O[0]e-^^[*l (3.9) 

The contribution of the (translationally invariant) classical vacuum to (O) is simply 

classical vacuum 

(O) < O[0o] . (3.10) 
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The contribution to (O) of the instanton 0inst) is obtained from 

(O) ^ ^ — r . (3.11) 

/I?[^](e-i^[*"+^l+e"5^[<^--+^l) 

This expression is rather symbohc, since we have not written the integral over the 0-mode 
of the instanton. Since Sq < S'inst, we have So/g Sinst/g, ioi g —>■ 0. Thus the leading 
term of (3.11) is given by (3.10), and the subleading one is the contribution of the instanton, 
which (up to exponentially small terms) reads 

(O) [ d™z„ Vdet(/i(^)) (o[0i„stN] - 0[4>o]) e4(^[*-l-^M) 

/ Detq^--[0i„3t]//ig] V^^^ 
H Det[S"[0o]//i^] j • ^ 

One can check that the /ig dependence disappears (remember that the moduli metric h 
depends on /io) as long as there is no scale anomaly coming from UV-divergences in the 
ratio of the two determinants of the Hessians. 



2. Large orders of perturbation theory and instantons 

We now recall the basic argument which shows how the large orders of perturbative series 
obtained by functional integrals are related to instantons. 

We assume that the observable 3(5') has a series expansion for small positive coupling con- 
stant g and is in fact an analytic function of the coupling constant (7 in a complex neighbor- 
hood of the origin away from the negative real axis (i.e. for \g\ small enough, |Arg(5f)| < vr), 
but with a discontinuity along the negative real axis (|Arg((7)| = tt). 

Its asymptotic series expansion is written as 

00 

3(5) = E^/^^'- (3-13) 

fc=0 

The large order (large /c) asymptotic behavior of the coefficients 3fc can be estimated by 
semi-classical methods. Indeed, using a classical dispersion relation, 3fc can be written as a 
Mellin-Barnes integral transform of the discontinuity of 3(5) along the cut (see figure 1) 

3' = = (-lf/'"|^J-'-M3(-J + ..)-3(-j-i.)IU, 

= (-l)M -lg^^~Hm[3{-g + iO^)] (3.14) 

Jo 

with C a counterclockwise contour around the cut (3 is assumed to be real for real g > 0). 

For large positive k this integral is dominated by the small g behavior of the discontinuity, 
where semi-classical methods are expected to be applicable. Indeed, it turns out that for 
small real negative g, the discontinuity of 3{g) is dominated by the contribution of a complex 
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Re(g) 
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FIG. 1 Contour integration in the complex coupling constant plane for the large orders asymptotics 



instanton (the real part of 3(5') is still given by the contribution of the real classical solution). 
Therefore the small g behavior of Im[3(5)] is of the form 



Im[3(-5 + 0+)] = C\g 



fQ-A\9\- 



(3.15) 



with A a number (corresponding to the instanton action Sjnst); a (positive) constant given 
by power counting (in QM and standard local field theories a = 1), C a number related 
to the determinant of the Hessian operator of fluctuations around the instanton, and /3 a 
constant related to the number mo of zero-modes of the Hessian (in standard local field 
theories [3 = 1+ mo/2). (See appendix C for details on the polymer case.) 

Given (3.15), the integral (3.14) can be calculated: Changing variables from gtox:= 
we obtain 



3fc 



air 



^^(fc-/3)/«g-Ax 



X 



1 C fi-k 

f — r 

TXOL 



a 



(3.16) 



r is Euler's gamma function. Note that since (3.15) is valid for small g only, this result is 
valid for large k. Using Stirling's formula T{n) ^ n'^e~"'^^/27rjn, this amounts to 



-1/ 



[Ae 



l3-k 



c_ 



2Txa 



[i + o(i/r)] 



and for a 



3. = {-lfk^[Aef-^k-^^^[l + o{l/k*)] 



(3.17) 



(3.18) 



19 



It is an alternating asymptotic series with a Borel transform with non-zero radius of con- 
vergence. 

B. Instanton for the SAM 

We are thus interested in the analytic structure of the partition function and the corre- 
lators of the SAM model for small negative coupling constant b 

b<0 , b-*0. (3.19) 

In particular we are interested in the discontinuity along the negative real axis. As shown 
in [14] this can be done more easily by first rescaling the fields and the size of the manifold 
with b in an adequate way. 



1. Complex rotation and rescalings for coupling constant and fields 

We consider a finite manifold M. with internal size L defined as 

Vol(A^) = . (3.20) 

We are interested in the model for small complex coupling constant b, and more precisely 
in the discontinuity of the observables along the negative real axis (6 < real), where there 
is a cut. 

We denote the argument of the coupling constant b hj 9 

e = ATg{b) . (3.21) 

To reach the cut at negative b from above or below amounts to taking the limit 

e ±Tx , \b\ fixed . (3.22) 

We now rescale the internal coordinate of the membrane x and the field r with the size L of 
the manifold and the modulus \b\ oi the coupling constant 

1 D 2-D (2-D)D 

X \b\~L~-x , r \b\2(D-e) 1 2(D-e) f (3.23) 
so that we now deal with a rescaled manifold Ms with internal size and internal volume 

size{Ms) = \b\-T^L-T^^ , Vo\{Ms) = . (3.24) 

Similarly we must rescale the auxiliary fields p and V as 

p ^ \b\-^L-^p , V \b\~i^^L-^^V . (3.25) 

The purpose of these rescalings is that as the original coupling constant b goes to 0, the 
effective theory for the auxiliary field V becomes simple. Indeed it appears that both terms 
in the effective action S[V] now scale in the same way, as will be detailed now. 
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a. Coupling constant: Let us denote by g the inverse of the volume of the rescaled manifold 

g is the (dimensionless) effective coupling constant of the theory, which is real and positive 
and goes to with |6| as long as e < D. 

h. Partition function: The original partition function (for the manifold Ai) becomes for the 
rescaled theory involving the manifold Aig 

= I V[V] exp(^-FMAV] + ^/^(O') • (3.27) 
Due to (3.26) both terms in the exponential scale as Vol{Ais) = Vfi'- 

c. Functional measure: The functional measure over the rescaled field V is now normalized 
so that 



-DiV] exp(^— jV{ryj = 1 . (3.28) 

d. Correlation functions: The moments for the gyration radius of the manifold become in 
the rescaled effective theory 

(2-D)(2d+fc) „ „ 

J n J 1-2 

= e-2'^ f f |ri-r2|'=y(ri)r(r2) 



Hence 



R'^l = e-^'' l"-^'' g"-T^' / / |ri-r2|'=K(ri)\/(r2) . (3.29) 

Jri J r2 



L is the internal extension of the original manifold Ai. This has the correct dimension L^^^*^ 
with z/q = (2 — D)/2, since [i?^'^^] = [r]'^ and [r] = Note that there is no additional 

phase for 6 = zhvr. 

2. Semiclassical limit and the effective action S[V] 
Now come the crucial points: 
1. As long as 

e < D 
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taking the semiclassical limit 6 — amounts to taking both the small coupling limit 
(yf — > in the rescaled theory and the thermodynamic limit (infinite volume) for the 
rescaled manifold 

6^0 ^ ^ ^ and Vol (Ms) oo 
for the rescaled manifold. 

2. In this thermodynamic limit the free energy F^J^] becomes proportional to the 
volume of the manifold 

Fms[V] = yo\{Ms) £[V] + ■■■ . (3.30) 
The free energy density S[V] is defined as 

in the limit where the size of the manifold is rescaled to oo, and its shape kept 
fixed. In this limit, the manifold Ms becomes locally a fiat D-dimensional Euclidean 
space M.^: 

Ms ^ '^^ . (3.32) 

The free energy density S[V] is independent of the size and of the shape of the manifold 
and it is enough to compute it for the infinite fiat manifold. 

3. Moreover - and this is an important point - as long as we are interested in the con- 
tribution of potentials V such that the manifold is "trapped" in V (namely such that 
the free energy density S[V] < is negative, i.e. such that there is a "bound state" in 
V) the neglected terms + ■ ■ ■ are expected to be exponentially small in 1/g. 

4. Finally, since 

VoUMs) = - (3.33) 
9 

in the limit g ^ the functional integral takes the standard form 

(3.34) 



Z{b) I V[V] exp 



-S[V] 
9 

where g is given by (3.26), the measure is given by (3.28) and the effective action S[V] 
for the field V is given by 

S[V] = £[V] - — / (3.35) 



S[V] is the free energy density for an infinite fiat manifold trapped in the potential V, 
and is given by (3.31). 
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3. The functional integral for negative b and the instanton 

We are interested in the imaginary part of the partition function Z(b) for b along the 
negative real axis, that is for 

9 ^ ±71 . (3.36) 
In this limit the effective action S[V] for the rescaled theory is real 

S[V] = S[V] + ljy^ (3-37) 
and the measure over V is also real, since it is normalized such that 



V[V] exp 



1 

2g 



1 . (3.38) 



It is now the standard measure for a real white noise with variance g: 

{V{n)V{r2)) = 9 6{n-r2) . (3.39) 
Thus we can chose for integration measure over V the standard measure over real V{r) 



The instanton is a non-trivial finite action extremum of the action S\V] and was 

found in [14]. The saddle-point equation is 

= ^ = V{r) + {p{r))y , (3.41) 

where p(r) is the manifold density at point r 

1 r ,D. 



and from now on we drop the index at A^^. {■ ■ ■)v refers to the expectation value for the 
phantom manifold trapped in the external potential ^(r), that is with the action 

J^d^x Q (V.r) V \/(r)) . (3.43) 

The "classical vacuum" is = (free manifold). The instanton ^"^'^t is a configuration of 
potential which is negative (potential well ^(r) < 0), spherically symmetric, with V ^ 
as |r| — > oo. The solution of the instanton equation and its properties have been studied in 
[14]. 
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C. Contribution of fluctuations around the instanton 

1. Instanton zero modes 

The Hessian matrix (second derivative of the action) is 

The instanton has d translational zero modes, corresponding to the position of the center of 
gravity ro of the instanton. Thus the Hessian has d zero modes 

V^^'" = VaV''''' ; S"[V] ■ l^f™ = . (3.45) 

According to the previous analysis, see Eq. (3.2), the metric on the instanton moduh space 
M = ds^ = hab drgdr^, is 



l-ab 



^a^y.eroy.ero ^ 5ab I d'^r ( V r"'^*) ^ (3.46) 



2ng J " 271 gd 

[using rotational invariance). Therefore the measure over the instanton position rg is 



271 g d 



1 d/2 



Hence the contribution of the instanton to the partition function will be (depending on 
whether 6 = Arg(6) = ±7r) 



instanton 

Z{b) < C± / d'^ro 



27Tdg 



1 <i/2 

2 e-^-^l^l [det'{S"[V])Y^^ (3.47) 



with C± a simple factor (usually 1 or an integer for a real instanton) giving the weight of 
the instanton in the functional integral. 

One might also expect zero-modes associated to the rotational invariance of the theory. 
Such modes would indeed appear for a non-rot at ionally invariant instanton solution. As it 
will turn out, the instanton is rotationally invariant, such no such zero-modes exist. 



2. Unstable eigenmode 

However, as expected for a theory with the wrong sign of the coupling and as shown in 
[14], the instanton has one unstable eigenmode V~{r). Thus the Hessian has one negative 
eigenvalue A~ and its determinant is real but negative: det'(S'") < 0. Therefore we expect 
that the factor C± will be complex. 

In fact, as this is the case for the instanton in the local 0^ field theory, the real part of C± 
is not unambiguously defined, but depends on the resummation procedure used to define 
the contribution of the classical saddle-point ^ = in the functional integral (this is known 
as the Stokes phenomenon). However, the instanton gives the dominant contribution to the 
imaginary part of the functional integral, and one can show that 

d/2 ^ 

e-^-^l^l \det'{S"[V])\'~^ (3.48) 



lm[Z 



instanton 

< D± 



d^ro 



27rdg 
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with the weight factor 



^2- 



(3.49) 



This result can be obtained by a more precise analysis of the respective position of the 
integration path and of the instanton solution in the space of complex potentials ^(r) G C 
as 6 is rotated from to ivr, using the steepest descent method. This is shown in Appendix 
B. 



D. Final result for the instanton contribution 

The final result for the imaginary part of the partition function at negative coupling is 



lmZ{b) 



^2 / d\ 



2Txdg 



-1 d/2 



e s 



^•'[''1 Idet' {S"[V])\ 



(3.50) 



depending on whether Arg(6) = ±7r. The infinite bulk volume factor J d'^rg disappears (as 
it should) in the normalized partition function ^ = Z/Zq 



Im3(&) 



^2^ 



d/2 



e 9 



\det' {S"[V])\- 



where C'(0) defined in (2.9). One must remember that 

g = (|6|L-)^ 



(3.51) 



(3.52) 



2-D 
2D 



and that r is in fact the dimensionless rescaled field T = r (|6|L^) 2(0-6) _ ((^i/,^) 
defined in (3.23). We thus obtain for the discontinuity of the partition function Z{b) for a 
marked manifold with a fixed point (as defined by Eq.(2.13)) 

-1 d/2 



Im Z{b) 



9 



1 



27id 



(vvy 



e-i-^'^l \det' {S"[V])\~-' . (3.53) 



For the A'^-point correlators TZ^^\ri, ■ ■ ■ ^fN'^b) defined by (2.21) the result is more compli- 
cated since the r^'s are rescaled in the process b ^ g. However he result takes a simple 
form for global quantities such as the moments of the radius of gyration of the manifold 



7^gi = (R^^l) defined by (2.28) 



>(k) 



1 

^ 2 



T+fc 



g D^"- 2D Q g 



X 



-1 d/2 



Idet' (5"[\/]) 




|ri-r2|'=y(ri)V(r2) 



(3.54) 



E. Large orders 

In the rest of this article, we shall denote for simplicity 

6 = , D = det' , Z = log(D) , W 



2Txd 



d/2 



(3.55) 
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If no UV divergences were present at e = 0, the final result at e = would be 

lmZ{b) = T^L^^^ |6|5^e"W®2n|D|^^ . (3.56) 

Using the the arguments of Sect. III. A. 2, in particular the dispersion relation (3.14) and 
(3.18), the large-order asymptotics for the perturbative expansion of Z{b) 

oo 

Z{b) = Y^Zkb'' (3.57) 



k=0 



would be (e = 0) 



Zk ~ (-l)'T U;-- —L-^''W\^\~-^G^-'' (3.58) 

y Z — L> J ZTX 

or equivalently (e = 0) 

Zk ~ {-lfv(k-2--\ —L-^^W\T)\-"^&+i-\ (3.59) 
\ 1 ) 1t{ 

indicating that the Borel transform of Z{\)) has a finite radius of convergence 6. Of course 
the instanton normalization 2n|S)|~2 depends also on d. 



IV. UV DIVERGENCES AND RENORMALIZATION 

We now discuss the UV divergences in the determinant factor for the instanton, and how 
they are renormalized. We remind the reader that at one loop in perturbation theory, for 
< e < D there is a divergence associated to the operator I (super-renormalizable case); for 
e = two divergences associated to the operators (Vr)^ and (5'^(r(x) — r(y)) (renormalizable 
case). For e < the theory is not renormalizable. The model is always considered for D <2 
and e is given by 

e = 2D--{2-D) . (4.1) 



A. Series representation of the determinant for the fluctuations 

The Hessian matrix S" is given by (3.44). We rewrite it as 

S'' = - O,,,, , 1,,,, = 5'^(ri-r2) (4.2) 



M 



where V is the instanton potential V^™*^*. O can be rewritten, using translational invariance 
X ^ X + xo when Ai M^, and the saddle point equation for the instanton potential V 



conn 



d^x(5'^(ri-r(0))5'^(r2-r(x))>; 

d^x [{6'{r, - m)6\r, - r(x))>^ - {6\r, - - r(x))>^] 

/ d^x [{5\r, - r(0))5'^(r2 - r(x))>^ - V{r,)V{r,)\ . (4.4) 
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Let us already note that such an integral is IR finite, since from clustering we expect that 
at large distances 

((5'^(ri - r(x))5"'(r2 - r(y)))™™ = C(exp(-|x - y|m)) when |x - y| ^ oo , (4.5) 

where m is the "mass gap" of the excitations for the manifold trapped in the instanton 
potential V. 

We have seen that the operator S" has d zero modes V^'^'^° oc VaV^^^, which, as discussed 
in section (III. A), are eigenvectors of O with eigenvalue Aq = 1, and one unstable eigenmode 
V~ , which is an eigenvector of O with eigenvalue A_ larger than 1. For convenience, we 
normalize its L2 norm to 1. Let us denote P° the projector on the zero-modes, and P~ the 
projector on the unstable mode 

= J2^^'%r^)V:-ir,) = , P',, = r_(ri)l^-(r2) (4.6) 

X(vv) 

and P the sum 

P = po + p- . (4.7) 

Apart from these eigenvalues, is is easy to see that all other eigenvalues of O are smaller 
than 1, but positive. Indeed, from Eq.(4.3), O is a positive operator, since for any /(r) 



f-O-f ^ 



VoliM) 



/(r(x)) 



> . (4.8) 



To compute the determinant of the fiuctuations we treat separately the negative and zero 
modes from the rest. We write the logarithm of det'[5"] 

£ = log (det' [5"]) = log(l-A_) + trp-P)log(I-O)] . (4.9) 

The first term is the contribution of the unstable mode (it has an imaginary part), the 
second term is the contribution of all other modes with < A < 1. In this last term we can 
expand the log and obtain a convergent series 

00 _ 

£ = log(l - A_) - V-Lfc , Lk = tr[(I-P)0'=] = tr [O^'j - ci - A^ , (4.10) 

k=i ^ 

provided that each term is UV finite (that is the trace is well defined). 

We now show that only the first two terms k = 1 and k = 2 are UV divergent, and require 
renormalization. 



B. UV divergences 

1. UV divergences in r and in x space 

UV divergences in the determinant are expected to come from the high momentum eigen- 
modes of S". If we consider a potential V = V^^^ + Vy, with K> a high momentum fiuc- 
tuation, we expect that a phantom manifold trapped in V will feel only weakly the small 
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wavelength variations of V, so its free energy S[V] will depend only weakly on V^. The 
other term J d'^rV'^{r) will be dominant in the variation of the effective action S[V]. As a 
consequence, high momentum eigenmodes of S" will have eigenvalues close to 1, that is will 
be eigenmodes of O with very small eigenvalues A — > 0. 

Therefore UV divergences will come from the contribution of the numerous eigenvalues 
of O close to 0, that is from the divergence of the spectral density po(A) of the operator O 
at A = 0. We shall show that po(A) diverges as \^^^~^, and that 



tr 



is UV divergent ii e < D 



tr 



UV divergent if e = 



(4.11) 



and that higher powers tr[0^] {k > 3) are UV convergent. 

The tr[-] amounts in our representation to an integral over r in bulk space M'^. UV diver- 
gences will occur as short-distance singularities in r space. We shall also see that to analyze 
the UV divergences it is more convenient to come back to the equivalent representation of 
O in X space (internal manifold). 



2. trh 



This term is given by 



tr[0] 



d^r 



and is UV divergent for e < because we expect that 



2(g--D) 
2-D 



as 







(4.12) 



(4.13) 



The crucial point (to be proven later) is that the short-distance behavior of Orr' for a manifold 
in the background potential ^(r) does not depend on the details of the potential V, and is 
given (at leading order) by that of a free manifold in a constant potential (^(r) = Vq). We 
can compute explicitly Orr' in that case and find Eq. (4.13). 

Using (4.4) we can rewrite tr[0] as an x-integral over the manifold A4, and integrate 
explicitly over r, with the result 



tr[0] 



d'^r 



d^x 

M 



d^X 
M 



(5^(r-r(xo))<5^(r-r(x))>^-y(r)^ 
)-r(xo))>^ - fv{rr^ 



(4.14) 



It contains the integral of the correlation function 

(<5'^(r(x)-r(xo)))^ 



(4.15) 



for a phantom manifold (i.e. without self-interaction) trapped in the instanton potential V{r). 
The choice of the "origin" xq is arbitrary, since (4.15) depends only on x — xq (translational 
invariance in Ai = M^). 

This integral is IR convergent, as can be seen from Eq. (4.5). UV divergences occur if 
the x-integral is divergent at short distances on the manifold, i.e. for |x — xqI 0. The 
correlation function (4.15) is very similar to the 2-point correlation function which appears 
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at first order in the perturbative expansion of the self-avoiding manifold model, and more 
precisely for the normalized partition function 3(^) 

m = 1 - ^ / (^'(Kx) - r(y))>„ + f / (5^(r(x) - r(y)) <5^(r(x') - r(/))>„ + 0{b') . 

(4.16) 

One therefore expects that the renormalization group counter-terms at leading order, which 
subtract the leading order UV-divergences in (4.16) are also sufficient to render (4.14) finite. 
That this is indeed the case will be shown below. 



3. trh 



Similarly, starting from (4.4), we can rewrite tr[0^] in term of two "replicas" of the 
manifold, labeled Aii and A^2; fiuctuating independently in the same instanton potential V 
(without interactions). If we denote ri(x) and r2(x) the r-fields for the two replicas, we have 



tr[0^ 



d'^r / d'^r' 



d-^r / d'^r' / d^x / d^y 

Ml J M2 



(5'^(r-ri(xo))5V-ri(x))>^-y(r)F(r') 
X [{5\r'-r,{yo))nr-r,{y)))^-V{r')V{r)\ 

[ d^x / d^y I (6'{ri{xo) - r^iy)) 6'{ri{x) - r2(yo)) 



(V^(ri(xo))V^(ri(x))>^ - {V{rM)V{r,iy)))^ 



1 2 



V(r) 



(4.17) 



The choice of the origins xq and yo on the two manifolds Mi and is arbitrary. 

This integral is IR finite by the same arguments as those for tr[0]. UV divergences are 
only present in the first correlation function 



5'(ri(xo)-r2(y))5'^(ri(x)-r2(yo)) 



V 



(4.18) 



very similar to the correlation function which appears at second order in 3{b), see Eq. (4.16). 
We shall see that UV divergences occur when 



X — i> xo , y — ^ yo simultaneously 
while the other terms (V(r(x)) V(r(xo)))^ are not singular. 



(4.19) 



4. trm, A; > 3: 



We can similarly write the higher order terms. At order k we need k copies Aia of the 
manifold Ai, fiuctuating in the same instanton potential ^(r). The most UV singular term 
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in the x-representation of tr[0'^] is 

[ d^x(j]5'^(r„_i(o,„i) - r,(x„)) (4.20) 

(where we identify a = with a = k), that we can represent graphically as a "necklace of k 
manifolds" . The reference points Oa on each Ada can be chosen arbitrarily, for instance fixed 
to the origin. UV divergences occur when all pairs of points (Oq, Xq) collapse simultaneously 
on each Aia- These terms are in fact UV finite for e = 0. 




FIG. 2 Diagrammatic representation of the UV divergent correlation functions at order k = 1 (one 
manifold), k = 2 {2 manifolds), and /e = 3 (3 manifolds). 



C. MOPE for manifold(s) in a background potential 

In [12, 13] the UV divergences of the self-avoiding manifold model have been analyzed 
using a Multilocal Operator Product Expansion (MOPE). This formalism was developed to 
study the correlation function of multilocal operators of the form (2.19), 

i 

where the expectation values (■ ■ ■ )q are calculated for a free manifold model (V = 0). We 
show here how this formalism can be adapted to deal with expectation values (■ ■ ■ )y for 
manifolds trapped in a non-zero background potential V{r). 



1. Normal product decomposition of the potential V 

In order to compute easily expectation values of operators in the background potential 
U, we shall use the normal product formalism already developed in [14]. 

For simplicity we consider a potential V{r) spherically symmetric (as the instanton po- 
tential) with its minimum at r = 0, of the form 

oo 

U(r) = J2^_i^^T ' ^^1 = ^0 > mo = "bare mass" . (4.22) 

n=0 
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We may (at least formally), compute expectation values of operators ■ ■ in perturbation 

2 

theory, starting from the Gaussian potential ^(r) = -^r^ and expanding in powers of the 
non-linear couplings {vk, k > 2}. This perturbation theory involves Feynman diagrams with 
massive propagators l/{p^ + '^o)- more convenient to resum all tadpole diagrams and 
to deal with an expansion of the potential ^(r) in terms of normal products : (i'^)":^- The 
normal product : [ ] :^ with the subtraction mass scale fi is defined by the global formula 
(expanded in k, it generates all operators which are local powers of r) 

.gikr.^ = e^'G,/2^kr ^ ^423) 

where is the tadpole amplitude evaluated with the propagator of mass /i, 

G - f^^_ - r n^2 

Thus we rewrite the potential V{r) given in (4.22) as 

oo 

n') -Y^i^^i'^--^ ■ (4.25) 

n=0 

The mass scale m used to define the normal product :■ ■ - '.m is defined self-consistently from 
V so that it coincides with the "renormalized mass" in (4.25) 

^1 = . (4.26) 

This gives a self-consistent equation for m in terms of V{r) (or its Fourier transform V{k)) 

All other couplings Qq, g2, gs, etc. in (4.25) are then uniquely defined from the potential V. 
We rewrite V as 

V^(r)=^o + ^:r^:+t/(r) ; t/(r) = || : (r^ : : (r^)^ : + • • • (4.28) 

and we shall treat the non-linear terms f/(r) as perturbation. The expectation value of a 
(multilocal) operator 0(xi, . . . , xk) can be expanded as 

(0(xi,...,x^))^ = ^L^ / ^d^^^(0(xl,...,x^)f/(zl)■■■^/(z^))~"^^ (4.29) 

N=0 ' 1=1 

where (■ ■ ■ ) is the expectation value in the massive free theory [U = 0). 

In this new perturbative expansion there are no tadpole diagrams. This makes the di- 
agramatics much simpler. In addition many simplifications occur in the limit d —>■ oo, as 
already noted in [14]. 
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2. MOPE in a harmonic potential 



First we consider the case of a potential quadratic in r, which is especially simple. The 
potential reads 



„2 



m 



y(r) = Vo + -^r = go + —■.x-:r,^ , mo = m ; vo = go 
The field r is still free but massive with mass m and the propagator is 



2 



G'm(x- y) 



d^p e 



ip(x-y) 



(27r)^ p2 



1 

2^ 



27r|x 



D-2 
2 



-ft"p-2 (m|x — y| 



(4.30) 



(4.31) 



where K^, is the modified Bessel Function. 

It is simple to study the short- distance limit of products of local and multilocal operators 
in this massive Gaussian theory, using exactly the same ideas and techniques as for the free 
massless case (m = 0) developed in [13]. 



a. OPE for the massive propagator Gm- We express the short-distance expansion of multilocal 
operators in terms of the expansion for the massive propagator^ 

G^(x-y) = co{D) m^^2 _ ^^p) |x-yp-^ + ci{D) jx-yp - di{D) |x-y|^"^ + ■ ■ ■ . 

(4.33) 

The coefficients Cq, Ci, do, di, are finite as long as -D < 2 and are given by 

co(l') - g§ . c.(D) - . MD) ^ . MD) ^ ^ . (4.34) 

Note that 

1 2 TT"? 

do{D) = — with Sd = , , . = volume of the unit sphere in . (4.35) 

(2 — D)Sd r(iJ'/2) 

This expansion follows itself from the OPE for the product of two r fields in the massive 
theory, which reads 

EyPl wP2 
y .yPir"vP2r'':o , (4.36) 
Pi! P2! 



^ The expansion is easily obtained from the proper-time integral representation of G,„(x), by expanding the 
integrand in to get the analytic terms in m^, and in y? to get the analytic terms in x^: 

(^i d.e-e-(P^-'^)^ = ^^ d.e-™^^.--/V^V(4.) (4.33,) 



x r X 



47r 2 IGtt 2 

1 



0(to'') + non-analytic terms in w? (4.32b) 



r ( 1 I ) |xp + 0(x'') + non-analytic terms in |xp (4.32c) 
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where the coefficient d{\x — ypm^) has an (asymptotic) series expansion in |x — ypm^ 

d{\x — y|^m^) = do + di\x — y|^m^ + d2\x — y|^m^ + ■ ■ • 

and where the normal products : ■ ■ ■ :o with respect to the zero mass means that the operators 
V*rV*r are defined through dimensional regularization (see below). 



b. MOPE for 5'^{ri — r[) and tr[0]: We first consider the short-distance expansion for the 
operator 6'^{r(x) — r(y)), which enters in tr[0]. Using the definition (4.23) for the normal 
product we can write it as 



^'(r(x) - r(y)) 



d'^k 
(27r)° 

d'^k 
(2^ 



,ik(r(x)-r(y)) 



-k2 (G™(0)-G™(x-y)) . gik(r(x)-r(y)) . 



(4.37) 



The last bilocal operator is regular at short distance (when x — * y) and can be expanded in 
X — y as 



ifc(r(x)-r(y)) . 
■ -m 



l{z) - hah (x^' - ym^" - yn : V^r'^V.r^(z) 



(4.38) 



where z = ^ and the subdominant terms are of order 0{\x — y|^) with higher derivative 
operators. We insert (4.38) into (4.37) and integrate over k to obtain 



<5'^(r(x) - r(y)) = (47r)-5 [G„(0) - G,n(x - y)]"^ x 



X 



]l(z)-^ (x^'-y^)(x--y-) .^^.^^,.y ^ 
'^'^ 4 [G^(0)-G„(x-y)]-''^''''^'^'^- + 



(4.39) 



We now use the short-distance expansion (4.33) of the massive propagator Gm(x — y) and 
insert it into (4.39) to obtain 



5'^(r(x-r(y)) = (47rdo)-^|x-y| 



-2D 



d Ci 



1 + -— ^m'^|x— y|^ — ^m^|x— y 



D 



d di 



2d, 

Sab (x^'-y^)(x'--y'-) 
4 do|x-yp-^ 



2dn 



:V^r"V.r^(z):„, + 



In (4.40) we can regroup the two terms of order |x — y|^ ^ as 
{A7rdo)-'/'\x-yr^ 



dci D^^^ 1 (x^-y^)(x'^-y'') ^ ^ .. 
--fm^]l(z) - — ^ ; ' ^ :V^rV.r(z) 



2dn 



4dn 



|x-y| 



l(z) 



(4.40) 



(4.41) 



Note that the OPE (4.40) is a relation between operators, and is valid for any choice of the 
mass m used to define the normal product. Thus the term (4.41) can be rewritten as the 
normal ordered operator : VrVr:o with subtraction mass fi = 



(47rdo)-'/'|x 



-D 



1 (x^ -y'^)(x^ -y'' 
Ado |x-yp 



:V^rV;.r(z):o 



(4.42) 
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Indeed we have the relation 

■.V^r^Vlr:^ = iV^V^rig - 2 6'^' 6,,m'' c,{D) 1 . (4.43) 

Since this relation will be crucial to prove renormalizability, let us show it explicitly. From 
the definition of the normal product we have 

-.V^^^'Vlr:^ = V^r^V.r'' - (V^r'^ V.r^)^ 1 (4.44) 

for any m, hence 

: V^r"V,r^:„- : V^r'^V.r" :o = - ((V^r«V,r^)„ - (V^r'^V.r") J 11 . (4.45) 
The r.h.s. is easily calculated using the OPE (4.33) for the propagator Gm itself, since 

{r''{xy{y))ra = 5"''G^(x-y) . (4.46) 

This yields 

(V.r'^V.r^L - (V.r'^V.r^), = ^|i7 " " ^o(>< " v)] 

= -2 5"^5^,m^ ci(D) . (4.47) 

Note that the massless propagator Go{x — y) is IR divergent but the IR divergent term is 
constant (independent of x — y) and disappears in (4.47) because of the x derivatives. Hence 
we obtain (4.43). 

Thus we have obtained the first three terms of the MOPE for the 6 operator in the U = 
background 



5'^(r(x) - r(y)) = {A7rdo{D))-^\x - y 



\t-2D 



-m |x — y| + 



4(4 - D) 

- TT {A7rd^{D)r^^^i) |x - y|^-^-2(x^ - y^)(x'^ - y^) : V^rV.rig 



t{z) 

(4.48) 



The same argument can be used to construct the higher orders of the MOPE. They involve 
higher dimensional operators of the form Op = : V^^rV^rV^^r- ■ -iq (since the operator 
(5'^(r(x) — r(y)) is invariant by translation r ^ i'+ the Op must contain only derivatives Vr, 
that is Pj > 0, and by parity in r the Op must be even in r). They give subdominant powers 
of |x - y| of the form m''^\x - y|-2i?+2/c+E, fe-i+i?/2)_ 

Finally let us stress that the two first terms of the MOPE (for D < 2) are the terms of 
order |x — y]*^"^^ and |x — y\'^~^ and that they are the same as for the MOPE for the free 
membrane, that is for m = 0. This will imply that the (one-loop) UV divergences (single 
poles at e = D and e = 0) due to this MOPE in the massive theory (self-avoiding manifold 
in a harmonic confining external potential) are canceled by the same counterterms as for 
the free theory (self-avoiding manifold with no confining potential). These counterterms are 
proportional to the operators 11 and (Vr)^. 



34 



c. MOPE for 5'^{ri — r2)5'^(ri — V2) and tr[0^]; The reader familiar with the techniques of [13] 
will see that the same arguments can be used to construct the MOPE for general products 
of local and multilocal operators in the m ^ 0, U = background. 

Let us concentrate on the MOPE for two 6 operators, which enters in tr[0^]. We are 
interested in the short- distance expansion (x ^ xq, y ^ yo) for the product of two bilocal 
operators 

5'(ri(xo) - r2{y))6\n{x) - r^) , (4.49) 

where ri and r2 belong to two independent manifolds Aii and A42- As above, we write the 
5's as a Fourier transform of an exponential and reexpress it in terms of normal products 

5^(ri(xo) - r2(y))5^(ri(x) - = I ^^^e^(^^[^^(xo)-.(y)]+^.[.(x)-.(yo)]) 



(2vr) 



d'^kid'^ks .gi(ki[ri(xo)-r2(y)]+k2[ri(x)-r2(yo)]) . g-kik2[G„(xo-x)+G„(y-yo)F(k2+k2)G„(0) 



(2vr) 



2d 



Note that there are no cross-terms, as those proportional to (^(x — y), since x and y belong to 
different manifolds, thus i'i(x) and r2(y) are uncorrelated. We now keep the dominant term 
for the OPE when x — > xq and y — > yo 

.gi{ki[ri(xo)-r2(y())]+k2[ri(x)-r2(yo)]) . _ . gi((ki+k2)[ri(xo)-r2(yo)]) . _|_ . . . 

(the neglected terms contain subdominant V^r's), rewrite this term as 

.gi((ki+k2)[ri{xo)-r2{yo)]) . ^ gi((ki+k2)[ri (xo)-r2{yo)]) g(ki+k2)2G„(0) 

and integrate over ki and V.2 to obtain 



(27r)2'^ 



gi(ki+k2)[ri(xo)-r2(yo) . gkik2[2Gm(0)-G,„(xo-x)-G,„(y-yo)] 



(4.50) 



(2vr) 



2d 



gik[ri(xo)-r2(yo)]g(fcV4-fc")[2G„(0)-G™(xo-x)-G,„(y-yo)] 



with k = ki + k2 and k' = (ki — k2)/2. The leading term is obtained by dropping the factor 
of /A in the second exponential (the neglected terms give subdominant ^'^"•'(ri — ^2) terms). 
This allows to do the integrations explicitly 

5"(ri(xo) - r2(y))5"(ri(x) - r2(yo)) - {Atx)-'"^ [2G„(0) - ^^(x - xq) - ^^(y - yo)]""/' 

X(5'^(ri(xo)-r2(yo)) . (4.51) 

From the short-distance expansion (4.33) for Gm(0) — G'm(x) the most singular term when 
both X ^ xq and y ^ yo is 

5'^(ri(xo)-r2(y))5^(ri(x)-r2(yo)) = ^ {A^d^D))^'^ '5'(ri(xo) - r2(yo)) + ■■■ • 

(4.52) 

Thus we have obtained the leading term for the MOPE in the harmonic background ?7 = 0, 
m 7^ 0. 
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This leading coefficient given by (4.52) is the same as for the free membrane {V = 0). 
The same calculation can be done for the MOPE of two 5's on the same membrane, and we 
get (at leading order) a MOPE with the same coefficient 

Pi i2— I |2 D~\ — 

5'^(r(xo) - r(y))5'^(r(x) - r(yo)) = ^ (4jj^y)J/°2 ^ ^'(K>^o) - r(yo)) + • • • • 

(4.53) 

This implies in particular that the (one-loop) UV divergence (single pole at e = 0) due to 
this MOPE in the massive theory (self-avoiding manifold in a harmonic confining external 
potential) is canceled by the same counterterm as for the free theory (self-avoiding manifold 
with no confining potential). This counterterm is proportional to the bilocal operator 6{r—r'). 



d. MOPE for higher order terms and tr[0'^]: The same analysis can be performed for the 
product of three in particular 5'^[ri(xo) — i'2(y)](5'^[r2(yo) — r3(z)]5'^[r3(zo) — ri(x)], which 
has to be considered for the quantity tr[0'^]. It shows that the leading singularity when 
X ^ xq, y ^ yo, z zq is given by the same MOPE as in the free theory, with the same 
leading coefficient. No additional UV divergences arise. The same result holds for higher 
order products of (5's. 



3. MOPE in the anharmonic potential 

We now generalize this analysis to the SAM model in an anharmonic confining potential. 



a. General discussion: The perturbative expansion involves now interaction vertices given 
by the expansion of the local potential U{r). For D < 2 and as long as no bilocal 6{r — r') 
operators are inserted this perturbation theory is UV finite. The only UV divergences that 
occur when D — >• 2 are given by the tadpole amplitudes Gm, but they are subtracted by 
the normal product prescription : ■ ■ Thus as long as the normal ordered potential V is 
finite (i.e. its coefficients gi = m^, g2, g^, etc. are UV finite) the "vacuum diagrams" are 
UV finite. Since we deal with a massive theory no IR divergences are expected. 

Now we have to consider insertions of the bilocal 6'^{r — r') operators, and thus to look 
for instance at 

N 

d^d^ n d^^^ (^'(^W - ^(y)Mz,) . . . Uizj,))'—'-' . (4.54) 

i=l 

The UV divergences which may occur when |x — y| 0, while the other distances remain 
finite, have already been analyzed with the MOPE in the harmonic case. We have seen 
there that when some Zj come close, no UV divergences occur. The only dangerous case is 
when some z's, x and y come close at the same rate. Thus we must study the short- distance 
expansion of a product of local operators (the t/'s) and of multilocal operators (the (5's), 
in the massive theory. This short-distance expansion can be studied by the same MOPE 
techniques as above. Let us ffist give a simple explicit example. 
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FIG. 3 Two contributions to tr(0) at first order in perturbation theory, associated with the 
insertion of one : r(z)^ : 



h. Example: To be explicit, we first regard as an example the simple case of the contribution 
to tr(0) given by one of the terms of (4.54) with only one U{z), and more precisely one 
quartic term : i'(z)^:. The arguments for higher powers in r or higher orders in perturbation 
theory will be identical. Following (4.14) the crucial term to calculate is 



f d^x / d^y(<5'^(r(x)-r(y)) :r(z)^:X 

M Jm 



(4.55) 



Applying Wick's theorem we can decompose it in terms of multilocal diagrams such as those 
depicted on Fig. 3. More explicitly this term can be written as 



d^x / d^y 
M Jm 



d^k 



^gik[r(x)-r{y)] .r(z)4:) 



d^x / d^y 
M Jm 

1 



d^k 



,»k[r(x)-r{y)] 



>^((k[r(x)-r(y)] r(z))^)' 



X 



[^^(x - z) - Gm{y - z)]' 



yo\{M) Jm Jm Jm J [6*^(0) - ^^(x - y)] 



i?(x,y) / [g^(x - z) - GUy - z)]' _ 

x,yeM JzeM [GmiO) — Gmi^ — V)] 

We now derive an important bound. First of all, due to the triangular inequality 

(r(y) - r{z)f < (r(x) - r(y))^ + (r(x) - r{z)f 



leading to 



Gm{0) - Gmiy - z) < 2G„(0) - G„,(x - y) - G„(x - z) 



^^(x - z) - Gm{y - z) < G„(0) - G„(x - y) 



(4.56) 

(4.57) 
(4.58) 

(4.59) 



An analog relation is valid with x and y exchanged, resulting in a bound for the absolute 
value. The r.h.s. is thus also positive and we get the bound for the ratio 



Gm(x - z) - ^^(y - z) 



G™(0)-G„(x-y) 



< 1 . 



(4.60) 
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We now want to show that the counter-terms remains the same. Using (4.60), we can write 
the bound 



M 



d^y 



< 



< 



M 



x,yGA4 



d^k 



,ik[r{x)-r(y)] 



:r(z)^:> 



5(x,y) [GUO)-GUx-y)f 

d{y^,y)[Gm{0)-GUx-y)f 



X 



G^^{x-z) -Gm{y-z) 

zeM I Gm{0) — G'm.(x — Y) 



X Vol(A<) 



(4.61) 



The latter bound is already enough to show that no additional counter-terms proportional 
to the elastic energy are necessary. It would also be sufficient for the perturbation expansion 
of tr(O^). However, we can do better and show that there is no divergence at all. To do so, 
we now estimate the integral over z. Two domains of integration have to be distinguished: 



S : 


x+y 
^ 2 


< ax 


- y| 


C : 


7 x+y 

^ 2 


> ax 


-y| 



a is chosen large (to be specified below), but finite. The integrals over z are bounded by 



Gmjy- - z) - Gm(y - ^ 
G'^(0)-G™(x-y) 



1 4 



< 



zG-S 



G^{x - z) 



^^(y 



1 4 



+ 



G'„(0)-G'„(x-y) 



zee 



1 4 



G'„(0)-G„(x-y) 



(4.62) 



Using (4.60), the first term, is bounded by 

"(^^(x-z) - G,niy- 



zes 



G'^(0)-G„(x-y) 



< 



/ l<(a|x-y|)^. 
Jzes 



In domain C, analyticity of the propagator allows the bound 



G„(x-z) 



Gm{y 



G™(0)-G™(x-y) 



< ai 



(x - y)VGm{y- 



Gm.{0) ~ Gr 



y) 



< 02 (m|x - y|) 



(4.63) 



(4.64) 



We do not give a rigorous proof here, but it is clear that a should be sufficiently larger than 
1 (say 10), which allows to establish a value for oi, itself depending on a, but saturating for 
large a. The constant a2 is chosen in order to bound VGm{x — z) by its maximal value on 
AI, which has to scale with m by power-counting in the way given above. 
We are now in a position, to put everything together. 

The integration over the distance s := x — y (which contains the possible UV-divergence) 
can now be written for small s as follows (we drop all constants for simplicity of notations) 



— s X s 2 X s 
s 



2(2-D) 



X 



for S 

S4(D-1) ^ for £ • 



(4.65) 



The factor of s comes from the integration measure; s 



-d ; 



is the leading UV-divergence 



in g^(x, y); the next factor s^*^^ is the short- distance scaling of [Grn(O) — Gm(x — y)] , and 
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the remaining factors have been estabhshed in (4.63) and (4.64) respectively. Using e = 
2D — ^-^(i, this can be rewritten as 

ds , rs2(2-^) , for 5 , 

4 for/:- ^'-''^ 

As long as D < 2, all integrals are UV-convergent in the limit of e — * 0. Thus no additional 
counter-terms are needed. The only possible UV-divergence is when first taking D 2 
before e — > 0. Note however, that this divergence only effects the contribution to the 
free energy (proportional to the counter-term 1), but cancels in all properly normalized 
observables. 



c. General analysis: We now consider the MOPE for the operator with one 6'^{r — r') and 
P = J2Pi fields r 



N 



0(x,y,z,) = 5'{r{x) - r(y)) J] :r^»(z,):^ (4.67) 

4 = 1 

when the + 2 points x, y, Zj o. The generating functional for these operators is 
1=1 ^ ^ i=i 

^\ ik(r(x)-r{y))+Eq.r(z,) _ (4.68) 



(2vr) 



d 



-k2[G„(0)-G„(x-y)] + | Y: q>qiG„(zi-z,)+iE kq^[C?m(x-zO-G„(y-z,)] 
X e * 



Expanding the normal ordered operator in x, y and z, using the short-distance expansion 
for the propagator Gm and integrating over k we get the MOPE. We see that in this MOPE 
for (4.67) local operators appear, of the form 

A = r^^V^rV^^r ■ ■ ■ V^^r , < M < P , vj > . (4.69) 

The dimension of the operator (4.67) is dim[0] = e — 2D + P{2 — D)/2, while the dimension 
of (4.69) is dim[A] = {M + Q){2 - D)/2 - R where R = ^rj > Q. Hence the coefficients 

j 

in the MOPE 

0(x,y,z,) = ^C^(x,y,z,;m)A(o) (4.70) 

A 

scale as 

C^iSx, Sy, Szf, m) ~ C° (5x, ^y, ^z,; m) + ■ ■ ■ , (4.71) 

where 

= dim[0] - dim[A] = e - 2D + {P - M){2 - D)/2 + Q{A - D)/2 + {R - Q) . (4.72) 

There will be short-distance UV divergences if the integration over the -|- 1 independent 
positions x, Zj is not convergent. This occurs if 

D{N + 1)+LJ<0 e + L)(iV_i) + (p_M)^— ^+gi^ + i?-g < . (4.73) 
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The case = has aheady been studied. When > 1, since P > M and R> Q we see that 
as long as D < 2 the condition (4.72) is satisfied only if e = 0, = 1 and P = M. Let us look 
at what this last condition means. P = M means that all the r^' in O appear in A, namely 
that no combination of propagators of the form [Y[ Gm(x — Zj) — Gmiy — Zj)] Y[ G{zj — z^) 
appear in the coefficient C of the MOPE, which therefore depends only on x — y. In other 
words, this particular coefficient comes from the product of two independent expansions 

1. The = MOPE 5'^(r(x) - r(y)) ^ \x-y\'~^^l 

2. The trivial OPE H : ^^'(zi) : ^ : r-^: with coefficient 1. 

i 

and contains no connected diagram with propagators connecting any of the Zj's to x or y. 
Thus this apparent divergence is not real. It is part of the A^ = leading divergence at e = D 
and disappears in the connected expectation value (■ ■ ■5(r(x) — r(y))?7(r(z)) ■ ■ ■ ^™nnectGd_ 

All other coefficients of the MOPE have scaling dimension u, which satisfy the inequality 
(4.72). No additional UV divergences occur beyond those which appear already for the free 
manifold and the manifold in a harmonic potential, even when e = 0. 

The same argument can be developed when there are two 6 operators and several ?7's. 
One can show that when considering the short-distance expansion of 

(5^(r(xo) - r(yo))(5'^(r(xi) - r(yi)) JJf/(x,) JJt/(yj) , xi, x^ ^ xq and yi, y. ^ yo 

i j 

(4.74) 

bifocal operators are generated by the MOPE. Power counting shows nevertheless that no 
additional UV divergence appears beyond those already studied for xi xq, yi yo while 
all other distances remain finite. 

This is sufficient to prove (at least at one loop) that the counterterms which make the 
SAM model UV finite at e = also render the SAM in a confining potential UV-finite, as 
long as D < 2. 



d. The limit D ^ 2: It is interesting to notice that there is a potentially divergent term 
when e = and D ^ 2 which corresponds to 

Ar = l, g = i? = 0, 0<M<P arbitrary . (4.75) 

(We have already seen that the case M = P is not relevant). This fact is not unrelated to 
the following observation. In the MOPE (4.48) for the single bilocal operator 6{r — r'), the 
third term, which is a subdominant term in the MOPE of the form 

5'^(r(x) - r(y)) ^ m^\x - y\'-''^+H 

is not UV divergent ii D < 2, but when D ^ 2 it becomes of the same order as the divergent 
term 

5'^(r(x) - r(y)) ^ |x-y|^-^ -.{Vrf: 

and is potentially dangerous when D ^ 2. Since this term depends on m, it depends linearly 
on the potential V^(i'), like the A^ = 1 terms that we consider here. It would be interesting 
to study this more. 
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Since Q = R = 0, this means that there are no Vr involved in the MOPE and we are 
only interested in the terms of the MOPE of the form 

5'^(r(x) - r(y)) :rP(z): ^ :r™(o): , 0<m<p. (4.76) 

It is quite easy to compute the corresponding coefficients. We find 



^-^Kx) - r(y)) :e"^W:^^ {in)-"/' [G^(0) - G^(x - y)]-^/%-^^^^M^ro^ :e-W:^ 

(4.77) 



hence at short distances 

5'^(r(x) - r(y)) : e"^^^) :^ ^ {Att 4)"''/' |x - yr'"" Q-^Hi.,y,z) . ^ar{.) ^4 ^3) 



with the function H(x,y,z) defined as 

12-D _ |y_z|2-D)2 



x-z 



|x-y|2-^ 



i/(x,y,z) = ^ ',. (4.79) 

or, after averaging with weight exp(a2j/4) 

5'^(r(x) - r(y)) :e^^'W:^ 

^ (47r c/o)"^/' |x - y 1^"'^ [1 - JdoH{x, y, z)]"''/' : g^^' W[i-^'^o^/(x,y,z)]-i ^4_gQ) 



D. Renormalization 

1. Explicit form of the UV divergences for the determinant det'(5"[F]) 

From the definition (4.14) of tr[0] as an x integral and the MOPE (4.48) for 6'^{r — r'), 
we see that the x integral (4.14) has short- distance UV divergences if e < D. The usual rule 
of dimensional regularization 

d^x Ixl"" = Sd — - — + finite terms (4.81) 
D — a 

implies that tr[0] has an UV pole at e = D, proportional to the insertion of the identity 
operator 11, i.e. 

tr[0] = Co -^—jj Wv + regular terms at e = (4.82) 
(of course (1)^^ = 1), with the residue Co given by 



d 



Co = CoiD, d) = Sd [AndoiD]]" . (4.83) 

Sd is the volume of the unit sphere in M.^ and dQ{D) = 1/(2 — D)Sd the coefficient of the 
first subleading term in the OPE of G{x); they are given in (4.34). 

Using dimensional regularization, tr[0] is analytically continued to < e < D. The next 
term in the MOPE gives the UV divergence at e = 0, hence a pole given from (4.48) by 

tr[0] = Ci (: (Vr)2 :o)y + regular terms at e = (4.84) 
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with residue 

Ci = Ci(D, d) = -Snj^^^ [ATTdoiD)]"'^' ■ (4.85) 
Similarly, tr[0^] has an UV pole at e = 0, given from (4.52) by 

tr[0^] = Ca ^ ('^'^(>'i(xo) - r2{yo)))v + regular terms at e = (4.86) 

with residue 

Here ri and r2 are associated to two independent copies Aii and Ai2 of the infinite fiat 
manifold A4. Thus we have 



5'(ri(xo) - r2(yo)))^ = / d'^ {6'{r,{xo) - r))y{S'{r,{yo) - r))y = / d'^r[(p(r)) 



where p(r) is the manifold density in bulk space. Using (4.10) and the discussion of section 
IV.B, we see that the logarithm of the determinant of the instanton fiuctuations £ = log(S)) 
has a UV pole at e = given by 



£ = log(det'[5"]) = ^ (-Ci((Vr)2>^ - | | d'^r [(p(r))^]^) 



(4.88) 



where £ms is the UV finite part of £, obtained by subtracting the UV pole of £ at e = 0; 
hence the "MS" (for minimal subtraction) subscript. 

2. Renormalized effective action 

We now study how the perturbative counterterms modify the effective action S[V] used 
in the instanton calculus. For this purpose, we now repeat for the renormalized theory the 
transformation S[r] S[V] and the rescalings performed for the bare theory in Sect. II. D 
and III.B. 



a. Renormalized original action S'ron[i']- 

The renormalized action for the SAM model is according to [13] 

^.en[r] = ^ [ (Vr(x))^ + "^^^^ 1 1 - r(y)) (4-89) 

bj- is the dimensionless renormalized coupling constant and p is the renormalization mass 
scale. At one loop the counterterms Z{bj.) and Zf,{bj.) are found to be 

Z{h) = 1 - 6,^ , Zhih) = 1+brl-- (4.90) 
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with Ci and C2 the same residues as those obtained above in (4.85) and (4.87). We first 
rewrite the renormahzed action as the bare action 5'[r] plus the "one-loop counterterm AiS'[r]. 



^ren[r] = S[r] + A,S[r] 

r Cl 1 



Ai5[r] 



e 2 



Sir] 



(Vr(x))^ 



, , (Vr(x)) 



bin' C2 
4 e 



2 , ^r/i' 



6M^) - r(y)) 

x,yeM 



(4.91) 



5'(r(x) - r(y)) 



Note that (Vr)^ =: (Vr)^:o +(i 5^(0)11 and that in dimensional regularization 5^{0) = 0. 



b. Renormahzed effective action 5rcn[^]- 

We repeat the transformation of Sect. II. D to pass from the action S[r] to the effective action 
S[V] for the effective field V{r), keeping AiS'[r] as a perturbation. We thus arrive at the 
representation for the renormahzed partition function Zi.en(&r) 



V[r] exp(-5,en[r]) 



V[r]V[V] exp 



2(Vr)^ + nr) 



J V[V] exp (^-Fm[V] + ^ Jv'^ (exp (-Ai5[ 



(4.92) 



We now perform the same rescalings and the same rotation in the complex coupling-constant 
plane as for the bare theory (see section III.B): 



\b,\fi'L^)''-'x 



r {\br\fi'L^) X , e = Arg(6,) ^ ivr . (4.93) 



Starting from a finite manifold Ai with size L (volume L^), we end up with a rescaled 
manifold Aig with volume Vol^Aig) and renormahzed effective coupling g,. 



9v 



yo\{Ms) 

The functional integral becomes 



(4.94) 



ZUbr) = J V[V] exp (^-FmAV] + |£ ^l^') (exp (-A;5[r]))^ (4.95) 



KS[r] = b. 



Cil 
e 2 



(Vr) 



4 e 



(4.96) 



0. Since 

1 



As in section III.B, 9 = Arg(6r)- We are interested in the semiclassical limit b^ 
this limit is a thermodynamic limit, where the volume of the manifold Vo\{Ais) = 9r ^ ^ 
it is natural to assume that clustering takes place (since for the instanton configuration the 
manifold is confined in the potential V). We may thus approximate the contribution of the 
counterterm by 

(exp {-/^',S[r]))y = exp (- {A\S[r])y) (4.97) 
up to terms exponentially small in Qj-. The last expectation value is 



b,Vo\{M, 



i^((Vr(o)yO,H-^| 



<5''frfo) 



(4.98) 
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Now we easily check that 

hVo\{Ms) = h/g, = e'^ {gl^'iiL)'' (4.99) 

and that when e = it reduces to e'^ = 0{1). The first expectation value in (4.98) ((Vr)^)^ 
is of order 0{1). The study of the second expectation value is slightly more subtle. We 
write 




(4.100) 



From clustering we expect that what dominates is the large-|x| regime where 

^^iKr(o)-r(x))^^ = (e^'^(°%(e-^^^W>^ = {mv {p{-^))v , (4.101) 

and where /)(k) is the Fourier transform of the manifold density p(r), see (3.42). So we finally 
obtain 

jy{r{o)-r{x))^ ^ g?Jo\{Ms) j (p(k))v (p(-k))y = jm))vf (4-102) 

also of order 0{1). (4.100) contains an UV- divergence when x and this will give a 
double pole when e — > in (4.98), but this divergence is of order 6i.Vol(7Vls)5'r — Qv This is 
in fact a two-loop divergence that we do not have to consider here. 

The final result is that we can rewrite the renormalized functional integral (at one loop) 

as 

^rcn(M = j V[V]e^^(^-JS[V]-^'g?{^JiLy^^^S[V]^ (4.103) 

with S\y] the bare effective action (3.35) and AiiS[V] the one-loop counterterm for the 
effective action 

Ai5[V] = - |i ((Vr)^>^ + I Y /(P(r))/ . (4.104) 
This amounts to state that the renormalised effective action iSrcnf^] at one loop is 

S,,^[V] = S[V]+^'gf^{^iLr'^^S[V], (4.105) 
with S\V] the original bare effective action (3.35), and AiiS[V] given by (4.104). 

3. 1-loop renormalizability 

It is now easy to show that the renormalized action for the SAM model which makes 
perturbation theory finite at one loop makes also the determinant factor for the instanton 
D = det'(5"[y''^^*]) UV finite at e = 0. 

a. Instanton contribution in the renormalized theory 

If we evaluate the renormalized functional integral around the instanton saddle point 
ymst saddle-point method, we see that the contribution at one loop of the instanton 

in the bare theory (in (3.53) and (3.54)) 



e4^[^]|det'(5"[V])r^ = e4^[^]-ii^'=(^) 



(4.106) 
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is replaced in the renormahzed theory by 

e-i'5[^]|det'(5"[y])r^e^'"*('^^)"'^i'^t^] = e~^siv]-iReis,..^) ^ ^^^^^^^ 

where the "renormahzed trace-log" of the instanton-fluctuations' determinant i2ren = 
log(Drcn) is simply (from now on we set 6 = ivr) 

£ren = - 2 (^g? fiL^ ^ Ai5[V] . (4.108) 

h. Limit e — > and UV Rniteness. 

From Eq. (4.104) for the counterterm and Eq. (4.88) which gives the UV poles of £, one 
easily checks that 2rcn is UV finite when e — 0. It is given in this limit by 

iiren = iiivis- (^^log^r + log(/iL)^ B whcn e = 0, (4.109) 

where £ms is the UV- finite part of £, as defined in Eq. (4.88), and the coefficient B is (minus) 
the residue in (4.88) 

B = Cl((Vr)2>^ + |^r(r)^ (4.110) 

(We used the instanton equation (p(r))y + V{r) = to simphfy the last term). 
Finally it is shown in Appendix E that for the instanton potential V we have 

{{vrr)y = -d{i-^y' G , Jv{rr = 2{i-^y' &, (4.111) 

where © = iS[l^] is the instanton action. Hence for e = we have 

B = (-dCi + C2)6 . (4.112) 

c. UV pole at e = D 

A similar calculation shows that the counterterm which subtracts the perturbative UV pole 
in also subtracts the leading divergence for the instanton. This justifies our use of 
dimensional regularization to deal with this divergence. 

E. Large orders for the renormalized theory 

1. Asymptotics 

From these results we can easily obtain the large-orders asymptotics for the renormalized 
theory at e = 0. The semiclassical estimate (3.56) for the discontinuity of the partition 
function Z{b) becomes for the renormalized partition function 2i.cn(^r) 

ImZ,en(&r) = T ^ i^"''' |6r | ^ e"^® 2IJ |Dren]~^ 

1 4 1 (4.113) 

= 16,1^+* (/iL)te-^®2IJ |Dms]"^ 
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with 2)ms = exp(£Ms)- The large order asymptotics for the renormahzed partition function 

oo 

ZUbr) = JZ^r&f (4.114) 



are 



2-D 2D 



k=0 



^ L~=^^(/iL)i2n |Dms]"^ (4.115) 



27r 



and the analog of (3.59) obtained by using d/2 = 4/(2 — D) — 2 at e = 0. 



2. Discussion 

From these semiclassical estimates we expect that the Borel transform of the renormahzed 
theory still has a finite radius of convergence, given by the instanton effective action &. We 
also see that as in ordinary QFT, renormalization at e = implies a dependence on the 
renormalization scale /i, an anomalous dependence on the size L of the manifold (anomalous 
dimension) and an anomalous power dependence in the renormahzed coupling constant g^. 
These anomalous dimensions are given by the factor B, which combines the perturbative 
anomalous dimensions Ci and C2 with the instanton action &. 



V. VARIATIONAL CALCULATION 



In [14] we used a Gaussian variational approximation to compute the instanton V^^r* and 
its action 5™f . Moreover we showed that the variational method was a good approximation 
for the instanton in the limit d ^ 00 (for fixed e), and the 0-th order of a systematic 1/d 
expansion. We computed explicitly the first correction in the 1/d expansion, and showed 
that for the instanton action S]^^ it was finite when e — > 0. 

We apply the same strategy here to compute the fluctuations around the instanton, 
namely the determinant factor 

I, = det'(5") = det'(ll-0) , 0..„ = -j^,gM_ (5.1) 

We first recall briefly the principle of the variational method. Then we present a direct 
calculation of D using a variational estimate for O. We show that this method does not 
treat properly the fluctuations and thus the UV divergences. We then present a calculation 
of D based on the variational method and the reorganization of the perturbative expansion 
at large d already used in [14] and in section IV. C. 



A. Variational approximation for the instanton 

We first briefly recall the variational approximation developed in [14]. We use a trial 
Gaussian Hamiltonian iftriaj''] of the form 



-f^trial[l'] 



d^X 



M 



ro)M(r - ro) 



(5.2) 
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where the variational parameters are the position of the instanton Tq and the variational 
mass matrix M = (Mat) (a symmetric real dx d matrix). The variational approximation for 
the free energy of the manifold Ai in the potential V is 

F,.,,[V] = min[F,ar[V^;M,ro]] , F^,,[V;M,ro] = Ft,,,i + {H - Ht,,,i) h,,,^, ■ (5.3) 

-Atrial = ~ 111 [/ ^[i"] exp (— iJtriaifr])] is the free energy for the trial Hamiltonian, and is a 
function of M only (translational invariance). We are interested in the limit of the infinite 
fiat manifold Ai —>■ M^, and we consider the free energy densities 

^var[l^] = ^^-j^FvarM , ^var[V^;M,ro] = ^^-jl^F.^^ (5.4) 

Obviously 

S^,,[V] = min[^,ar[V^;M,ro]] . (5.5) 
£'var[V"; M, ro] can be written in terms of the Fourier transform of the potential ^(r) 



V{p) = j d''re-'P^l^(r) . (5.6) 

and in [14] is given as 

f.„|V;M,r„l = lL|p|ltr(M°/'^) + / V'(p) e--l* , (5.7) 

where G = {G"'^) is the "variational tadpole" matrix, defined as 

G = G(M) = / 1— = LlizJllMlf-i . (5.8) 

^ ^ J (27r)^ F + M (47r)^/2 ^ ^ 

Extremization of (5.7) with respect to the variational parameters M and ro for fixed V gives 
the two equations for the the variational parameters M = M[y] and ro = ^CqIV] as a function 
of the potential V 



d-^p 



^p,V(p)eiP^«-^P«P . (5.10) 



(27r) 



Inserting these solutions in (5.7) gives £^var[^] = ^var [^^MfV], ro[y]]. Now, extremization of 
the variational effective action 

S.ar[V] = S,,,[V] + ljv' (5.11) 

with respect to variations of ^(r) leads to the equation for the variational instanton V^^^^, 

Cf(r) + (<5''(r-r(a:o)))^^^^^^ = 0. (5.12) 
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The variational instanton is rotationally invariant (as expected), so the associated mass 
matrix M^ar = M[yJ^f ] and the tadpole matrix G^ar = G(M[V;;^f ]) are constants times the 
unit matrix 1, 



r(i-f) 



(47r)^/2 



(5.13) 



(5.12) imphes that the variational instanton has Gaussian profile, and (5.9) gives Mvar as 
the solution of 

2M.ar(4vr)'^/2G.,/+l = 1. (5.14) 



The variational instanton is a Gaussian well (centered at Tq), its width is given by 



Kr(p) 



-ipro--%iJ^p2 



(27rG', 



The variational instanton action was found to be [14] 



5, 



inst 
var 



5var[K 



inst"! 
var J 



Cvar^var ( "'" ^ ^ 



(5.15) 



(5.16) 



B. A poor man's direct variational calculation of the instanton determinant D 

1. The approximation 

We have to compute the determinant of the fluctuations around the instanton solution 



V 



inst 



s\v\ 



In section V.A, we have calculated the instanton solution in the variational approximation 
Vy™/*. A flrst approximation for SD is to replace it by 



det', 



V 



d^^[\/] 
d\/[r]d1/[r'] 



(5.18) 



but this is still diflicult to compute. A further approximation is to replace this by 

d25.ar[l^] 



2). 



det',/ 



d\/[r]d1/[r'] 



(5.19) 



since we have seen that ^varl^^] for a general potential V is easy to calculate. 

This flrst and simple approximation (5.19) is presented in details in this section. We 
shall see from the result that it misses important features of the true result, especially the 
UV- divergences due to the fluctuations, which are expected as we have discussed in section 
IV. In the following section V.C, we will therefore calculate (5.18), which seems to be more 
appropriate. 
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2. Reduction to a finite dimensional determinant in variational space 
In order to calculate (5.19), we start from (5.7), and we need 

dF(r)dV(r') ^''^^ • 



We use 



Thus 



d 



dro d dm d 



dF(r) dV{r) d1/(r) 9ro dV{r)dM 



dro d£^ 



dM d£^ 



d£^ 



dV{r) dV{r) dV{r) dr^ dV{r) dM dV{r) 
since due to the saddle-point equations 

d£ d£ 







and 



. 



(5.20) 

(5.21) 
(5.22) 

(5.23) 



The second derivative is 
d^£ 



d'^£ 



dro f^^^v 



d'^£ 



dV{x)dV{x') dV(x)dV{x') dV{x')dV{x)dxo dV {x') dV {x)dm 



dro d'^£^ 



d'^£ 



dy(r') dV{x)dxo d1/(r') dV{x)dm ' 

since the explicit dependence of £^var on ^ is linear. Using the saddle-point equations 
we obtain 



d d£^, 
dl/(r) (9ro 
d 











dro S^f^v 



dM 



dV{x)dxo d1/(r) 9ro9ro d1/(r) 9ro9M 



d'^£ 



dro f^^^v 



dM d'^£^ 



dV{x) dm " dV{x)dM dV{x)dxodM dV{x)dMdM 
Eqs. (5.24) to (5.26) lead to (attention to the counter-intuitive sign) 



(5.24) 
(5.23) 

(5.25) 
(5.26) 



dV{x) dV{x') 



/ dM \ 

dV{x) 
dro 



V dVix) ) 



&^£ 



d'^£ 



(9M9M dmdxQ 
d'^£ d^£ 



\ aro^M (9ro9ro / 



dV(r') 
dro 
V dV(rO / 



(5.27) 



with (remind that everything is evaluated at the saddle-point) 

1 



£ \V 



inst jyj 
var 5 ^^^^ 



)d/2 



det(All + G)-^/V^^o^ 



(5.28) 



The quantity A is defined as follows: 



A ■= G 



inst 



(5.29) 
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i.e. it is the same as G, defined in (5.13), but always taken at the variational instanton. 
Thus when varying V, and thus M and G, only G changes, but not A. 

The determinant to be calculated is (the prime indicating that the zero-modes are omit- 
ted) 



det',/ 



det',/ 



d\/(r)d1/(r'). 



det',/ 



+ 



( 



dM 
dV{x) 
drp 
V dV{x) ) 



( 



dV(x)dV{x>)_ 
d'^E \ ( 



d'^E d'^E 



dV{x') 
drp 
V dV{x') ) 



(5.30) 



Now we use the cyclic invariance of the determinant^ to reduce the above expression (5.30), 
which is the determinant of an integral kernel operator over M"^, to the determinant of a finite 
dimensional matrix, acting on the space of the variational parameters Tq [d dimensional) and 
M (d X (i-dimensional): 



det: 



d\ 



\ 



dM 
dV{x) 
drp 
V dV{x) ) 



( 



dM 
dV{x) 
drp 

V dV{x) ) J 



&^E 



S^E 



dmdm dmdxQ 

d'^E d^E 



\ dx^dm dx^dxQ I . 



det: 







( dm 


\ 




( dM 




> 




( d'^E 


d'^E 


\ 


11- < 


r d'^p 


dV{p) 






dV{-p) 








9M9M 


dmdxQ 




J i27TY 


dro 






drp 








d'^E 


d'^E 








V dV{p) 


I 




I dV{-p) 


I 


> 




\ dx^dm 


dx^dxQ 


) 



(5.31) 



1 is the corresponding d[d + l)-dimensional unit-matrix. In fact the variational mass ma- 
trix parameter space is d[d + l)/2 dimensional, since one has to consider only symmetric 
mass matrices M. However in our calculation is is simpler to consider the (i^-dimensional 
variational space of all real matrices M. 



3. The calculation 

We now evaluate the elements of the matrix. First of all, due to rotational invariance 
and parity of the instanton, the off-diagonal blocks of the two matrices {□} and (□) vanish 



d'^E 



dx^dm 
d'^p dm dxQ 
{2TxYdV{p) dV{-p) 



(5.32) 
. (5.33) 



^ If X is a n X TO matrix and Y a m x n matrix, and det' denotes the product over non-zero eigenvalues, we 
have the general identity det'[l — XY] = det'[l — YX], although the first determinant is the determinant 
of a n X n matrix, and the second one the determinant of a to x m matrix. 
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The second relation will be explicitly checked below. As a consequence (5.31) takes block- 
diagonal form, leading to the factorization of the determinant as the product of the deter- 
minants over each diagonal block 



(1) 



Svar' = (5.34) 

aet'll- f^^^^P^] (5.35) 

[ Jp dV{p) dV{-p) dMdM\ ^ ' 

det'h- [ 4^^^^^] . (5.36) 

Second, we shall see that the second block, relative to the zero-mode collective coordinate 
ro, is also 0. Indeed, we shall show that 

dro _ dro S^i^var 



(2) 



dV^(p) dV{-p) drodro 



1 



so that 



(2) 



det' [0] 



1 . 



(5.37) 



(5.38) 



Thus it remains to compute the determinant of the first block, involving only dependencies on 
the variational mass M. Using (5.28) and the matrix derivative rules gathered in Appendix 
D, we find 

_^!^ = [2(2 + D)E - d(2 - Dm , (5.39) 

with E the projector on symmetric matrices and P the projector on the unity matrices 



^ij,kl = 2 i^ik^jl + ^il^jk) 



ij,kl 



(5.40) 

(5.41) 



Next, we calculate ^y^- Using Eq. (5.9) and varying V yields 

5W^[V] = - 



^"^■^ p'p^ 5V{p)e'P''>e-^p'p'^^^ 



1 

+2 



d'^p 



(27r)' 



p'p^ V{p) e'P^° e-^P'P''^'^ p'^p^SGki 



(5.42) 



2 M 



Using that at the saddle-point 6G — ^ ^ ^ 
d'^p 



6W^ [V] 



{2nY 
d'^p 



p'p^SV{p) e*P^°e"^P'^-^ 



and V{p) from Eq. (5.15), we obtain 
^pp^e P P^-^ 



p'p'6V{p) e*P^Oe-^P'^ + ^-— ^ 



(27r 

6Mm (rfP,,, + 2%,fcz 



kl 



This leads to 



- — —dF5M - ^^-^5M 



(2vr) 



^p*p^<5K(p) e*P^°e 



(5.43) 



(5.44) 
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and finally upon varying 5V 

_^(i_^,P_i±^E)^p«pe-e-l-=^. (5,45) 

This can be inverted (in the subspace of symmetric matrices) as 

( —E - J^dC^-D) \ ,p,o lyA ^5_4g) 

V 2 + D {2 + D){A-2d + 2D + Dd) )^^^ ^ ' 



5V(p) 



Next, we need ^^[^ ® trnfm ■ saddle-point equations, or more explicitly 



looking at Eqs. (5.45) and (5.39), the following combination is relatively simple: 



5M d'^S^^r A2-D 



p® pe*P^«e"2P ^ , (5.47) 



6V{p) dMdM M 4 
and after (Gaussian) integration over p we obtain finally 

f d'^p 6M SM 92^varo.„o.„ 2-D^ 2d(2-D) _ , , , 

/ ® ^(9M9M = E + ^- -P . 5.48 

J {27tYsV{-p) 6V{p) 2 + D {2 + D){e + 2-D) ^ ' 

The first block determinant (5.35) is therefore the determinant of the following operator 
acting on the d{d + l)/2 dimensional space of d x d symmetric matrices 

- det' (^E ^^(^ - ^) A (5 49) 

Since in this space the projector E reduces to the identity, while P is the projector on the 
1-dimensional subspace generated by the identity, it is easy to see that the operator has 
d{d + l)/2 — 1 eigenvalues equal to 2D/ (2 + D), plus one eigenvalue equal to 2D /{2 + D) — 
2d(2 - D)/{2 + D){e + 2-D) = -2{D - e)/(e + 2 - D). Hence the final resuh is 



d(d+l) 



'—1 

e + 2-D \2 + D 



- ^..^ - ) . (5.50) 



4. Terms associated with the zero modes 

Before discussing this result, we calculate the other entries of the matrix (5.27), associated 
with the 0-modes. First we vary Eq. (5.10) with respect to SfQ and the corresponding SV{p): 

I (0j^^(P)^Pe-^^''VP^'' = / (0^^(P) pe-^P^VP^nP^ro) . (5.51) 
Deriving with respect to SV{p) and evaluating at V^nst yields 

■I j T j 

= -Mi.-^ = -M6i,-^^ . (5.52) 
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This gives 



d\/(p) 



M ^ 



Combining Eqs. (5.46) and (5.53) checks (5.33). 

We now calculate the determinant of the lower block, for which we need 



(5.53) 



d-^p dr^ 



(27r)^dV(p)dV(-p) J (27r) 



^p^p^e-P^ 



M 



ik 



as well as 



dr^dcQ 



(27r) 



F(p)e*P^0e-5P-'^-P 



d'^p 

(27r)^ 



p p-'e ^ 



(5.54) 



(5.55) 



where we used that the first term of S^ar in (5.7) does not depend on ro, as well as the 
instanton at the saddle-point from Eq. (5.15) and the mass from Eq. (5.14). Hence the 
second block matrix, relative to the zero mode rg, is identically zero. This is not surprising. 
Therefore 



d-^p dr^ 



dr^- 



(27r)'^d1/(p)d1/(-p) J 9r*9r^o 



8 



ik 



(5.56) 



and indeed the determinant (5.31) is the contribution of the d translational instanton zero- 
modes. 



(2) 



det; [0] 



(5.57) 



5. Discussion 

We now discuss our result (5.50) for Svar' in our simple variational approximation. We 
see that D^ar' is finite and negative for e < D, thus we recover the unstable mode with a 
negative eigenvalue for S". However we see that for e = 0, D^ar' is still finite, while we expect 
from our general argument that D will have UV divergences. Thus our approximation does 
not properly take into account the short- wavelength fluctuations around the instanton, and 
renormalization, which is important when e — > 0. 

Finally it is interesting to look at the behavior of D^ar' in the limit d ^ oo, e fixed. We 
find for the logarithm of Dvar'? 

£w = log(S)w) - (l - ^) = ^(^) (5-58) 

as expected from the variational approximation. However, as we shall see later, the better 
approximation £var and the exact solution £ behaves respectively at large d as 



(5.59) 
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C. Expansion around the variational approximation and \/d expansion 

1. The large-d limit 

A better method to compute D is to start from (4.3) 

O,, = / {^r, - r{o))6\r, - r(x)))>;r. (5.60) 

J X 

(o is an arbitrary point on = MP) and to make a perturbation expansion around the 
variational Gaussian Hamiltonian -fftrial- Since the problem is invariant under translations, 
we chose for V the instanton centered at the origin (tq = 0). m will denote the variational 
mass (M = w?) and Gm the variational tadpole = (47r)"^/^r((2 — D)/2)im?~^ . m is 
solution of (5.14), that we rewrite 

2m^Gm = {^nGm)-"^^ • (5.61) 
a. Large-d limit and the variational approximation: 

The first crucial point used in [14] is that when the variational instanton potential (5.15) is 
written in terms of normal products relative to the variational mass m, it takes the simple 
form 

Kaf (0 = -{47rG^)-^/^:e-^:^ = ~2m^ G^-.e-^ -.m (5.62) 

that we rewrite as the variational trial potential |m^r^ plus a perturbation U (r) as in section 
IV.C (see (4.28)) 

Krf(r) = -2m^G^]l+^:r^:„ + f/(r) , Uir) = -2m^G^J2^\-^) ^(^')"^. 

(5.63) 

and to treat U{r) as perturbation, see (4.29) and Fig. 4. 



V 













FIG. 4 The variational instanton (black) and its approximation by a harmonic potential (grey) 
(here for D = 1, d = 4). Note that the curvature of V{r) is the quadratic term before normal- 
ordering, whereas in the variational approximation the quadratic term after normal-ordering ap- 
pears. 

The second point is that in the limit when 



c? — s> oo , e fixed 



(5.64) 
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these perturbative terms are subdominant (of order 1/d) with respect to the leading term 
obtained by replacing V by the trial harmonic potential ^ : r^:^- This is seen by rescaling 
X and r in units of the variational mass m (as described in detail in Appendix G), so that 
m — >• 1, and the propagator Gm(x) becomes G(x) = Gi{x) 

Gra{x) ^ G(x) = G„=i(x) = (27r)"f K£_2(|x|) (5.65) 

and the tadpole amplitude Gm becomes 

G^ co{D) = G(0) = (47r)-^/2r ( ) ~ when d ^ oo , e fixed . 

\ 2 J 47r 4 — e 

(5.66) 

Co(-D) is noted C in [14]. When we shall not deal with the explicit dependence on D of cq{D) 
we shall denote it simply by cq. 

The variational instanton potential becomes (see Appendix G) 

V^\r) = -2col+ --.r'-. + Uir) , f/(r) = - 2co g ^ (^— j :(r2)":, (5.67) 

where the normal product : ■ ■ ■ : refers to the normal product with respect to the unit mass 
m = 1, i.e. :■■■:=:■■■ 

Since Cq ~ d, in perturbation theory, the 2?T,-leg vertices carry a weight rf^"" and closed 
loops carry a weight d (summation over bulk space indices). Counting the resulting factors 
of d for each graph, as in the large- iV expansion for vector models, only "cactus diagrams" 
with tadpoles survive in the large-d limit. However within our normal product scheme, there 
are no tadpoles. Therefore for any observable at large d we can replace 

(Observable) yinst = (Observable)^ + subdominant terms in - , (5.68) 

where {■ ■ refers to the expectation value with respect to the trial variational action 

= + (5-69) 

For the same reason, as shown in [14] , at leading order in 1/d, the variational instanton, 
solution of 

Kaf(0 + (5'^(r-r(o))>^ = (5.70) 
is a good approximation for the exact instanton V^'^^^, solution of (3.41): 

r-*(r) = Kr(r)(l + 0(l/ci)) . (5.71) 

The first correction of order 1/d was computed in [14]. Finally the action for the variational 
instanton was found to be 

SZ' = (l - -i) coiD) . (5.72) 

If we rescale the effective coupling constant g (or equivalently the initial coupling constant b) 
in terms of the variational mass m, 

g rrP g i.e. h m^'" b (5.73) 
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the instanton action becomes 

e 



SZ' = [^-^) Cop) = 0{d) . (5.74) 

This rescahng will the done at the end, but for the time-being, we keep the explicit mass 
dependence. 

b. Large-d limit for O; 

For our problem, in the large-d limit, we shall firstly approximate the Hessian O in the 
exact instanton background, with kernel Orir2 given by (5.60), by the Hessian O^'^'' in the 
variational instanton background, with kernel O™,!^ given by 



jvar 
rir2 



{6\r^-r{o))S'{r,-r{x))y;:i (5.75) 

^ var 

and then approximate this O™ by its large-c? limit O™' , with kernel 

[6'{r, - r{o))6'{r, - r{x)))2'''' • (5.76) 



rir2 



This will be the leading term of a systematic ^ expansion, which can be performed along 
similar lines as in [14]. 

0™r2 can easily be computed, since we now deal with a massive free theory. It is even 
easier to compute its Fourier transform 

0™J^'_^ = g-i(kiri+k2r2) Qvar _ / ^^ikj r(o)gik2r(x) \^ _^gikir(o)^ {^C^'^^''^'*-' ) 



, , (5.77) 

/ ^g-kik2G,„(x) _ 
J X 



= g-(k2+k2)G„(0)/2 

where Gm(x) is the massive scalar propagator (4.31). Note that we have Gm = G'm(O). 



c. Zero modes: 

In order to compute D, we must take into account the translational zero modes of S"= l — O 
and the projector Pq onto the subspace of zero modes. According to section III.C, these zero 
modes are the partial derivatives of V^^^, V^^"^" = daV^^^, and from section IV. A (see (4.6)) 
the projector is 



a 

I>2 



(with the constant Cq such that Pq = Pq). In the large-ci limit we may approximate Pq by 

p^var 

Po™^. = c'o5^9,y4f(r09„K;r(r2) (5.78) 

a 

and since V^ar* is a Gaussian function, P™'' is easily computed. We obtain for its Fourier 
transform 

Po^l = -^e-^(^?+^i)«'"(r (5.79) 
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Finally, in the large-order formulas such as (3.53), to the instanton zero-modes is associated 
the weight factor 



2IJ = g~T^ 



2Txd 



inst \ 2 



d/2 



In the \aige-d limit this gives 



2n™'- = g-D 



2Txd Jf 



(VK 



-1 d/2 



inst N 2 
var / 





r 21 


d/2 


2 


d 


m 


9 ^ 
2tx 


D 


_2ti_ 





d/2 



(5.80) 



so that 



log (2n™'') = 0{d) 



2. Large-d calculation of £ 

a. Series representation for £; 

We now apply these results to the computation of the determinant, or rather of its logarithm 

£ = log(det'[5"]) . (5.81) 

Since S" = 1 — O and since det' subtracts the zero modes, we can write 

£ = tr[log(Il-Q)] ; Q = O - P„ . (5.82) 

Note that £ has an imaginary part since S" has one negative eigenvalue A~ . We expand the 
log as 

£=-^-tr[Q^] . (5.83) 

k=l 

As we shall see, further simplifications occur in the large-ci limit. In this limit we can 
approximate Q by Q™'' given by 

where O^^''' defined by (5.76) is the Hessian —£" at the variational instanton, computed in the 
variational approximation, while F^^^^ defined by (5.79) is the projector on the zero-modes 
of S" in the variational approximation. Therefore we approximate £ by 

(5.85) 

fc=i " 



■tr 



b. "Beads" and "necklace" diagrammatic representation: 

Starting from (5.77), (5.79) and using the fact that J^Gm{'>^) = l/rri^, we can write the 
kernel of Q™' as 



Km = e 



-(k2+k2)G„{0)/2 



f [e-k,k,G„(x) _ ^ ^ kik2G„(x)] , (5.86) 

J X 
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and expanding in kik2, we get a simple diagrammatic representation for Ql^^l^ ^ snm of 
"watermelon" diagrams 



n=2 



(5.87) 




Each line represents a propagator G^- No internal Ai momentum flows in the diagram, the 
p's are external momenta relative to the embedding space M.'^. In this series the term n = 
is removed by the fact that O is defined by a connected correlator in (5.60); while the term 
= 1 is removed by the projector onto the zero modes Pq- 
Now we consider the trf 



tr 



] in (5.85). Each trace is given by 
{2ttY (27r) 



d ^ki,-k2Vk2,-k3 



var 

kfc,-ki 



(5.^ 



Thus £™''' can be represented as a sum over "necklace" diagrams made out of the "beads" 
of (5.87). The integration over the k's can be done explicitly and gives a decomposition of 
the form 



tr 



i=l 



ni,...,nfe>2 ■ " '■■ 



, (5.89) 



where Pni,...,nkid) is a polynomial in d (the bulk space dimension), with integer coefficients, 
given by the average 



Pnu...,nM = (-klk2)"i(-k2k3)"---(-k,,ki)«^- 



(5.90) 



with the normalized Gaussian independent variables kj G M.'^, i.e. k^-k^ = 6°'''6ij . The 
polynomial P can be computed by Wick's theorem. Typical configurations are: 





Note that the first and last points are identified. Let us denote by the total number of 
lines = ^rij in the diagram. From (5.90) the P's are non zero if and only if the n^'s are 
either all even, or all odd. 

• If A; = 1 this is always true and Pn{d) is of degree n in d. 

• If A; > 1 and the n^'s are even, N > 2k and the degree of P{d) is N/2 > k. 

• If A; > 1 the rij's are odd, N >3k and the degree of P{d) is 1 + (A^ - A;)/2 > k. 
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c. Large-d power counting: 

We now look at the behavior of these terms when c? — >■ cx), e fixed. First we rescale everything 
in units of the variational mass m, 



X — i> x/m 



2-D 

p — i> pm 2 



(5.91) 



that is set the variational mass m to unity in our calculations, since the tr[Q''"] are dimen- 
sionless quantities. We refer to Appendix G for the details on this rescaling. Then we note 
that the propagator G{x) for x 7^ given by (4.31) is of order 0{1) when d 00 

G(x) = (27r)"T|x|^i^^(|x|) ^ :^A^o(|x|) = 0{1) , 



while G{0) is of order 0{d) since 

G'(0) = (47r)-fr(2-i) 



1 



d 



2tt2-D 47r(4-e) 
Thus the integrals /„ given by (5.89) are of order d"" 

= o{d-n , 



Oid) . 





rG(x)i 


L 


[g(o)J 



and the term associated to the /c-bead necklace [?t,i,?7,2, ■ ■ ■ , ri^] in the decomposition (5.89) 
is of order 

where N = Y^Ui. 

• If /c = 1, we have seen that degree [P] = A^, and all the terms are of order d. Therefore, 
if the series over the n's converges (we shall discuss this later) 



tr 



0{d). 



(5.92) 



If A; > 1 we have seen that there are two cases. For even necklaces the n^'s are all 
even and degree [P] = N/2 > k so we obtain a term of order d ~~ < d^ We note that 
the most dominant terms are those with = 2k. These are the [2, 2, ■ ■ ■ , 2] necklaces 
whose beads contain 2 links (chains of bubbles). 



(5.93) 

For odd necklaces, the nj's are all odd and degree[P] = 1 + (A^ — k)/2 > k, while 
A^ > 3A;. This gives a term of order d^ ~ < d^~^ <^ 1. The conclusion is that (as 
long as we can sum the necklace series) the k > 1 terms are of order 0{1) 



k > 1 



tr 



'\ k 



) 



0(1) 



(5.94) 



and that the dominant contribution is given by the chain of bubbles. 



59 



d. Final result: 

All the k > 1 terms in (5.89) are subdominant with respect to the k = 1 term. In the large-c? 
limit, £ is of order 0{d) and can be approximated by 



tr 



Oil). 



(5.95) 



We shall check this result with explicit calculations. As we shall see, the summation of 
the necklace series is not completely obvious, and is impaired by the UV divergences of the 
theory. Let us also note that the imaginary part which comes from the unstable eigenmode 
of 5" = 11 — O is an effect of order 0{1) (since it is associated with one single eigenvalue). 



D. Explicit calculations at large d 



1. tr 



and its large-d limit for e > 0: 



We first consider the leading term tr [Q™'], given by 



tr 



tr 



d 



(5.96) 



a. x-integral representation: tr [O™'] is easily calculated from (5.77) and (5. 



tr 



-k2G™(0) 



,k2G„(x) 



{27tY J {27rY 

The k-integration is Gaussian and gives, using the equation for m (5.61) 



tr 



1 



^^(0) 



(5.97) 



Since tr [O™'''] is dimensionless we can set the variational mass m to unity m = 1, in the 
r.h.s. of (5.97) (see Appendix G). Using the explicit form (5.65) for the propagator G{x) 
and integrating over the x angular variables via f d^x = So dxx^~^ with x = \x\ we 
obtain 



tr 



-,2-D 



r(f) 



dxx 



D-l 



1 - 2 



X 

l2 



2z_D Kd-2(x) 



(5.98) 



Let us first consider this integral for finite (but a priori large) d, and study its convergence. 



b. IR convergence: 

At large x the integral is convergent. Indeed the massive propagator is exponentially de- 
creasing as G{x) ^ exp(— |x|). For finite d, and thanks to the —1 that comes from the 
subtraction of the disconnected part, the integrand in (5.97) is also exponentially decreasing 
at large x. 
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c. UV divergences: 

The small-x behavior of the integral (5.97) has in fact aheady been studied in Sect.(IV.C). 
It was shown that this behavior is governed by the MOPE (4.48) and is related to the UV 
divergences at one loop of the model. The integrand in (5.98) behaves as dxx'^^^^^Co 
with Co given by (4.83). We thus recover the expected UV divergence at e < D, which is 
proportional to the insertion of the operator 11. This UV divergence appears in the series 
representation (5.89) of tr [O™'' ] as the onset of the nonsumability of the series^. Indeed, 
this series is 



tr 



n=l 



2G(0)^:^^4 with Pn = d{d + 2) ■■■{(! + 2n -2) 



(5.99) 



and /„ = JJ^G{x)/G{0)]"' ~ n ^-d at large n. It is easy to check that the series (5.99) 

behaves as Xln'^"^^^"^ convergent only if e > -D. 

Since the model is defined for e < D by dimensional regularization, the analytic continu- 
ation of the integral (5.97) is its finite part (in the sense of distribution theory). Therefore 
tr [O^'^'^'] is defined for e > by 



tr 



2G(0) X f.p. 



1 - 



G(0), 



- 1 



(5.100) 



or equivalently by the resummation of the series (5.99) by a zeta-function prescription. 

For e = the integral has another UV divergence, which is canceled by the (Vr)^ coun- 
terterm of the renormalized theory. We shall discuss this point later. 



d. large-d limit: We can now take the limit of (5.99) when 

d ^ oo , e > fixed - 

Since in this limit 



G(x) ^ ^i^o(|x|) 



G{0) 



1 d 
4:71 A - e 



we obtain 



tr 



d- — ^ f.p. / dxx [e(^-'^^oix) _ ^ 0(1) 



(5.101) 



From the short distance behavior of the 2-dimensional propagator Kq{x) ~ log(l/x), the 
last integral is 

roc poo 

Ti(e) = f.p. dxx [e(^-^)^°(") - l] = / dxx [e^^-^)^"^") - 1 - x'^+^j 
Jo Jo 

and is UV finite for < e < 2. Thus we recover that tr [O™'] = 0{d) in this case. 

Ti(e) has a single pole at e = 2, as expected. It is UV divergent when e — ^ 0. This will 
be studied later. 



^ This nonsumability has of course nothing to do with the large-order behavior we are after. 
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2. tr 



and its large-d limit for e > 0: 



We now perform the same analysis for tr [Q^]. We have, using (5.61^ 



tr 



(2^(2^ 



-(k?+k2)G„{0) 



(5.102) 

Setting 771 = 1 and performing the k integrations we get 



d 



G(xi)+G(x2 ) 
2G(0) 

G(xi)G(x2) 



- 2" 


d 
2 


1 - 


rG(xi)i 


2" 


d 
2 


1 - 


rG(x2)i 


2" 






.2G(0)_ 








.2G'(0)_ 





1 - 



g(xi) 
2^(0) 



1 - 



2G{0) 



(5. 



03) 



This integral is IR and UV finite as long as e > 0. When e = we recover the UV divergence 
when both xi and X2 0. 

Now in the large-ci limit, e fixed, since G{0) ~ d and G{x) ~ 1 we can expand the [■ • ■]^'^^^ 
and get 



3d(d + 2) 
16G(0)2 



o{d-'] 



(5.104) 



This expansion is not valid for xi and X2 = when e = and this gives the 1/e UV pole 
(coupling constant renormalization) , but this is an effect exponentially small in the large-d 
limit. We have thus checked the fact that 



tr 



Oil) 



(5.105) 



3. tr 



and its large-d limit for e > 0: 



Calculation of higher powers can be done along the same line. We get 

k' 



tr 



f \ k 

1--) +oid-'). 



(5.106) 



VI. 1/d CORRECTIONS TO THE LARGE d LIMIT 



In this section we study the first 1/d correction to the variational solution, which was 
shown to be valid for large d, e being kept fixed. 
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A. 1/d diagrammatic 

We first recall in this subsection how is constructed and organized the \/d expansion, 
following the ideas of our first paper [14]. We have performed the rescaling (5.91) so that 
the variational mass m is set to unity. This rescaling is detailled in Appendix G. We denote 
by Co the normalized tadpole amplitude^ and the integration measure over ci-momenta k is 
now normalized so that we have 




1 , co = (47r)-^/^r((2-D)/2)=r ) , /e-^''=°=2co (6.1) 



m 

The exact instanton potential is in these units of the form 

°° 1 /-I \ " 

^'"■V) = 2c„EF;j(5^j (6-2) 

n=0 ^ ^' 

where now the normal products are defined with respect to the unit variational mass m = 1, 
i.e. 

• • •m=l 

and the coefficients of order 1 in the large-d limit, and are found to be 

/Xn = -1 + ^ + 0{d-^) (6.3) 

in the large-ci limit, with (5„ = 6n{D,d) given by a self-consistent equation that we recall 
later. We remind the reader that if we set the //„ = — 1 we recover the variational instanton 

p^inst 



order n 



2n — vertex 



V 



couplings co{d/2 — 2^o) 
0{l/d) co{d/2 + 2-25i/d) 



ii. 

d 



-IJ-2 


M3 


-At4 


2co 




(2co)^ 


1 


-1 


-1 


2co 


(2co)^ 





FIG. 5 Self energy (n = 0), mass (n = 1) and interaction (n > 2) vertices and couplings in the U 
expansion (the symmetry factors l/(2"n!) for the vertices are not written). 



The perturbative diagrammatics is obtained by writing 

r-*(r) = lv^ + U{r) (6.4) 



^ Co is denoted C in [14]. 
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and treating [/ as a perturbation. The corresponding 2n-vertices and couplings are schemat- 
ically depicted on Fig. 5. The last line represents the couplings which have to be kept at 
order 1/d. The propagator is the usual bosonic propagator with unit mass G{x). The one- 
loop tadpole graph is absent since it is subtracted by the normal-product prescription. The 
external r-space indices a = 1, ■ ■ ■ , ci flow along the closed lines as in a standard 0(n) model. 

FIG. 6 The chains of bubbles in the large-d expansion 



It was shown in [14] that the diagrams can be reorganized in a l/rf expansion by summing 
all the chains of bubbles, as depicted in Fig. 6. More precisely, the propagator for the chain 
is given by the geometric series (in Fourier transform w.r.t. x space) 



l+/X2/-5(p) 



n -1 



(6.5) 



p being the D-momentum flowing through the chain, and B{p) the bubble amplitude (one- 
loop diagram) 



= ^(P)= l^'^'Gixf 
Jp 

^ 1 arcth {p/^/4 + f') ^ i arcsinh {p/2) 



(6.6) 



For zero momentum, we have 



vr 



when D = 2 



2 — D 



(6.7) 



In practice we also have to consider the chains with n > 1 ot n > 2 bubbles. They are 
depicted as follows, with the associated amplitude H^^\p) and H^'^\p) 



^---o+oo+ooo 



l+,.-i^(p) 



(6.8) 



-H- 



oo+ooo 



4co 



Acq 



(6.9) 



At the diagrammatic level this reorganisation of perturbation theory is very similar to what 
is done in the expansion for the (linear or non-linear) sigma models, where the bubble 
chain is the propagator for an auxiliary a fleld, and the interaction involves only rra and 
cr'^ {k > 3) terms. The analytic structure of the perturbation theory is nevertheless quite 
different, in particular for the UV and IR divergences of the theory, as already discussed in 
[14], and as we shall see below. 
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vertex 



coupling 



"M2 
2co 



1 

/2^(-/^2)3/2 



1 
2co 



coupling at order ^ 



1 



1 
2co 



FIG. 7 The vertices contributing to —y"^^^(^r) and their couplings in the large-d-reorganized per- 
turbative expansion. 



After this resummation the new vertices with their couplings are depicted in Fig. 7. The 
crucial point is that in the limit d ^ oo, e fixed, since D — > 2 the tadpole coefficient cq 
diverges as d so that 

/i27^^ > -8n(l-^] = 0{1) (6.10) 



■2co 



4 



and the bubble propagator H(p) is of order 0{1), while the vertices are of order 1/^/d, 1/d, 
etc. It was shown in [14] that only a finite number of diagrams contribute to a given order 
in 1/d, and explicit calculations where done at the first non-trivial order. 

With these notations we have found in [14] that at order 1/d the following diagrams 
contribute to the expectation value of the exponential (or vertex) operator 



CO 




(k^ 



(6.11) 

The symmetry factors of the diagrams are not written, they are respectively 1/2, 1/4, 1/2, 
1/4, 1/4 and 1/8 for the diagrams in (6.11). No r-space indices fiow through the unclosed 
line^ diagram is (taking into account the couplings and the symmetry factors) 




^(l + /ii) /"g(x)2= 1—^(1 + ^i)co 



This is different from the foUowing graph considered in [14], whose ampUtude differs by a factor of d d. 
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Similarly for the last diagram 



1 -1^2 

8 2cn 



GixfGiyfHix-y) 



The exact instanton saddle-point equation, which is (once again) 

V{k) + (e^'^(°))^ = 
fixes the /in's. In particular /ii is given by 

d L 



(6.12) 



(6.13) 



and using (6.12) and (6.11) at order l/rf we get the equation for /ii (i.e. 5i) which reads 
diagrammatically 



+ 



2_co 
d 



( 





V 






-A 





k6 



^0 



, (6.14) 



where and mean the average value of k^ and k^ respectively with the Gaussian weight 
e~'^ Since 



k4 



d{d + 2) ^ f d 
(2co)2 V2q5 



k6 



d{d + 2){d + 4:) [ d 



(2co)= 



2cn 



(6.15) 



are of order 0{1) we recover that /ii = — 1 + 0{l/d). 



B. The Hessian O 

We now show how this method to construct a l/rf expansion can be applied to compute 
the matrix elements of the Hessian S" and of the associated operator O. We start from the 
expression for O in momentum space 

6,,k, = / (e^^^^(°)e^^^^W>7" = I (e^^^^(°)e'^^^W>^ - (e'^^^(°)>^ (e''^^^^))^, (6.16) 

>/ X «/x 

and we use our perturbative rules to expand the e.v. (■ ■ ■ )y in l/d. 



1. O at order 1: 

At leading order C(l), we get (5.77) that we can represent as a sum over diagrams with 
n > 1 propagators between o and x, integrated over x 




(6.17) 
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where the integration over x and the symmetry factor 1/n! of the graphs are imphcit. We 
have introduced here an additional diagrammatic notation, which will be very convenient in 
the following discussion. 

2. A diagrammatic representation for the vertex V^(k): 

The circles in the last graph are a symbol for the factors which depend respectively on 
ki and k2 and are attached to the vertices o and x. More precisely, the circle represents the 
exponential and each line entering into the circle represents an additional (multi- 

plicative) factor ik, with an external space index a carried by the line. Thus the following 
picture, a circle with n external lines, represents the factor 



1 




' n 



3. O at order l/d: 

Now we make the perturbative expansion and keep the diagrams which contribute to O 
at order 0{l/d) only. We find that only (!) 21 different (classes of) diagrams contribute 




(6.19) 
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4. Diagrammatic for the mass renormalization and V: 

Moreover, Eq. (6.14) for the mass renormalization /ii may be rewritten at order 0{l/d) 
with our notations as 



(6.20) 



While Eq.(6.11) for V reads 



-V{k) 




(6.21) 



Note that in Eq. (6.19) the first 4 diagrams of lines 1 and 2, and the 4 diagrams of line 3 can 
be absorbed into a mass shift m = l m = l— ♦in the leading contribution represented 
in (6.17). The same mass shift absorbs the diagrams 2-5 in Eq.(6.21) for V. 



C. The zero-mode projector Pq 

We now compute the projector onto the zero-modes 

iki\/(ki) ■ ik2t/(k2) 



okik2 



(6.22) 



1. Po at order 1 

We have already seen that at leading order in the 1/d expansion J^^ k^y(k)^ = d and since 
G{x) = G{0) = ^ = 1 so that with our diagrammatic notations 



(0) 



okika 



-kik2)e-''^«/2e-''^o/2G'(x) 



(6.23) 



Thus the projector Pq subtracts the one-line diagram in O (see Eqs.(5.86)-(5.87)). 



2. Po at order 1/d 

We can now compute explicitely the first correction in l/rf to Pq, using (6.21) for V. It 
is easy to see that the numerator in (6.22) gives all the diagrams in the lines (rows) 1 and 
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2 of (6.19) with n = 1 line between the two points o and x. 

(-ki ■ k2)t/(ki)\/(k2) = O O + 





+ 



+ 




+ 



(6.24) 




+0 



■ vj tj 



+ 




+ 



The denominator is computed from the exphcit form of V given by (6.20)-(6.21). In fact it 
is easy to see, using (6.11), that 



i [ k2y(k)' = l-2^ + 0(l/rf2) = i-^ 



(6.25) 



and since G{Q) = 1 we can write o- 



-£) X + = o 



, etc. We obtain that only (!) 



16 diagrams contribute to Pq; the final result is 




+ 



+ 



+ 




+ 




-0 + 



+ 




-0 + 





cr +o 




+ 



+ 0(l/rf') 
(6.26) 



We note that the denominator gives all the one-line reducible diagrams (with only a single 

line joining o and x) with a tadpole-like graph attached to the line. The diagrams O — ^ ^ — O, 

O and G " although one-line-reducible, are not contained in (6.26). 
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3. An all order argument relating Pq with tadpole graphs 



A simple general argument shows that the denominator in (6.22), given at first order by 
the tadpole diagrams depicted in (6.25), is given at all orders by the tadpole diagrams with 
two truncated external legs attached to the same vertex. In diagrammatic language we shall 
show that 



(6.27) 



1 - o 

The l.h.s. of (6.27) is nothing but 
1-^- 




5 / '^'VW 





i(V"(r(o))>,. with V"(r)^j:^^ (6.28) 



(the two r derivatives pick two legs out of the vertex ^(r)). We can rewrite it as 

V"{r{o))= Jv'\r)6{{r-r{o)) = ^(-k^) i/(k) e^^^(°) (6.29) 

and using the exact equation (6.12) for the instanton potential, we get 

{V"{r{o)))y = [{-k')V{k){^'^^°% = [ k'V{kf (6.30) 



Q.E.D. 

(6.27) implies that P,j will contain all the tadpole chains with tadpole graphs attached at 
the o and x end-points, of the form 



Pn 




(6.31) 



D. Final result for Q = O - Po 

As a consequence the subtracted operator Q = O— Pq is given at order 0{l/d) by the same 
diagrams as those depicted in Eq.(6.19) for O, with the simple restriction that the diagrams 
in the first two lines of (6.19) must have at least 2 lines joining the two end-points (n > 2), 
and that the diagrams in the third line must have at least one non-dressed line joining the 
two end-points {n > 1), while for the diagrams in the fourth line, there is no additional 
restriction (no constraints on the number n of simple lines, ?7, > 0). Using equation (6.20) 
for /ii we can rewrite it as a sum over only 12 graphs (instead of 21!) 




(6.32) 
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E. The determinant D 

We now compute the log of the determinant of instanton fluctuations 

oo ^ 

£ = log(2)) = tr log(]l - Q) = - 5^ T tr(Q') 



k=l 



(6.33) 



1. Diagrammatic representation of the trace 

With our rescalings (see Appendix G) each trace still reads 



ki kfc 



(6.34) 



and with the representation for the kernel Q, we have to compute integrals over k of the 
form, 

/"e-^'^»(-ik'^^)---(-ik'^-)(ik^i)---(ik^") = (-l)"^(2co)^-"^ 5^ 5- ■ ■ ■ <5- . (6.35) 

•'k 



pairing 



Using Wick's theorem we can represent each term by pairing of lines between the left Q and 
the right Q, as already discussed when we introduced the diagrammatic necklace represen- 
tation for £. This is depicted below 



— ^ li'i'^s entering | | n li 



lines exiting 



For instance for m = n = 2 



and for m = 3, n = 1 



,m — 71, ,1 m-\-n \ '■ ^ 

(-l)^(2co)i-^^. 



(6.36) 
(6.37) 

(6.38) 
(6.39) 



The vertical dotted line indicates that no A^-momenta p flows through the vertex, since 
each Q is attached to a different replica of the manifold Ai. With these graphical notations, 
if we represent the kernel Q by the "bead" 



(6.40) 



tr[Q^'] is represented by the k-bead necklace (with periodic boundary condition between the 
left and right dashed vertical lines) 




tr[ 



tr 




, tr[' 




(6.41) 
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2. tr[Q] 

We first consider the term k = 1. We have already seen that at leading order in l/o? 

tr[QW] = ^q[°1, = ^jre-^'^°(e^'«W-l-k2G(x)) = 0{d) (6.42) 
(Q(o) = 0(°) - Po^°^). This can be depicted graphically as 

oo 

tr [QW] = ^ (2co)^- 

n=2 

oo 

n=2 

and one checks easily that the first graph is of order 0{d), the second and the third of order 
(9(1), since each closed loop carries a factor of d, and there are periodic boundary conditions 
between the left and right vertical dashed lines. 

It is easy to see that the trace of the first order correction is of order 0{1) 

tr[Q(i)] = tr[0(^) -Po^^))] = C(l) (6.44) 
and that at this order it is given by the following 12 diagrams 




(6.45) 

This corresponds to a specific, but complicated analytical expression, that we do not write 
here. 

Finally it is quite easy to check that higher order diagrams that contribute to the term 
of order (9((i~^) of Q, will contribute to the terms of order 0{d^~^) of tr[Q]. 





Oil/d) 



3. tr 



We have seen in Sect. V.C that tr 



was of order 0{1). Since 



we 



could have expected that the next order correction Q*-^-* would contribute by a term of order 
0{l/d) to tr[Q^]. We shall see that this is not exact, but that there are nevertheless a lot 
of simplifications, and that a simple subclass of diagrams contributes at order 0{1). 
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In fact there are simply two beads which contribute at leading order to tr[Q^]. These are 



(6.46) 




More precisely, the only 2-bead necklaces which are of order 0{1) are 



(6.47) 

A careful but not difficult analysis shows that each of these diagrams is of order 0{1), and 
that all the other possible diagrams are of order 0{l/d). 
The first diagram contributes by 



■d'{2cor'BiOf 



4c, 



-5(0) 



(6.48) 



(we have taken into account the different contractions of the vertices of (6.38) which give 
this diagram). The four last one give the square opf a single bead amphtude 



(6.49) 



with the single bead amplitude 



dB{0) 

4:Ca 



1 + 



iBjO) 
4co 



1 + /X2 



dB{0) 
Acq 

(6.50) 

where we used (6.5) and (6.6). 

We thus see that the chain of bubbles contributes already to tr[Q^] at the leading order 
0{1). In the large-d limit, e being fixed, since 112 —1 and dB{^)/AcQ — 1 — e/4 we get 



tr 



2(l-^)'+(J-l) +0{l/d). 



(6.51) 



4. tr 



k>2 



The same analysis can be done for the general term tr[Q'^]. Here also the only diagrams 
that contribute to order 0{1) are those of (6.46), and more precisely those with the beads 
of (6.50). It follows that at leading order 



tr 



+ 0{l/d) . 



(6.52) 



We shall comment later on the meaning of the pole in 1/e. 
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5. Summation of the log series 



We see that, except for more complicated graphs coming from the k = 1 and k = 2 terms, 
the whole series (6.33) for £ = log(2D) contains the series — J2k>i ~ ^)'' which can be 
resummed formally as a logarithm, so that (the second term compensates for the missing 
term in the sum giving the log) 



-tr 



(1) 



'4 






2 




- - 1 


+ 




+ log 


[2-11 


e 




L 4J 




e 



Oil/d) . 



(6.53) 



This last series is not convergent if e < 2 and the argument of the logarithm is negative, 
hence 2 has an imaginary part ivr. 

In fact this is not surprising, and is a feature of the model, since we have in fact recovered 
the unstable eigenvalue Amin = 1 — A_ of S" of the Hessian 5", which indeed gives an 
imaginary part ±7r to 2. We show this fact in the next section. 



F. The unstable mode 

It was shown in [14] and in Sect. III.C.2 by general arguments that as long as < e < D 
the Hessian S" [V^'^^^] has one single negative eigenvalue Amin < 0, corresponding to the mode 
of unstable fluctuations around the instanton configuration. 

In appendix F we derive a variational estimate for an upper bound for this Amin- This 
estimates is given by Eq. (F16) and becomes in the large-c? limit 

This is precisely the argument of the log in (6.53). 

Here we show that this is not a coincidence, and that the variational bound A™^ is 
saturated in the limit d ^ oo, e finite, so that the infinite series of necklace diagrams, with 
beads made themselves out of chains of bubbles of (6.50) reconstructs precisely the logarithm 
of the unstable eigenvalue log(Amin)- 

4 

Amin = 2 when d ^ oo, e fixed . (6.55) 

e 

To obtain this result, we shall simply take the following ansatz \E'_ for the unstable eigenmode 

and show that at leading order 

(Il-Q)*_ = (2-4/e)*_ + 0{l/d) . (6.57) 
Let us first compute Q(°)\l/_ 

$-(ki)= / Q£U$4k2) = ^e-^?-/^ / e-^i-k^ / [e'^^^^^M - 1 + kik,G(x)] 

Jk2 ' ^ Jk2 Jx 

= -2co e "i"^"/^ y ( e * '^o ' ^ ^ ^ ' > 



4 eg 2coJ 2co 



(6.58) 
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In the limit d oo since Cq ~ the dominant term is 
Note that we may represent graphically and Q\l'_ by 




(6.60) 



(where the little handle represents Sat) and that the dominant contribution (6.59) at large 
d corresponds simply to the diagram 

q(o)^_ ^ Q^^^^lo^ 5(o)-L^^_. (6.61) 

Note that the rightmost little loop is just o = d/2. 

We can now compute in the large-rf limit the contribution of Q'^^''\l/_. It is given a priori 
by all the diagrams of (6.19) inserted into (6.60). However a careful but easy analysis shows 
that the only diagram which contributes finally at leading order 0{1) is the chain of bubbles, 
that appears already in (6.50) 

4 / e\2 

= - 1 \E'_ when d ^ oo 

e V 4/ 

Now if we consider the graphs that appear in the higher-order terms Q'^''^ of the 1/d 
expansion of Q, one can see that when applied to they also give only terms of order at 
most 0{l/d). Hence we have 

= Q(o)^„ + Q«^_ + C(l/rf) (6.63) 
and combining (6.59) and (6.62) we obtain (6.57). Q.E.D. 



(6.59) 



(6.62) 



G. The zero-mode measure 

Finally, we have to compute the 1/d correction to the weight 22J for the collective- 
coordinate measure for the instanton. According to (3.55), this weight is given by 

W = g""" 

and using the explicit form for V , and in particular (6.25) we get (for D ^ 2) 

g^e-^''^'\l + 0{l/d)) , (6.65) 

where 5i(e) is the coefficient for the mass correction at order 1/d defined by (6.3). 5i(e) is 
of order 0{1) in the large-d limit, its exact value is given by the self-consistent equation 





75 



(6.13)-(6.20). The large-c? limit for Si{e) was already obtained in [14]. It is given by the 
integral 



4-e 



dpp 



log 



1 - 



1 - 



1 - 



4J 



Ap) , Jip) 



which is convergent as long as e > 0. 



2 arcsinh(p/2) 

PV^l +p74 
(6.66) 



VII. THE LIMIT e = AND THE RENORMALIZED THEORY 

We are interested in the renormalized theory in which the UV divergences have been 
subtracted and the limit e ^ has been taken. We have already discussed in Section IV the 
UV divergences and how they are renormalized. Here we discuss this limit in more detail 
and its interplay with the large-c? limit. Our main result is that the 1/d expansion is plagued 
by IR divergences when e = 0, so that the limits d ^ oo and e ^ do not simply commute. 
As we shall see in our discussion, this does not mean that our instanton calculus does not 
make sense at e = 0, but rather that when e = the large-c? limit is of a different nature 
and contains non-analytic terms in d such as logarithms of d. 



A. Minimal subtraction schemes 

To study the renormalized theory at a given dimension d we must first specify a renor- 
malization scheme. We shall use the minimal subtraction scheme (MS) such that the field 
and coupling-constant counterterms in the original action subtract the poles at e = (see 
Eqs. (4.89)-(4.90)). In fact the definition of a MS scheme requires some care. Indeed, since 
e = 2D — d(2 — D)/2 depends both on D and d (the manifold and bulk space dimensions) the 
limit e — > to construct the renormalized theory of a D-manifold ind = dc{D) = AD /{2 — D) 
dimension can be taken in different ways. These different limits corresponds to different 
renormalized theories which differ by a finite renormalization of the field and the coupling 
constant, i.e. these limits correspond to different renormalization schemes. 



1. Definition of the MS-D and MS-d schemes 

1. MS-D scheme: A first scheme is to work at fixed manifold dimension D = D^ and to 
take the limit d dc{D^) = AD J {2 - D^). Then e = (4 - d){2 - D^)/2. This allows 
direct comparison with the field theoretical calculations for SAW and polymers, since 
for Dc = 1 and e = e/2 where e = A — dis the parameter of the standard Wilson-Fisher 
expansion. 

2. MS-d scheme: Another scheme, more natural for 2-dimensional manifolds {D = 2), 
is to fix d = dc and to take the limit e ^ by varying D. In this case e = (2 + 
dj2){D - D,) with D, = D,{d,) = 24/(4 + 4). 
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2. Relation between the schemes 



In both schemes we take as counterterms 



br- 



(7.1) 



e • ■ 2 e 

and the relation between the bare fields r and coupling constant b and renormalized ones r,- 
and bj. is 

r = Z^^h, , b = brfi'ZbZ'^/^ . (7.2) 
We see that both e and d appear explicitly in the second relation for h. At one loop it gives 

1 C2 - rfCi 



b = ^"br 



br: 



In the MS-d scheme the last term gives 

C2 - rfCi = C2{D^, dc) - 4Ci(£'c, 4) 
while in the MS-D scheme it gives 

2Ci(D„4 



C2 - dCt = C2{Dc, dc) - dcCi{Dc, 4 



(7.3) 



(7.4) 



(7.5) 



2-Dc 

We see that renormalization in the MS-D and the MS-d schemes with the same subtraction 



mass scale jj, amounts to a finite coupling-constant renormalization 



^MS-d 



Cl 



2- Dr 



(7.6) 



or equivalently that the MS-d subtraction scale fius-d and the MS-D subtraction scale /ims-d 
are related by 



log 



/^MS-d 



Cl 



(7.7) 



.A'ms-dJ 2 — D C2 — dCi 
Let us also note that we recover the combination of counterterms C2 — dCi that appears in 
the result (4.112) for the coefficient B (defined by Eq. (4.110)) of the UV pole in - for the 
effective action S[V] (see Eq. (4.104)) and for £ = tr log 5" [V] . 



B. Variational mass subtraction scale 



Now for simplicity and in order to study more easily the large-ci limit of the renormalized 
theory we shall work with the normalizations of Appendix G where x and r are rescaled as 
X JTivarX, I" ^var*^^~^^^^ and the coupling constant b is redefined by 6 ^ 'm^a.r^~^h so that 
the variational mass is now set to unity (mvar = 1) in all the calculations. Since the rescaling 



of the coupling constant amounts to g 



-D 



g, this last rescaling amounts to choosing as 



subtraction scale a multiple of the variational mass (/i — * /imvar) in the renormalized theory. 

In this normalization the field and coupling-constant counterterms (as defined in (4.91)) 
Cl = Cl (-Dc5 dc) and C2 = C2(-Dc5 4) in the action become (see Appendix G and in particular 
Eqs. (G30)-(G31)) 



Cl 



-S 



D 



2D 



Co 

do 



1+1 



2Sl 



T[D/i2-D)f 



{2- D)^T[2D/{2- D)] 



Co 

do 



(7.8) 



77 



The logarithm of the renormahzed instanton determinant 2,^ is still given by (4.108) 



Vir) 



(7.9) 



We have seen that in the large-d limit (e fixed), the first counterterm is of order one Ci = 0{1) 
while the second one is exponentially small, C2 ~ 0{exp{—d)). For our discussion of the 
variational approximation and of the large-d limit we only have to consider the wave-function 
counterterm Ci = Ci(-D, dc{D)) = Ci{Dc{d), d) which is given explicitly when e = by 



Ci 



4 (47r 



,D/2 



D V[D/2] 



m-D)i2] 

-r[(D-2)/2], 



2+D 
2-D 



{ATI 



2d 
I 4+d 



dV[d/{A + d)\ 



-r[-4/(4 + d)] 

r[4/(4 + d)] 



(7.10) 



C. Renormalized theory in the variational approximation for finite d 

We first consider the renormalized instanton determinant in the variational approxima- 
tion, but for finite embedding space dimension c?, following the lines of Sect. V.D. We thus 
approximate £ = logdet'[5"] by £(°) = -tr [Q(o)] (as defined by Eq. (6.42)). This gives, 
after integration over k and using (Gil), 



-tr 



(0) 



-k2 



CO 



,k2G(x) 



1 - \eG{x) 



d—2co 



G(x) 

Co 



(7.11) 

To renormalize consistently £ we must take for the condensate ((Vr)^)y in the counterterm 
in (7.9) its value in the variational approximation 



((Vr)^)y - ((Vr)^)_, = d I ^ 



-dcn 



(7.12) 



(we use dimensional regularization) , and neglect the coupling-constant counterterm C2, since 
there is no coupling-constant renormalization in the variational approximation. Thus we 
obtain for the renormalized log in the MS-D scheme 



£(0) 



with 



and 



lim 

e-*0,D fixed 



B„ 



£(°) -dco (gff^L 



D 



log(^r) + log(/xL) 



32 



(2-D)2 



r[(2- J)/2] 
-r[(D-2)/2] 



4 

2-15 



^MS-D 



lim 

e-+0, D fixed 



^ '^MS-D 

(7.13) 
(7.14) 

(7.15) 



Integrating over the angular degrees of freedom of x we can rewrite 

£(0) 

as 

G{x 



£(0) = d- 2cq{d)Sd 



J = f.p. 



dxx 



D-l 



Co 



- 1 



(7.16) 
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where the finite part prescription "f.p." deals with the short-distance divergence at x = 
still present when e > 0. The UV divergence of 3 comes from the short- distance behavior of 
the propagator G(x), obtained from (4.33) 



G{x) =Co-doX^ ^ + 



Co 2 



do 



2{A-D) 

This implies that the integrand in 3 behaves at small x as 

d 



X 



0{X^) . 



(7.17) 



J ~ / dx 


'co' 


2 




Jo 


_do_ 







_d_co 
AD do^ 



A{A-D) 



X 



e+l-D 



0{x 



(7.18) 



The first term gives the UV pole at e = -D, which is subtracted by dimensional regularization, 
and is dealt with by the f.p. prescription. The second term gives the UV pole at e = 0. 
The third one gives a non-singular pole at e = -D — 2, but will be important in the large-d 
limit. Now we use the explicit result (7.8) for Ci, which implies that we can rewrite the 
counterterm in (7.15) as 



dc( 



-2cnS 



0<^D 



Co 

do 



1+ 



dc(D) 



[ dx x" ^ 



(7.19) 



Thus this counterterm cancels the pole at e = 0, but we must notice that since we use the 
MS-D scheme, there is a slight difference between the coefficient of the x*^"^ term in (7.18) 
and (7.19): the first one contains [cq/ dQ\^^'''/'^ and the second one [cq/ d^Y^'^"^'^^^'^ . Since 
d = dc{D) — 2e/(2 — D) this gives a difference of order C(e) for the residue of the poles at 
e = 0, hence a term of order 0{1) in the limit e — ^ 0. We carefully rewrite the expression 



(7.15) for £j^g_j3 
o(o) 



as 



^MS-D 



lim {d+'^ 

e^O,D fixed I AD 





2+D 






— e ~ 


"co" 


2-D \ 


1 - 


'co" 




_do_ 


e 









2coSd 3' 



3' 



POO 

/ dx 








d 
2 

- 1 




Jo 






Co 







Co 

do 



X 



~D-1 



_d_ 

w 



Co 

do 



1+1 



X 



^(1 -x) 



(7.20) 



6'(1 — x) = lifx<l,Oifx>lis the Heaviside step function. This integral representation 
is a priori valid for e > but is now convergent if we take the limit e — 0. We can 
interchange this limit and the small x integration and obtain, using dc{D) = AD/ [2 — D) 
and Sd = 1/(2 - D)do, 



- d 



dx 



Co 



i2-D) 





2+D 




'co 


2~D 

log 


"co" 


_do_ 


_do_ 



Co 



Ddn 



3' 



X 



D-l 



1 - 



G(x) 

Co 



2D 

' 2-D 



Co 

do 



2D 

2-D 



X 



2-D 



Co 

do 



2+D 

2-D 



x^^e{i -x) 

(7.21) 



This last integral over x is UV and IR convergent as long as -D < 2. It can be computed 
numerically. 



For D = 1 we have cq = 1/2, do = 1/2 and G{x) = e '^1/2, and we obtain 

-2 



£ = 4 - 2 J , 



dx 



[1 



1 - x^"^ - x^^en - x) 



C(o) 
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1 

~2 
(7.22) 



Details of the calculation: We compute the integral 



£ 4 - 23' , 3' 



dx 



[1-e-"]] ^ - l- x-^ ~x-^9{l-x) 



First we put a regulator e, and notice that the last term is here to subtract the pole in e 

1 /"^ r _ -2 

3' = lim5(e) , 5(e) = f. p. / dxx^ [l - e""'^ 

£^0 e Jq L'- 



Now 



a(e) = ^(n+l) / dccx^e-"" = ^(n + l)r(e+l)n^^-i=r(e + l)(C(£)+C(£ + l)) 

Ti=l n=l 



Now we use 



and obtain 



C(0) = -^, C(l + e) = i+7E + 0(£) , r(l + e) = l-7Ee + 0(e') 
2 e 



a(e) = --^ + 0(£ 
e 2 



hence / 



1 



But note that Mathcmatica® 5 gives the result directly. For D ^ 1 this integral representation 



d 





2 


4 


6 


8 


10 


12 


14 


16 


20 


^(0) 
■^MS-D 


2 


4.01 


5 


6.60 


9.16 


13.0 


18.5 


25.9 


35.6 


62.7 



TABLE I 



allows easy numerical integration. This gives the following results, presented on Table I and 
Figure 8. Finally, we see on the numerical results that -Cms-d diverges when ci — > oo (i.e. 
when D —>■ 2). As we shall see later, it behaves as 



o(o) 



(2eT^)-^ d^ 



(7.23) 



and this asymptotic behavior is reached as soon as d ^ 20, as shown on Fig. 9. 

Details of the calculation: To compute the integral 3' numerically, it is more convenient to separate the 

integral over x e]0, 1] and over x €z [1, oo[ 



/•l poo 

3' = 3',+3[, 3; = J dx[---] , y,= j dx[---] . 
For the second one we can integrate directly the first term, and explicitly the counterterms and get 



D 



y 7 



POO 




r G{x)^ 


2Z? 

- 1 






Co 





(7.24) 



(7.25) 



For the first one it is better to over-subtract it. in order to improve the integration at a; = and the study 
of the large-d limit, and write 



X . -1 - 



D 



D (2-I?)2(4-D) 



Co 

do 



+ 7: 



(7.26) 
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FIG. 9 £(o)(d) (black curve) in Log-Log coordinates compared to its large-d asymptotics (grey 
curve, straight line). 



y 



/ dx 






2D 
2-D 


'co 


2D 
2-D 


'o 




Co 




da 





2- D 



2 + D 
2-D 



D 



{2-D){A-D) [do 



Co 



(7.27) 



Both integrals 3'[ and 3" are convergent for any D < 2 and have a smooth limit when D ^ 2. 
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D. Do the limit d ^ oo and the limit e ^ commute? 

1. An apparent paradox 

For e > the variational approximation L^^-* for L has a regular large-ci limit. We have 
studied it already in Sect. V.D.2. It is of order 0{d) and is given by the convergent integral 

= d(^l- f.p. dx X [e^^-^)^''^") ~ ^0 ' ^'^''^^^ 

which is the large-c? limit of the integral (7.16) {Kq{x) is the Modified Bessel function of the 
2nd kind). 

This integral can be computed numerically. To study its UV structure, we use the small 
X expansion for the 2D propagator 

Ko{x) = -log(x/2)-7E + ^x2(-log(x/2) + l-7s) + 0(x^) (7.29) 

The integrand in (7.29) behaves at small 

[e^^/2]'-^ (^x'~^ + ^^x'-^ (- log(x/2) - 7£; + 1) + C(x^+^ logx)^ - 1 . (7.30) 

The first term x"^"^ gives the UV pole at e = 2, and is subtracted by the f.p. prescription. 
The second term gives the poles at e = but the logx gives in fact a double pole, so 
that 

There seems to be a discrepancy between this calculation and the results of the previous 
section: 

• Here we take the limit d — * oo, then e — > 0; Sl^°^ has a UV pole oc d/e^ 

• Previously we took the limit e — > 0, then d — ^ oo; £,^°^ has a UV pole oc d'^/e 
Clearly the limits e ^ and d ^ oo do not simply commute. 



2. Resolution of the paradox 

This apparent paradox can be understood if we use the results of the previous section 
to study carefully how the bare quantity H*-*^-* = d — tr[0^^^] behaves when both e — > and 
d — ^ oo. is given by (7.16) and in that limit the dominant contribution is the integral 

3, and more precisely the terms of order x~^ when x — in the integral (7.18) for 3. There 
are two such terms, the dominant one of order x'^~^ (which will give the UV pole at e = 0), 
and the subdominant one of order These two terms combine so that in the large-d, 
small-e limit, we have 

+ with e = 2D--(2-D). (7.32) 
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If we take the limit e ^ with d finite and large (hence 2 — D small but non-zero) the 
first term is singular, and we recover the standard single UV pole, while the second one 
stays finite. Renormalization within the MS-scheme amounts to subtracting the first term 
and we recover 

o(0) _ _ f_ (^oo) 

- 2-D ~ 4D- 

All the other terms contributing to £rcn are at most of order d'^. Thus we recover the fact 
that the renormalized £ is of order 0{d^). 

If we now take the limit c? — > oo with e non-zero but small, we rewrite (7.32) as 

^ - e{e + 2-D) ~ e{e + 2-D) " ^^'^^^ 

and we recover the fact that the bare £, is of order 0{d) but with a double pole when e — 0. 
Thus (7.32) contains both (7.33) and (7.34). 



3. Discussion 

Of course in the full theory, it is the first limit (e — »■ then d — *■ oo) which must be 
considered to study the large-ci behavior of the renormalized theory. At the level of the 
variational approximation framework, from the previous calculations one can show that at 
large-d the variational renormalized log det Hms-d has a regular large-c? asymptotic expansion 
in powers of 1/d, 

-cLId = ty^ + iid^ + ly + ["d"" + ■■■ (7.35) 

with the [°'s real and finite. 

Indeed, setting e = 0, starting from (7.21), using equations (7.24)-(7.27) and the fact 
that Co/do is analytic in 1/d {co/do = 1 + 0{l/d)), one sees that the only possible non- 
analytic terms are the integrals 3'[ and 3'[. Now the mass= 1 propagator G{x) is analytic 
in 1/d ^ 2 — D, except at x = 0, where it has a log(x) singularity when D = 2. Ki x = oo 
it behaves as exp(— x) for any D and it is then easy to see that 3"^ is analytic in 1/d. The 
integral 3"^ given by (7.27) behaves at D = 2 as J^dzx logx and its n*^ derivative with 
respects to D behaves as da; x log" x and is convergent for any n. Hence we deduce that 
yi too has an (asymptotic) expansion in 2 — D ~ Q.E.D. 

Thus the variational renormalized theory at e = scales with d as d^ (and not as rf), 
but still is amenable by a l/c? expansion. As we shall see in the next section, the situation 
becomes more complicated when we deal with the corrections to the variational approxima- 
tion. Indeed, the perturbative expansion studied in Sect. VI is plagued with infra-red (IR) 
divergences at e = 0, in addition to the UV divergences, and we shall argue that this means 
that the renormalized theory contains non-analytic terms such as log((i)'s in the large-d 
limit. 



E. Renormalized theory: first 1/d correction and IR divergences 

1. The IR divergences at e = 

In section VI we have isolated the classes of diagrams in the expansion of the kernel Q 
which give a contribution of order 0{1) in £ = tr log[l — Q]. This analysis is valid provided 
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that e > 0. Indeed, as long as e > the individual diagrams are IR and UV convergent, and 
the summation over the diagrams in each class in also convergent. 

If we now take the limit e — 0, IR problems may occur when summing these diagrams. 

First we consider the diagrams that contribute to tr[Q*^^)], depicted in (6.45). In each of 
the 11 classes of diagrams in (6.45) the sum over the n lines joining the left to the right 
contribute to a similar sum as the sum considered in Sect. V.D.I at leading order, i.e. to 
integrals of the form 



-k^(co-G(x))_g-k2(co)j ^ 2co / ([l-G(x)/co]-''/'-l) X 



These integrals are UV and IR finite when e > (with a finite-part prescription to deal 
with the singularity at x = 0). When e = they are still IR finite (the x-integration is 
convergent at |x| oo since the propagator G(x) is massive, hence exponentially decaying 
at large distance). On the other hand, these integrals are UV divergent when e = (there 
is a singularity at x = which gives a 1/e UV pole) but this divergence is dealt by the 
renormalization procedure. 

Now the second infinite sum in the 11 classes of (6.45) is given by the "chain-of-bubbles" 
propagator of the 1/d expansion 

«(p)- 

given by (6.5) and depicted in Fig.6. Combining the results of sect. VI. A, in particular (6.5), 
(6.6) and (6.10), we easily obtain that in the limit d — > cxd, e finite, this propagator is 

r / e\ / , / N 2 arcsinh(»/2) 

H{p) = 1 - 1 - - J{p) with J p = 7===- = ^i? P 7.36 

L V 4/ J P a/I +pV4 

(we use the notations of [14] for J(p), -B(p) is the bubble amplitude (6.6) at D = 2). For 
large p the UV behavior of H is 

Hip) l + U-ey^ + ■■■ as p oo (7.37) 

and does not raise additional UV problems. For small p its IR behavior is 

H{p) ^ as p . (7.38) 

As long as e > 0, the IR behavior of H is that of a massive scalar field with effective mass 

meff = v^S^ , (7.39) 

to be compared with the variational mass rrivar = 1- However, when e = 0, this propagator 
becomes massless tjIcs = and since we are dealing with an effective theory in two dimensions 
{D = 2), IR divergences occur! Indeed, in the diagrams of (6.45) there are two sources of 
IR divergences: 

1. Firstly, the mass shift depicted by ♦ is given by the solution of (6.20) which involves 
tadpole diagrams with the propagator H{p) at zero momentum p = 0, which gives 
potential powerlike IR divergences ♦ oc i/(0) since 

4 

H{0) = - . 
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2. Secondly, both the tadpole diagrams in the r.h.s. of (6.20) and the other diagrams in 
(6.45) contain internal loops with the propagator H{p). Integration over the internal 
loop momentum gives logarithmic IR divergences since 

/^(P) = |^log(l/e) + --- 



If we now consider the diagrams which contribute to tr[Q'^'^^], k > 2, depicted in (6.47), 
they also contain the zero-momentum propagator H{0) = 4/e. Their amplitudes at large d 
are given by (6.52) and have a powerlike IR divergence in 

This IR problem was in fact first discovered by the authors in [14] , and its significance for 
the calculation of the instanton action studied. It was shown in [14] that these IR divergences 
exist for the instanton profile V{r), but cancel in the first 1/d correction to the instanton 
action iSjnst- As we shall see now, some partial cancellations of IR divergences also occur in 
the contributions of the fluctuations around the instanton, but the first 1/d correction ilren 
to the renormalized fluctuation contribution iiron is still IR divergent at e = 0. 



2. Cancellation of IR divergences in the mass shift 5i 



We first look at the mass shift 6i depicted by ♦ and solution of (6.20). We have already 
computed 6i in [14] and 6i is in fact IR finite when e ^ 0. We refer to sect. 6.5 and 
Appendix B of [14] for the details of the calculation, the final result being given by Eqs. (135), 
(137) & (138) of [14], i.e. the integral 



5i = ^27r^(-log[/7(p)]-(l-|) J(p) 



4-e 
2 

4-e 

2^ JO 

4 ' '■"^ 



dpp ( — log 



dv sinhf 







log 



1-i 



4/ sinhf 



Jip) 



(7.40) 



sinhf 



7.5583 ... for e = . 



This integral is IR and UV convergent even for e = 0, since it behaves at small p as J^l and 
at large p as p-^log^(p). 

To prove the IR finiteness of 6i it is sufficient to rewrite (6.20) as 




X 



Since no momentum flows through the first (vertical) H line we have 

while 



(7.41) 



1 4 

2co e 



d , ^ d 1 

-B(0) = 

2 ^ ^ 2 47r 
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is finite, and the last two diagrams contain log(e) IR divergences. However, one can easily 
check that these IR divergences cancel. Indeed the coefficient of the log is obtained by using 
that as long as /p /(p) is not itself IR-divergent, 

/ /(p) = / H{p) X /(O) + infrared convergent terms 

J p J p 

3 

= — ln(l/e)/(0) + infrared convergent terms . (7.42) 

This means that any H line is to be replaced by ^ log(l/e) and treated as if no 

momentum flows through it. We thus obtain 



while for the second graph 



= log(l/e) I + Oil) . 




The coefficient of the log(l/e) IR divergences is therefore zero, since 

Thus eq. (7.41) for (5i is of the form 

(5i = ^ + C(l)] ^ 5i = C(l) when e = . (7.43) 



3. IR divergences in Q^^^ 

We now perform the same analysis for the first l/c? correction to the Hessian Q, Q*-"^\ 
calculated in Sect. VI. The Fourier transform of Q*^^-* is given by the graphs of (6.32). For 
reasons that will become clear later, let us separate Q*-^-* into 4 parts 

Q(l) ^ Q(la) + q{16) ^ q(1c) ^ Q(ld) (^7_44) 

Q^^'*) is the sum of the graphs which contain the mass shift ♦ 
Q*^^^^ is the sum of the graphs which are "really irreducible" 




(7.45) 
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Q^^"^-* is the sum of the 4 graphs 

Q(-) = c^^—^ + Q) + (^--^ + (^--^ (7-47) 

and Q*^^'^) is the single remaining graph. 

Q^^') = . (7.48) 

Q*^^"-* is IR finite since the mass shift ♦ is IR finite. 

Q(la) _ Q^i^ _ (7 49) 

Q(i^) and Q*^^'^) have a logarithmic IR divergence in log(l/e). By the same argument as 
above, the coefficient of the IR divergence is obtained by removing the H propagator in the 
graph, so that 




(7.50) 



(7.51) 



Similarly for Q*^^'^) we have 

Q('^) = log(l/e)D(i^) + Oil) (7.52) 

ZCo 271 

with 

= + Q) + P , o + Q , • (7-53) 

Finally Q'^^'^^ has an IR pole in 1/e, since the bubble propagator carries zero momentum, 

4. Partial IR cancellations in tr[Q(^)] 

Now come the IR cancellations in tr[Q^^'']. We first consider the second term. We notice 
that the diagrams in (7.51), which contribute in are obtained by two mass insertions 

in the diagrams which contribute to Q^^\ 



(7.55) 
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and more explicitly, since a mass insertion corresponds to (minus) a derivative with respect 
to (where m is the variational mass in the propagators) 



1 



2 dim? 



^ikir(o) „ik2r(x) 



>: 



d 



1 ^2 



m=l 



2 d{rn? 



\2 ^ki.ka 



(7.56) 



m=l 



The effect of the derivative on the propagators from o to x is easy to understand. The 
tadpoles with one or two mass insertions are generated by the derivative acting on the circle 
(6.18) at o or x. 

The kernel D^^^^ is clearly non-zero, but it is traceless for e = 0. Indeed, 



tr 



1 ^2 



2 d(rn? 



■tr 



m=l 



and tr [Q^°^] scales with the mass m like 

tr [QW] = m^-^ tr [Q^] . 
Therefore, since in the large-c? limit, D = 2, we have 



tr [D^^")] 



e(2 - e) 



tr 



(7.57) 



(7.58) 



(7.59) 



and formally^ tr[D(^'')] = when e = 0. 

Similarly we can compute the IR coefficient D^^^^ for the 3rd term and its trace. We 
obtain easily the explicit result 



tr[D(i^)] = 47rcoe2(l-e/4) = de^A 
It is also zero when e 0. 
Details of the calculation: We have 



(7.60) 



a 







^ — o 






o — ^ 




— 


— 



-A^) (-k^)(-k^)(-kiMe-(^?+^i)-/^ 



2 I 4^' (-k?)(-kik2)e-(^?+^i)-/^ 



1 / 1 



2 U^) (-k^)(-kik2)e 



-(k2 + k2)co/2 



-kik2) e 



-(k2+ki)co/2 



we integrate over ki = — k2 to obtain the trace and we use the fact that 

d 



4:Cn 



47r(l - e/4) 



Q.E.D. 



® One must be cautious since tr [Q'''-'] has an UV pole at e = 0, so there is a mixture of IR and UV 
singularities, that we shall discuss later. 
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The final result is therefore that the IR divergence of tr[Q'^^)] comes only from the last single 
diagram in (7.48), which gives Q*^^'^)! It is a single IR pole in 1/e, since 



hence 



tr[Q(i'^)] = -fl-i) . (7.62) 



e V 4 



5. IR divergence and the unstable mode 

We now look at the IR singularity in the other terms of the expansion for £, 



oo ^ 

£ = trlog[]l-Q] = -5^-tr[Q^] . 



k 

k=l 

We have shown in section VI. E that the next terms tr[Q'''] are of order 0{1) and can be 
computed explicitly at that order, since they are given by the contribution of two diagrams 

only (see (6.46)) in Q, namely the single bubble diagram CC^^X' "which is IR finite, and 
the diagram (^^' "which is nothing but the IR divergent diagram of (7.48). Thus 



we see that all the IR divergences of £,^^\ i.e. the term of order 0{1) in the 1/d expansion 
of £ are contained in the last term of (6.53), namely in the summation of the log series 

£(1) log ^2 — + IR finite (but U"V divergent) term when e ^ . (7.63) 

Now we have shown in Sect. "VI. F that this IR singular log(2 — 4/e) is nothing but the 
contribution of the smallest (and negative) eigenvalue of the Hessian associated to the 

unstable eigenmode (dilation) for the instanton 

Amin = 2 - ^ < . (7.64) 

The conclusion of our analysis of the IR divergence of £ = logdet'[i5"] is that, at least at 
order 1/d, the IR divergence can be attributed entirely to the contribution of the smallest 
eigenvalue. This is in fact quite natural, since IR divergences must come from the large 
distance properties of the fiuctuations around the instanton configuration. 



6. A conjecture for the large-d behavior of the unstable mode 

This IR divergence in our large-d estimate of the negative eigenvalue Amin for the instanton 
Hessian does not mean that Amin is IR singular when e = 0, but rather that Amin does not 
behave in the same way when d ^ oo, depending on whether e > or e = 0. 

• If e > 0, we have seen that Amin = C'(l) when d ^ oo. 
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• If c? is finite and we take the limit e — ^ 0, we have also Amin = 0{1), as can be checked 
explicitly for the case d = 4,e = 0, where we recover the classical (pf instanton for the 
0(n = 0) model for SAW. 

• Therefore we expect that the IR pole in (7.64) means simply that when we first take 
e = 0, then d — > oo, Amin is no more of order but becomes infinite (Amin oo). 

In fact we conjecture that for the renormalized theory, Amin scales as d 

Amm(e = 0,c/) ~ C(rf) when d ^ oo (7.65) 
by analogy with the behavior of the leading term £,^°^ which is found to behave as 

- ■ - when (i — > oo then e — 

'd ' ^^-^^^ 

- ■ d when e — then d ^ oo 
e 

The first ^ being an UV pole, and only the last ^ ~ d being IR. Even if this form of the 
conjecture is not correct, it is clear that once again for the unstable mode the limits e — > 
and d — > oo do not commute. 



VIII. CONCLUSION 

In this paper we have shown how to compute at one loop the fluctuations around the 
instanton in the self-avoiding manifold model, and how this is related to the normalization 
for the large order asymptotics for the SAM model. We have shown that the perturbative 
counterterms which make the SAM model UV finite in perturbation theory do renormalize 
(at one loop) the instanton contribution. We have constructed a systematic 1/d expansion, 
and studied the first terms of this expansion and the interplay between the 1/d expansion 
and renormalization. 

Although we have obtained many results in this article, and checked at one loop the 
consistency of the instanton calculus for the SAM model, several points deserve further 
studies: 

• It would be interesting to get a better understanding of how to resum the IR diver- 
gences present in the 1/d expansion for the renormalized theory at e = 0, or to find 
another approximation scheme which does not suffer from IR divergences. 

• We have checked that the instanton factor obtained by our method is for D = 1 (self- 
avoiding walk) equal to the factor obtained by field theoretical methods. However, it 
would be interesting in this case to compare the approximate result that we obtain via 
the large-d limit with the exact result (as was done for the instanton action in [14]). 

• A practical application of the theoretical results obtained in this paper would be to 
compare our large-order asymptotics with our explicit calculations at 2-loop order 
for the scaling exponents for the SAM [22, 23]. Since the non-perturbative effects 
become small when d is large, it is expected (and checked numerically) that the 2-loop 
estimates for the critical exponents are reliable for large d. Such a study would help 
our understanding of the domain of validity of the 2-loop calculations, and perhaps 
suggest better resummation procedures than those used previously. 
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• For renormalized local field theories, in addition to instantons, other contributions 
occur in the large-order asymptotics, denoted renormalons. They are associated both 
to the short-distance behavior of the theory (UV renormalons) and to its large-distance 
behavior (IR renormalons). We expect that such effects occur also for the SAM model 
at e = 0, since for D = 1 it is equivalent to the (p^ theory, but it is not known how 
to treat these renormalon effects (if they are present) in the framework of the SAM 
model, which is a multi-local theory. 
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APPENDIX A: Measure and normalizations for the functional integral 

In this appendix we precise the normahzation for the functional integration over the fields 
and the treatment of the zero modes. 



1. DeWitt metric and measure for the functional integral 

We consider the free membrane model. The functional measure P[r] is normalized as 
follows. We start from the DeWitt metric G over the manifold configuration space C = {i'(x)} 

^^^^'^^^ = i/^^^'^'^^wr = §^\Mi. (Al) 

II ... II 2 is the L2 norm over Ai. This metric depends explicitly on an (arbitrary) normaliza- 
tion mass scale fio- 

The corresponding measure is defined (formally) by V[r] = nx'^'^''(^) V det G. This cor- 
responds to the normalization 

Jv[r]exp(^-^ J d^xr(x)2^ = 1 . (A2) 

With this normalization, a quadratic form A with kernel A"'^, i.e. (Ar(x))" = 
/d^yA%(x,y)r''(y), yields 

y'p[r]exp(^-i^d^x^d^r(a;)A(x,y)r(y)^ = det [A//.g] . (A3) 

To evaluate the partition function for the free membrane 



Zo = j V[r] exp (^-1 ^d^x Vr(x) 

we must treat separately the zero modes ro(x) = Tq of the scalar Laplacian Ax over M and 
the fluctuations 7 orthogonal to the zero mode, G{rQ,T) = 0. Let G^^^ be the metric for the 
collective coordinate Tq of the zero mode induced on the "moduli space" of minima of the 
action r(x) = Tq by the DeWitt metric 

G(5ro,5ro) = ^ Jji^x\6ro\' =: Gi?5ro"5r^ G'^^ = ^Yo\{M)6., . (A4) 



Hence the measure is 



2 

^Vo\{M) 



(A5) 



d/i(ro) = d\^det{G'^^) = d\ 

The integration over the modes 7 orthogonal to the zero modes gives 

(det'[-A/^^])''^/^ , (A6) 
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where det' is the reduced determinant, that is the product over the non-zero eigenvalues of 
the operator — A//iQ. Hence 



d/iM (def [-^/^^V\) 



with the partition function for the marked manifold Zq 



det' 



-A' 



27r Y 
Yo\{M) J 



-d/2 



(AT) 



(A8) 



2. Zeta-function regularization 

The det' requires UV regularization for its definition. We use the standard zeta-function 
regularization (see for instance [13] and [24]). 



log(det'[-A//ig]) = tr'(log[-A//i^]) = - 



(A9) 



s=0 

2 



where the zeta-function for the operator A = —A/fi^ is defined by the sum over the 
non-zero eigenvalues Aj 

Ca{s) = J2 (AlO) 

for Re(s) large enough, and by analytic continuation down to s = 0. tr' means the trace 
over the subspace orthogonal to Ker(yl) ( w.r.t. to the metric G). 

The operator —A scales with the internal size L of the manifold Ai as Therefore 
scales as 

as) = (Lfioras) , (ah) 

where ({s) is a scale invariant zeta function which depends on the shape of Ai but not on 
its size L. 

If there is no global conformal anomaly, ({s) is analytic around s = and C'(0) = 
C'(0) + 2 log(L/io)C(0). Moreover, for any such A, one has 

Ca(0) = — dim(Ker(A)) = — number of zero modes of A . (A12) 

Indeed, one can show that if A has no zero-mode, for instance A = —A + m^, then Ca(0) = 
tr(l) = (this is analogous to the celebrated rule 5(0) = in dimensional regularization), and 
if A has zero modes, Ca{s) = lim^^Q[(A+e{s) — A^e~*], therefore Ca(0) = lim£^o[CA+e(0) — 
A^] = — A^. The Laplacian A has one zero mode, and therefore 

C(0) = -1 , C'(0) = C'(0)-21og(L/io) . (A13) 

Using the fact that the size of the manifold is defined as 

L = Vol(Al)^/^ (A14) 

we obtain for the partition function 

-1 d/2 



Zo = L'^(2-D)/2 



2n 



(A15) 



The dependence on the mass scale /xq used to define the measure T>[t\ has disappeared, as 
expected in the absence of a conformal anomaly. 
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3. Conformal anomaly 

It is known that there is no conformal anomaly if 

1. D = 1 and the manifold has no boundary (closed loop). This corresponds to a ring 
polymer. 

2. D = 2, the manifold has no boundary and has Euler characteristics x = 0. This 
corresponds to a closed membrane with the topology of a torus (or a Klein bottle). 

3. D non-integer. The model is defined by dimensional regularization, as detailed in [13]. 
This is the relevant case for the e-expansion. 

If there is a conformal anomaly, ({0) ^ —1 and there is an additional power of L/io in the 
partition function which depends explicitly on the scale /iQ. For instance for D = 2 
(membrane) it is known that 

(^(0) = — 1 H — X 1 with c = 1 the central charge for the free boson . (A16) 
6 

Hence 

Zo oc {L^i^Y^I^ . (A17) 

APPENDIX B: Integration paths for the functional integration over V[r] 

In this appendix we discuss in more detail via the steepest descent method the functional 
integration over V{r) and the relative position of the instanton V^'^^* and of the integration 
contour over V, as the argument 6 of the coupling constant b varies in [— vr, vr]. This is 
required to treat properly the contribution of the unstable eigenmode of the Hessian iS"[y] 
for the instanton. 

For general 6 G [— vr, vr] we know from (3.28) that the functional integral for the rescaled 
potential ^(r) is normalized so that 

{g is real positive). The effective action for V is given in (3.35) 

p-ie r 

Se[V] = S[V] - —J d'^rV-(r)^ (B2) 

and for large V is dominated by the last term J V^. The steepest descent integration path 
for \^(r) in C is such that (at least for large \V\) 

Arg{V) = ^ (B3) 

(see figure). Thus it turns anti-clockwise from the positive real axis for 6 = — vr to the 
imaginary axis for 6* = to the negative real axis for 6 = ir. 

For general 6 the instanton V^^"^* is an extremum of iS6)[y]. For negative coupling {6 = 
±vr), the instanton is known. It is the solution found and studied in [14], V™^^ = V™^^; it is 
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FIG. 10 Integration path for V as a function oi = Arg(6) 



real and negative; it lies on the steepest descent integration path given by Eq. (B3). Let us 
start from the case 6 = ir, i.e. b lies above the discontinuity along the negative real axis, and 
look at what happens when 6^0. From the solution for the instanton at 6* = vr, V^^^{r), 
the instanton for general 6 < tt , Vg™^*(r), is obtained by analytic continuation from real r to 
complex r. Indeed, we know from [14] that for a general V{r), under a scale transformation 



y(r) Vxir) = A^V(Ar) 
the two terms in the effective action S[V] scale respectively as 

S[V,] = X^S[V] ■ 



\2-D 



(B4) 



(B5) 



If we assume that the instanton y°'^*(r), obtained in [14] for real r, can be continued 
analytically to complex r's, then it is enough to take instead of a real scaling factor A a 
complex phase factor 

A = e'" 



and to choose as phase 



to know that 



D 



2{D-e) 



with u = 00^ an extremum of ^^[y]. Therefore the instanton for 6 < tt is 



(B6) 
(B7) 
(B8) 

(B9) 



It is clear that this instanton Vg'^^^ is now a complex field configuration, since it involves both 
a "global Wick rotation" in r space and the multiplication by a global phase. In particular 
for ^ = (real positive couphng constant) the instanton is 



(BIO) 



96 



The same argument applies for 9 G ]— 7r,0]. If we start from the same real instanton at 
9 = —7T and deform it to ^ = we obtain another instanton, which is the complex conjugate 
configuration Vq^q of the instanton obtained by starting from 9 = n. 

How is Vg^^ located with respect to the steepest descent integration path over V^(i')? 
Rather than considering the functional integral over V{r) for real r's, it is more convenient 
to rotate the space coordinate r in the complex plane. This is equivalent to deforming the 
time contour in the complex plane when dealing with time correlation functions in finite 
temperature QM and FT. Consider as bulk-space coordinates r defined as 

r = e'^^o r (Bll) 

and make the change of variables in the functional integral for V 

V{r) V{r) = V{r) , f real . (B12) 

The functional measure becomes ^^[1^], the measure for V, and from (Bl) it is normalized 
so that 

J V[V] exp (jj^^'"^""' j d'^f \/(f)^^ = 1 . (B13) 
The steepest descent integration path for V is therefore the line with argument fi^ 

Arg(y) = Qj = = TT + (B14) 

(remember that e = 2D — d{2 — D)/2). In the new variable V the instanton differs from the 
original real instanton V^^^^ by a pure phase (see Eq. (B9)) 



Vr\r) = e^('^-^)^r'''*(f) . (B15) 

Since V^'^^^ is real and negative, its argument is tt and therefore Vg"'^*(f) has a fixed argument 
(independent of f) fijf*^* given by 

Arg (t/r'*(f)) = hf'^ = TT + (71-9) . (BIG) 

For 9 < TT, il^^^^ is larger than ilg 

Q^'^ > Qj a 9 < 71 and e < D . (B17) 

This means that the instanton lies below the integration path for V, see figure 11. When 
9 ^ TT the integration path becomes the real axis (with the standard orientation from — cxo 
to = oo), while the complex instanton V^'^^^ becomes the real (and negative) instanton 

With this result the steepest descent integration prescription for the unstable mode 
around the instanton at 6* = vr is fixed. We boldly denote by V this mode. The inte- 
gration path from to — oo has to start from = (the real vacuum, minimum of the 
action 5[y]), go on the real negative axis up to the instanton V™^^ < which is a local 
extremum of S\V] with action 5™^* = 5[V'™^*] > S\V = 0] = 0, then "turn right" (see 
figure 12) in the upper half complex plane in order for the action to continue to increase, 
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ft — TT 

vf — 71 
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ReV 


/ e<ji 




= 





FIG. 11 The integration path for V and the instanton y™*^* for 
have set e = 0). 



TT, < ^? < vr and 9 = (we 



Instanton 



classical vacuum V=0 



e=-7i 



FIG. 12 Steepest descent integration paths for the unstable mode for 9 = ztyr. 



while leaving the instanton below, then go to — oo. The first part of the contour (from V^'^^^ 
to oo) contributes only to the real part of the partition function Z for negative coupling 
(and is dominated by the classical vacuum V = 0). The second part of the contour (from 
— oo to ^"^'^*) contributes to the imaginary part of Z; in fact the dominant contribution to 
the imaginary part comes from half the Gaussian integral in the imaginary direction at the 
instanton 



dV 



and gives a factor 



*+ioo 



i ^|det 



-1/2 



(B18) 



(B19) 
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when compared to the full contribution of the Gaussian integral (the factor — i comes from 
the integration path, the factor 1/2 from the fact that we integrate the unstable mode from 
ioo to 0, not from ioo to — ioo). 

The same argument shows that for 9 — tt the instanton is above the real axis. Therefore 
for the unstable mode the steepest-descent integration path has to stay below the real axis, 
with factor 

i^ |det(5"[r°^*])|~'/' . 
These are the results used in Section III.C.2. 

APPENDIX C: SAW: D = 1 SAM versus 0(n = 0) field theory 

In this appendix we recall the "Laplace-De Gennes" equivalence between the zero- 
component 0(n = 0) </)^ field theory and the (weakly) self-avoiding walk model, which 
corresponds to the case D = 1 for the SAM model. The first part of this appendix (sect. 
1-3) is basically textbook material, recalled here to fix the notations and the normalizations. 
We then show that the standard instanton calculus for the 0(n = 0) model gives the same 
result as our instanton calculus for the SAM model in the special case D = 1. This provides 
an important check for the consistency of our method. 

1. Free Field and Brownian walk: 

The action for the scalar free field in ci-dimensional space is [note the factor 1/4, which 
is not the most commonly used '^) 

4(v.0)^ + —{<pr . (ci) 

The 2-points correlation function is 
Go(ri,r2;m2) = mn)<l>{r,))o = {n\ . , |r,) = / dL e-^^^r^le^^/^^) , (C2) 

and is the Laplace transform with respect to L of the heat kernel /('(ri, L) = (ri|e'^^/^|r2), 
which admits the random-walk representation 

ir(ri,r2;L) = (rile^^/^jra) = ! D[r] e"/^ (C3) 

J r(0) = ri 
r(i) = r2 



^ Two choices of normalizations are convenient for polymers: Here we use S(i[(f>] = \ {^r(t>f' + ^4>'^^ 
which corresponds to having the polymer action S'poiymcr — Xc 5 (^''(^))^- The other convenient choice 
is to use So[4'] = /r 5 (Vr^)^ + ^0^, which corresponds to S'poiymcr = Xc i (^''(^))^- This is the choice 
most often taken, see e.g. [4]. Here we employ the first choice, since we want to use the most convenient 
normalization for the polymer action. We also note that for both choices, e~ , with L the length of 
the polymer is the weight in the Laplace-De Gennes transform. 
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with r = dr/ds. To check the normahzation, use the semi-classical estimate for the small-L 
limit of the r.h.s. of (C3), K ^ exp (— |rp/2L) and check that 2dK/dL = AK. In particular 
at coinciding points 



dL e 



-Lm? 



LA/2 



ri) 



-A/2 + m2 

and the heat kernel at coinciding points admits the closed random-walk representation: 



(C4) 



V[r] e~JL2 

{0) = r{L) = ri 



\D=l:L 



(C5) 



It is the partition function for a closed 1-dimensional membrane (i.e. a closed polymer, or a 
loop) with length L, attached to the point ri. Similarly for the one-loop connected diagram 



- [(0(ri)V(r2)^)o - (0(rO^)o(0(r2)^)o] 



(ril 



-A/2 + m2 



k2) 



dL e 



-Lrn^ 



dL e-^'^'L-^ n'i\r,,r2 



dL, I P[r]5^(r(0) - ri)5'^(r(Li) - ra) e^^^ ^(O^d- 

(C6) 



D=1,L 



This means that the first derivative w.r.t. of the l.h.s. of (C6) is the Laplace transform of 

(2) 

the 2-point correlation function T^q for a free closed loop with length L. Similarly connected 
correlators of a product of 0^ operators are associated to A^-point correlation functions 
for the closed loop. 



2. SAW and 0(n = 0) field theory: 

It is well-know that this equivalence extends to the Edwards Model, defined with the 
normalizations as in (2.12). The 0(?7,)-invariant 0^ model is defined by the action 

= I d'^\{^r$y + ^(0 ^) + liPr (C7) 

with a n-component real vector field: 

= {0'^;a = l,n} (C8) 

t = is the squared mass, b is the coupling constant. 

The model is equivalent to the Edwards model of polymer with (weak) 2-chain repulsive 
contact interaction, as defined by the model of (2.12) for the D = 1 case. The equivalence 
holds thanks to the very same Laplace transform between correlation functions as in the 
free case (6 = 0). It is valid to all orders in perturbation theory, that is as an asymptotic 
series expansion for small b. The operator |0^(r) is represented by a (5-distribution, or more 
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formally^ 

> j d^x5(r(x) -r) . (C9) 



2 

For instance, for the 1-point correlators we have 

2 /I 



l'iS;i\2*(''' ; = /„ die-"Z(r01, (CIO) 



and for the two-points correlators 
n^o n \2 2 

etc. . . 



/ dLe-^*L-l7^(2)(rl,r2)|^ (Cll) 
Jo 



3. Instanton calculus and large-orders for the 0(n) field theory 

We now recall the principle of instanton calculus and large-order estimates for the scalar 
0{n) (j)^ field theory, following the standard references. This example is useful, since in the 
limit of n = it describes polymers, i.e. D = 1 membranes. The field is a n-component real 
vector field 0(r), = {(j)""; a = l,n). The action is the 0(n)-invariant 0^ action 

S[^] = I d'r^^Vjy + ^(0 ^) + (C12) 

t = is the squared mass, b the coupling constant. We are interested in observables O[0] 
which are local monomials in with degree do in 0, the simplest being the energy operator 
E 

E[ri] = (0)'(ri) ; degree ) = d^, = 2 . (C13) 
The expectation value for the observable O is given by the standard formula 

(O) = y"l?[0]O[0]e-^[^Y Jv[$]e~^^^K (C14) 

We are interested in the large orders of the perturbative series expansion in b. As we 
have seen, we have to use the dispersion relation in the complex-6 plane and consider what 
happens for small b close to the negative real axis, where b is complex and its argument is 
close to ±7r. Therefore we rescale the field 

= l^r^/V ; with e = Arg(6) . (C15) 

This gives 

5[0] = ; SoM = J d'r^-{Vr^f + 2) + ^(^ r (C16) 



^ Note the factor of ^. Intuitively it is there to compensate for the fact that the two fields of (j)^ can be 
contracted in two different ways. This also leads to a relative factor of 4 between the interactions (2.12) 
and (C7). 
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so that 



(O[0]) = \b\-^{OM)e,w = \b\- 



JV[if]0[if]e 



'SeM/\b\ 



(C17) 



For small positive b the functional integral is dominated by the constant classical saddle 
point 

0o{r) = 00 = (C18) 

constant and absolute minima of Se for ^ = 0. The functional integral is in fact well defined 
as long as — vr < 9 < +7i. Now along the cut at 6 < 0, that is for 6 = ±7r, another real 
extremum of the action Sg becomes important, the instanton 



(Pi{r - rojuo 



(C19) 



The instanton is characterized by its position tq in space, and its orientation uq in the 
internal ra-dimensional space [uo being a unit vector in R*^). v^j(r) is the real finite-action 
solution of the equation 

-Ar^i + t^i- = , (C20) 

which is rotationally invariant around the origin (i.e. depends only on |r|) and is non-zero 
except for |r| — >• cx) (or equivalently > 0, this is enough to define it uniquely). 

The contribution of the instanton in the functional integral is at one loop proportional 



to 



e W 



[Def [S'M)]] 



^1/2 



(C21) 



The measure for the collective coordinate ro (the position of the instanton) is easily obtained, 
since the metric is 



hab 

Hence the measure is 



1 



2n\b\ 



1 



27r\b\d 



llVf^illa Sab ■ 



d/i(ro) 



d-^ro 



1 



_27r|6|d' 

The measure for the internal coordinate Uq (the orientation) is 



(C22) 



(C23) 



dur 



2TT\b\ 



n-l 
2 



27r|fe| 



\0i\\l 



n-l 
2 



(C24) 



with duo the standard measure on the unit sphere iS„„i in M". For the 0(n) invariant 
observables which are of interest to us, the integration over Sn-i can be performed explicitly, 
giving the factor f2„ (the volume of the unit sphere in R") 



d/i(no) 



2 7r- 



T{n/2) 



(C25) 



Note that this volume factor f2„ vanishes as n when n 0. However, since 0(n) invariant 
observables such as (p^/n behave in the background of the instanton (pi as /n and are 
therefore of order the factors n and compensate to give a finite n — > limit. 
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Now the Hessian is a n x n matrix in internal space 



Hence for the instanton background 0i 
of the longitudinal operator 

SI' = 

times n — 1 transverse operators 



(C26) 



LPi uq the Hessian can be written as the product 
3 



-A + t--cf^ 



(C27) 



(C28) 



S" = SI' ® (5^')" ■ Note that SI' has d zero modes ipl^ = d^^pi and that S'l_ has one zero 
mode ip^^ = y:>i, so S" has d + n — 1 zero modes. Thus 

det'(5") = det'(S'/') det'(S'l')""' (C29) 

and in the n = limit 

det'(^" 



det'(5'/'^ 



I 71=0 



det'(5'"^ 



(C30) 



SI' has one (and only one) eigenvector ip^ with negative eigenvalue < —2t. Therefore 



det SI' < 9 



To understand which signe must be chosen for the square-root of the negative determinant 
[det(S"')]~^/^ ( +i or — i ?) we have to consider the steepest-descent integration path for the 
global (f variable. It has to go from Aigif = n ± 6/A to =F^/4, hence for 6 = +7r it is as 
depicted on Fig. C.3. 

Hence the instanton contributes to the imaginary part of an observable by a coefficient 



//C//M-1/2 



T^\det' {Sl')\ 



if 



e = ±7T 



(C31) 



The rest goes into the real part together with the contribution of the classical vacuum ipo. 

Putting things together, and using Eq. (3.12) we obtain for the imaginary part at 6 < 
of the e.v. of the 0(n) invariant observable O and for arg(6) = 6 = ±7r 



lm(0) 



det'iSM) {det'{Snv>^)]) 



det(5"[v9o])" 



27r\b\d 



2 

«ll2 



2vr|6| 



I<^i|l2 



71-1 

2 



ftr. 



d\ (O[v9.[ro]]-O[v9o]) 



(C32) 



^ This is true for rf < 4, for d = A there is no instanton solution for t > 0, for t = there is an instanton 
with an additional zero mode ipi ^ = (rV,. + corresponding to the scale invariance of the masslcss 
theory under scale transformation (p(r) Aiy9(Ar). The instanton at d = 4 is obtained from the instanton 
for d < 4 by taking the limit d ^ A, t (x A ~ d. 
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6 



(b) 







(a) 





(c) 



FIG. 13 Steepest descent integration path for the global (p variable for ^ = (a), 6 = tt (b) and 
9 = —TT (c). The black dot represents the classical vacua if = and the white dots the 2 instanton 
saddle-points ip = ±ipi. 



while of course 

Re((0)) = 0((^o) • 
In particular for n = and O a product of energy operators E defined as 



(C33) 



OiN = lim-(0(ri))' 
02[ri,r2] = lim i(0(ri))'(0(r2))' 



and omitting the i subscript for "instanton" , we obtain 



Im(Oi) = |6|-^e-W^' 



det'Sl' 



det'Sl 
det'Sl" 



. l-d 

27r — 



2 

«ll2 



det'S" 



(27r) 



d 



(C34) 
(C35) 



(C36) 



(C37) 



denotes the usual convolution product f * g{r) = f dr' f{r')g{r' + r). Note also a few useful 
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results 



^ 



d''r(V(^,)' 



(C38) 
(C39) 

(C40) 



4. Instanton calculus for the SAW model of polymers 

Since the Edwards model for SAW is the inverse Laplace transform w.r.t. t = m? of the 
0(n) model, instanton calculus must take into account this transformation and is (slightly) 
modified, as is explained here. 

The action for the SAW model is 



s[<p\ 



4^ ^^ 2^ 8 



k'2\2 



(C41) 



We have by inverse Laplace transform, for a closed polymer of length L 

r+ioo 1. 

2{r-L) = j |ie^'(Oi(r;t)) 
So we consider the effective action 



"+100 ^. 

L / — e^*(02(ri,r2;t)) 
^ 2i7r 



tL 



(C42) 
(C43) 

(C44) 



To factorize \h\ and L we must rescale both 0, t and r with 

0(r) = y2v5(r') , r=|6|^L^r' 

The action becomes 

= 2[|6|Lr^5[v^,r] ; 
and the effective action is now 



t=\h\d^L^T (C45) 



(C46) 



S[^,t] = S[ip] 



The effective coupling constant is 



S\f] 



'JcS 



\b\U 



2 

d-2 



(C48) 



instead of 
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The instanton is given by the saddle point equations 



and the Hessian is 



A , 

^ — if + Tif - = 



S" 



S" 
^ 



(C49) 



(C50) 



It can still be separated into its transverse part, which is — 1 times the transverse operator 



times the longitudinal part 

SI' = 



SI' V 



SI' 



A 



+ r - 3(p^ 



(C51) 



(C52) 



which has d translational zero modes, namely the 







smce 



\ J d^ip^ = 0. It is then easy to show, denoting by Pq the projector onto the kernel of S" 
generated by its zero modes, and defining the "inverse" of 5J" as 



S," 



(S/' + Po)"' - Po 



(C53) 



that 



10 



det'S;' = det(^'^'"+/° = det(5;' + Po) det ( - ^'(Sf + Poj^V) 



,/ c<// 





' 1 " 


('■ 


[sn 



But it turns out that 



[sn 



■ ^ 



2r 



[r- V + l)<^ 



(C54) 



(C55) 



Indeed, one can check that rVip as well as ip are orthogonal to ^/^^ and an explicit calculation 
shows that 

— + r-3(^2 (r-V + l)<^ = -2Tip . (C56) 

It follows that the additional factor in the determinant (which comes from the integration 
over r at the saddle point) is 





' 1 " 




[Si"\ 



■ V 



1 

27 



d-2 
At 



d-2 
At 



(C57) 



^'^ This is an application of the general formula for the determinant of bloc square matrices det 
det [A] ■ det [D - CA-^B] 



A B 
C D 



106 



r is > if c? < 4, so this factor is negative for 2 < d < A, but the integration path over r is 
also imaginary (it goes from — ioo to +ioo). Thus the integration over r gives the factor 



(C58) 



[coming from the measure y~) times 



27r 



OeS — 



1/2 



(C59) 



coming from the Gaussian integration. This gives 



d-2 



vr- 



-1/2 



\b\^L^ 



d-2) 



(C60) 



For the observable (p'^/n the integration over the zero modes gives 



n 



d _ -2 
d-2 



(C61) 



times the measure and determinant factors 





-1/2 


niv^in 


d/2 


det'Sf 


-1/2 


2 4-d" 
TT 6 d-2/,d-2 


1-d 
2 


iv<^ir 


d/2 








_27r6efr'^. 




det'5^ 






d 







-1/2 



Putting things together we get 



1 -d -2d 

ImZ(r;L) = T -L— [\h\UY~^ e 



(C62) 



-2 

(25-r)[|b|L^]^ 


det'S*/' 


1 

2 


nivv^iiq 


d 1 
2 /rf-2\"2 




det'S"!' 




nd 


V 4r J 



(C63) 

Remember that |6|L^ is dimensionless. We show below that with our normalizations and the 
equation for the instanton. 



T={A-d)S ; 
This simplifies slightly (C63) 

ImZ(r;L) = ^ 1 ^e-^"^"^)^ [I'-H 



||V<^f = 2dS 



(C64) 



-2 



dei'Sl' 



2S 



TT 



4:{4:-d)S' 

d-2 



(C65) 



5. d ^ 4 limit and scale invariance 

The 0(n) model at (i = 4 (e = 0) becomes scale invariant and the instanton has an 
additional zero mode associated with dilations. For the SAW model nothing special occurs 
when e •— > (as far as global scale transformations are concerned). Here we show in detail 
how the dilation zero mode is absorbed in the transformation 0(?7,) — > SAW, and the form 
of the instanton and of the large-order results at e = 0. 



107 



We first derive a few exact results. If ip is the instanton, solution of 
and the Hessian is SJ" = — y + r — 3^9^ then 



SI' 



SI' 



(C66) 
(C67) 



(rV + 1)^) = I ^'{rV + l)v. = 1 (1 + irV) = ^ / 



-2r {cp' 



(fi) = -2t {if 



[SI' 



r ((^1^) = r / 



The instanton action is 



S 



T 



Hence 



T 



{V^f = 2dS 



4S . 



(C68) 
(C69) 
(C70) 

(C71) 
(C72) 



Now remember that r is fixed by the normalization J if"^ = 1. In the limit d ^ 4, we must 
take the limit r — to get the finite-action instanton. The general solution of the equation 

-Alp = 2ip^ at d = 4 is 



2ro 



with corresponding action Sd=4 = -tt^ . 

3 



(C73) 



ro is the instanton size. The size is arbitrary for the massless d = 4 theory by scale invariance. 
For r ^ rQ, if satisfies the linearized equation —p" — + 2Tip = and is a Bessel function 



<^(r) cx ri-'^/2j^,/2_i(v^r) (C74) 

and is such that 

99(r)~Cr^"'^ for Tq < r < l/^/r ; (p{r) oc e''"^ for r > l/v^r . (C75) 

To match this behavior with the large- r behavior of the instanton at d = 4, the constant C 
must behave as 

= 2ro(l + C(e)) with e = A-d. (C76) 

For fixed rg, we must let r 0, but at which rate? The evaluate this, use the equation 
(C72) which tells us that r / y?^ ~ £Sd=4: when e — >• 0. When evaluating f ip"^ in this limit, 
it is easy to see that it is the contribution of the domain ro <C r <^ 1/ -^r which dominates 
the integral, so that 

/POO r^/y/r 
p^ = Qd drr'^-V(r)' ^ 2 7r2 / drr'^-i (C r^"'^) ~ 2 tt^ (2ro)Mn(l/rov^) . 
Jo Jro 

(C77) 



108 



Therefore r goes indeed to zero as — 4 according to 

e = 4-d~6rr2 ln(l/rr2) . 



(C78) 



Now if we consider the polymer, we have to keep its length L = 1 fixed, hence J (f"^ = 1. 
Then the instanton size rg has to vanish together with r as e — > 0. 



1 = '-S 

T 



T ^ e- 



27r^MT7i) ■ 



(C79) 



Finally, we are interested in the smallest positive eigenvalue of the Hessian 5"/' and the 
corresponding eigenvector As e — * we expect that A"*" — * and — > (rV + l)(f the 
zero-mode for scale transformations. In this limit A^ can be estimated as follows 



(^1 



(C80) 



But the l.h.s. is equal to 



2 — d 1 

- T^ivlv) 



1 

27 



4r 2r 
Using the asymptotics obtained for if in the e — > limit 

(/?(r) ~ 2 ro =^ ^/^^(r) ~ -2 Tq for ro<^r<^l/^/T 

we obtain 

(^/;+|(^) ^ -((^|(^) = -1 (ip^lip^) ^ = 1 for £ . 

Hence the smallest positive eigenvalue of i^" vanishes as e when — 4, as expected 

47r2 



A^ 



2r 



for 6—^0. 



(C81) 



(C82) 



(C83) 



(C84) 



Thus in the limit d ^ A the Hessian Sj" gets an additional zero mode, so that the zero-mode 
subtracted determinant det' is discontinuous at d = 4 (limrf^4det' [S"] ^ det' [limrf^4 5/']), 
but we can write in the semiclassical estimates 



det' [5/'] 



.4 A+-det'[S/tJ 



(C85) 



The singular factor in eq. (C63), which comes from the integration over t, is canceled 
by the A"*" in det'5/', as expected, since we cannot have IR divergences in the semiclassical 
estimate at o? = 4. We get the IR-finite result 



1 4 

ImZ(r;L) = ^ - L"^ e"^ 



det'S"/' 



det'5" 



16 TT^ 

9 



(C86) 



where the IR singular terms coming from the dilation zero mode have disappeared. The UV 
divergences are contained in the two determinants det' [5"']. 
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6. Comparison SAM versus 0(n) field theory for the coefficient of the instanton 



We are now ready to check that the determinant factor for the instanton obtained by 
our method (non-local SAM model) is equal to the coefficient C86 obtained by instanton 
calculus in the 0(n = 0) local field theory. 

We have already checked in [14] that for D = 1 the instanton corresponds to the instanton 
for the 0^ field theory. 

For D = 1 the free energy density S[V] is nothing but the ground state energy Eq of a 
particle with unit mass in the potential V, i.e. the lowest eigenvalue Eq of the Hamiltonian 
operator 

H = + V{r) (C87) 

acting on functions over M'^. We denote ipQ the corresponding ground-state wave function. 



S[V] = Eo 



^0 



1 . 



(C88) 



The saddle-point equation is (using first order perturbation theory) 



V{r) 



\Mr)\' . 



6V{r) '"^"'dVir) 
So it can be written as the non-linear Schrodinger equation -|- constraint 



(C89) 







Eq such that / iJjq 



1 . 



(C90) 



This is equivalent to the saddle-point equation for the polymer instanton 

1 



- Aip + Tip - -ip"^ = 



T such that / ip 



-T . 



by the identification 

V^o(r) = V'(r) Eq 

In particular for L = 2 this gives ipQ = 1^/2, r = r' and Eq = r/2. 
Using second order perturbation theory we have 



(C91) 
(C92) 



6^E, 



6V{r,)5Vir2) 



2^o(ri)(ri| 



E<,-H 



where as in a previous section the "inverse prime" of an Hermitian operator 
inverse of this operator restricted to the subspace orthogonal to its kernel 



Eq-H 



1 



Pq = |^o)(^o| . 



Eq-H + Pq 

If we denote by 0o the operator which multiplies any function ip by ipo 

ipQ : ip ^ tjjQtjj , 



(C93) 
means the 

(C94) 
(C95) 
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we rewrite (C93) as 







(ri|2V'o 



1 



6Viri)5V{r2) '''' ^"[Eq-H 
The second derivative of the effective action F is thus 

1 



V'o|r2) • 



r" = 1 + 2iJo 



_Eq — H 

If there where no problems with the zero modes, we could write 

aAa ^ 9 / r 1 1 / 1 _ det(g-Eo-2V;g) _ det(-A/2 - E, - 3V>g) 
+ [^^J n " det(if - E,) - det(-A/2-Eo-^o^) 
quite similar to the ratio of determinants 

det (.^^0 
det {S'D 



but the zero modes require some care. Let 



^0 



A = H-Eo = -A/2-Eo 
B = H-Eo-2tpl = -A/2 -Eo- 3'ipl 



(C96) 



(C97) 



(C98) 



(C99) 



(ClOO) 
(ClOl) 



V'o is the zero mode of A and since dfj^A = B, = d^ipo are the d zero modes of B, while 
W^/i = i^odiMipo are the d zero modes of T" (as can be seen by using VF^ = —d^V/2, or by 
direct calculation). 

We use the following simple result. Let £" be a hermitian operator. If E has zero 
modes, let n = dim{Ker{E)) and $j a basis of Ker{E) and Kq = Y^- |<I>j)($j| the projector 
on Ker{E). Let E be another hermitian operator such that its restriction to Ker{E), 
F' = KqFKq is invertible. Then 

det[E + eF] = e"det'(E) det(F') with of course det(F') = det [(<l>i|F|<l>,)] . (C102) 

Now we consider 

r;' = r" + ei , (cios) 

and obviously 

det (r;') =e^det'(r') . (C104) 

Rewrite 

r: = 1 + e - 2^Po rr r^ i^o (C105) 



'H-Eo + aPo 

(this does not depend on the real number a). Then, since we now deal with invertible 
operators, we have 

Po 



det (r") = det 



det 



((1 + e){H-Eo + aPo) - 2^o'(l " ^o)) 



1 



H-Eo + aPo 



det [(1 + e){H- Eq) - 2V;g + ((1 + e)a + 2V^g)Po] 
det [H -Eo + a] 



(C106) 



Obviously 

det [H -Eo + a] = a det' [H - Eq] . 
Now consider the (non-hermitian) operator in the numerator 

B, = {1 + e){H - Eo) - + ((1 + e)« + 2^o')^o , 
which is not very different from the hermitian operator 

C, = {l + e){H-Eo)-2tp', . 
If ip is a vector orthogonal to Ker{H — Eq), i.e. {ip\ipo) = (or Pqi/j = 0) we have 

5,^ = C^ij . 

So the only difference between and is when applied to ipQ 

B.iPo = (1 + e)aV'o ; C,^o = -S^Ajj • 
In a basis of the eigenvectors tpi oi H — Eq, B^ and have respectively the form 



B^ 



(1 + e)a bj 
d 



„ _ , a bj 



with 



a' = (1 + e)a a = -2 / bi = -2 ip^ipi dij = (1 + e) 



Thus 
while 



det (B,) = (1 + e)adet(rfij) 



det (C) = (a - 6 . rf-i . 6*) det(rf,j) . 
Now B = Ce=o has d zero modes, the = d^j^ipo, hence 

■{V^\H-Eo\V,) 



det (a) = e'^det'(5) det 



We have 



and since {H - Eo)Vf, = 2i)lV^, 



{V,\H-E^\V,) = 2 I ^''od^^oduiJo = ^ /^o|WorV> 



since tpo is invariant by rotation. Hence 



det 



{V^\H-Eo\V,y 
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It remains to calculate the coefficient a — b ■ d ■ b . For this we use the fact that 



a-b-d-^-U 

which follows from (C112). Now a simple calculation shows that 

5(r-V^Ao + ^o) = (-A/2-^o-3V'o)(P- V^o + ^o) = SEqV'o 
and since ipQ is orthogonal to the kernel of B we can write 

1 



2Er 



(r- VV^o + ^o) 



(C120) 



(C121) 



(C122) 



Therefore, integrating by part and using the fact that HV^oP = J "^o ~ ^ obtain 



1 



lim — 

e=o a-b- rf-i ■ ¥ 2Eo 
Hence finally 

det'[r"] = 



2E. 



2-d 
1^ 



-I d 



V^o(r- V^o+^o) 



det'[H - Eo-2'4!^] 
det' [H - Eo] 



2-d 

4:Ea 



. (C123) 



(C124) 



Putting this result into (3.50) we get (using the fact that V = —'iPq] 



ImZ(fe) = / d'^ro 



V||VV^o 
ird 



2ld/2 



2-d' 
AE. 



-1/2 



-vs 



det'B 



det' A 



'1/2 



(C125) 



to be compared to (C86). 

For this remember that we are dealing with rescaled fields and couplings (with tildes). 
So we go back to the original variables by rescaling 



It gives for D = 1 



b\L] 



1 , , 2d _ d{ti-d) 

lmZ{b) = T-\b\~^L 2(d-2) 
2 



2-D 

2(D-e) 




V ^ 


1 

d-2 J, 






"iiv^oir" 


d/2 


'2 - d' 


nd 




[ 4^0 J 



D D<L 



4-d 
" d-2 



(C126) 
(C127) 



'1/2 2 4-d 

|b| d-2L d-2 r 



det'5 



det'A 



-1/2 



(C128) 



This is the same as (C86) since 
E, = -r , B = S'l , A = S'i , i:o = f , F = i?o + ^ ^ ^o' = - 2)^ (C129) 
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APPENDIX D: Usefull formulas for derivatives of traces and determinants 



To compute the Hessian matrix (5.39) of the variational energy £^var given by (5.28) we 
need to compute the matrix derivatives 



tr[Mf/2] and -..^ { det 



dMdM 



D-2 D-2 



'1/2 



with 



I, = i(M + 



the symmetrized of the matrix M, at the special value 



(Dl) 

Mvarll- To compute these 



derivatives it is useful to define the matrix 



(ij) 



Ci ® Cj as the matrix which on the line i 



and row j is 1 and is elsewhere^^, so that for any matrix A = {Ajj} 

dk „ dK 1/^ 



(ij) 



dAij 2 



Using this we can compute the first derivatives 



d 



OA: 



-tr[A"] = atr 



d 



OA, 



det [A] = det [A] tr 



-11 



and using the important formula 

the second derivatives for the trace 

-tr[A:] 



il) Ojk 



OA., dA 



aia — \] 



A"' {SikSji + Sjk6u) 



(D2) 



(D3) 



(D4) 



(D5) 



=A1 



and for the determinant (n being the dimension of the matrix, so that det [AH] = A^' 



d d 
dAij dAki 



[det(A"l + A^)]" 



{2A' 



16^2 



[a^6ij5ki + «(2 - a){6ii5jk + kk^ji)) (D6) 



APPENDIX E: Instanton condenstates 

Here we show that if ^(r) is the instanton potential, and & its action, we have the exact 
identities 



(V^(r)), 



((Vr)^) 



V 



V(<) 



-2I1--J 6 



(El) 



(E2) 



where the e.v. ( )y refers to the auxiliary model of a free (non-self-interacting) manifold 
trapped in the potential V^(r), with action 



Sv[r] = lj{Vrr + V{r) 



(E3) 
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The first equality in (El) follows from the instanton equation of motion (p) + ^ = and 
from 

(V^(r(xo)))v = jvir) {6{r - r(xo)))^ = ^V(r) (p(r))^ , (E4) 

while the second equality comes from a simple result of [14], rederived in Appendix F, see 
Eqs. (F6) and (F8). 

The first equality in (E2) follows from the equations of motion for the auxiliary model 
with action (E3) and the instanton equation. If we made the change of variable 

r(x) Ar(x) (E5) 

in the functional integral we obtain (up to contact terms proportional to 5^(0) which vanishe 
in dimensional regularization, and which correspond to the normal-product definition of the 
composite operator (Vr)^ = : (Vr)^ -.q) 

{{Vr)X + {r-Vy{r))y = . (E6) 

Now we can rewrite this second term as 

{vVyir))y = JrVrV{r) (5(r - r(xo)))^ (E7) 
and using the instanton equation and integrating by part we rewrite it as 

(r-V.y(r))^ = -^r.V.nr)xy(r) = -\ jj-V. {V{rf) = \ jv{rf . (E8) 
Q.E.D. 

Then we use (El) to obtain the second equality in (E2). 



APPENDIX F: A variational bound for the smallest (negative) eigenvalue 

In this section we derive a bound for the (negative) smallest eigenvalue A_ of the Hessian 
iS" which is associated to the unstable mode. The basic idea is as follows: The instability is 
visible by studing a rescaling of r and correspondingly x, V and £. The unstable mode has 
a non-vanishing overlap with this dilaton, which leads to a variational bound. 

First of all, we recall the rescaling 

r — > rx = \r (Fl) 
X — > xx = A^x (F2) 

2D 

V{r) — ^ Vx{r) = \—V{\r) . (F3) 
Under this rescaling the two terms of the effective action scale as 

2D 

£[V] — > £[Vx\ = X^S[V] (F4) 
J^[V] J^[Vx] = X^J'IV] . (F5) 

Now we consider the full effective action 



S[Vx]=S[Vx\+J'[Vx\ 



(F6) 
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The saddle-point equations for the instanton inforce 

which imphes 

The dilaton-mode is 

Note that it does not matter, due to (F7), of how one exactly defines the dilaton: one could 
use A^d^/dA^ instead. 
On the other hand 

with 

Hr) = X^Vr\r) . (FIX) 



Expanding in eigenmodes, 



v=v 



= V ■ e,) A, (e, ■ ^) > A^in (V^ ■ ^) • (F12) 

mst f ^ 



Therefore we have the exact bound 



^instl 



Using 

A^V5. = ^^y + rVy(r) (F14) 



one obtains the still exact bound 



r] 



This bound can of course not be calculated exactly, if we do not know exactly the instan- 
ton potential V . However, we can use the variational approximation for V^"^*^* to calculate 
the r.h.s. of (F15) approximately. Using for V{r) the Gaussian V{r) = exp(— (all 
normalizations and the width cancel at the end), one obtains 

\ ^ Avar _ - HP - e) . 

Amin _ ^min (3 - D) (2/} - e) + 6^ ' ^ 
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APPENDIX G: Normalization w.r.t. the variational mass mvar in the variational and 
post-variational calculations 

In this appendix we discuss the rescahng used in the variational and large-d calculations 
of Sect. V, where all quantities are expressed in units of the variational mass scale m. 
This rescaling is in fact quite simple and natural, but it might become confusing in some 
calculations, so we present it here carefully and thoroughly. 



1. The rescaling for x, r and g 

The variational mass m satisfies the equation (5.14), which amounts to 

m^-' = 2co(47rco)'^/2 , (Gl) 

where Cq = Co(-D) is the tadpole 




Co = (47r)-^/2r((2 - D)/2) = V V . (G2) 

As in [14] we perform the rescalings x — x in i^-space and r — r in d-space, with 

x = m-^x , p = mp ; r = m(^"2)/2)^ , k = m^^"^)/^^ (G3) 

in order to set the variational mass to unity m — > m = 1. In the new units the instanton 
potential V is and its Fourier transform are rescaled a.s V ^ V_ with 

V{r) = m^V{r_) . (G4) 

In addition we also redefine the measure over r in ci-space (and the corresponding measure 
over k in reciprocal space) as 



d'^r ^ J with J = m^^^ J d\ = y d^r (G5) 
with / = m^^' / = m^^ / - — — . (G6) 



7k Jk J (2vr)<^ J (27r) 

With this new measure the definition of the Fourier transform ^ of ^ in d-space is changed 
into 

£(k) = je-'^V{r) ; V{r) = ^e'^£(ls) (G7) 

and using (G4) the rescaling for the Fourier transform V of the potential V is V ^ V_ with 

V{k) = £(k) . (G8) 

Finally since the functional integration measure T>[V] over V is normalised by (Bl) which 
involves the measure over r and the effective coupling constant g, (this is equivalent to 
state that the metric G{6V,6V) = {-e-'yAng) l6V{rf over the space of V configurations 
depends on g and the measure over r), the rescaling of the r-integration measure (G5) 
amounts to a rescaling of the effective coupling constant g ^ g with 

9 = m^'g (G9) 
or equivalently of the original coupling constant b ^ b with 

b = m^~'b . (GIO) 
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2. Consequences 

a. Normalization for integrals and distributions 

With these normahzations all powers of m disappear in the variational and post- 
variational calculations, but we have to be careful when we perform Gaussian integrals. 
Indeed the following Gaussian integral gives 

-k^co 



e-^=o = 2co, (Gil) 

where cq is given by (Gl). Indeed, we have using (G6) and the equation (Gl) for the 
variational mass 

= m^-^(47rco)-'^/2 = m-^(47rG'^)-'^/2 _ rn-^2m^Gm = 2co . (G12) 
One has also to take into account the fact that the Dirac distribution in r space is now 

5(r) = m^^'5'^{r) such that j^5{r) = 1 . (G13) 

b. Action and Hessians 

Once this is done, all the results for the instanton and the large orders still hold without 
any factor m, in particular (3.51)-(3.53). The effective action S for the potential V is rescaled 
into S_[V] given simply by 



S[V] =£[V] + ^ fv^ , and is such that S[V] = m^S[V] , 



(G14) 



as well as its functional derivatives S'[V] = m^S^lV], S"[V] = m^S^'[V], etc. The instanton 
equation (3.41) is still 

£(k) + {e'^^-^)y = (G15) 
and the Hessian is still given by (4.2)- (4.3), i.e. (in reciprocal space) 

5" =1-0 with l,^,,^ = 5(xi-X2) , = -S" i.e. JZ] = /(e^^^^^^V^^^^^))^" 

(G16) 

The logarithm of the Hessian is now 

oo ^ 

£ = £ = log det'[5"] = tr log [1 - Q] = - V - tr \Q^] (G17) 

rh 

k=l 

with in particular 

Finally the zero-mode measure factor 2IJ is rescaled as expected 

d/2 

(G19) 



27id 
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c. Instanton 



In particular, this gives the variational instanton potential (obtained by replacing the e.v. 
in the instanton potential ( )y by e.v. in the quadratic potential ( )m=i (G15)) 

Vj^^\k) = -e-^'^''/^ and by Fourier transform Zlafd) = -2co 2^^/^ q-lV{2co) _ (Q2o) 

The instanton potential expanded in normal products w.r.t. the unit variational mass (i.e. 
: :=: :^=i) reads 

m-2.„±^[=lj,„:ier: (G21) 



= ^ n 7 / (-k'^TV(k)e-^'° . (G22) 

^ rf(c/+l)---(c/ + n- 1) A ~ 

d. Renormalized quantities and counterterms 

Finally let us see how the UV counterterms and the renormalized action transform under 
this rescaling. In Sect. IV.D.2 the one-loop counterterm AiiS for the effective action S was 
found to be given by (4.104) 

A.SIV] = -^i((Vrf>,-^i^V(r)^ 

and the renormalised effective action 5rcn[^] was 

5re„[V] = S[V]-g7^ifiLr^A^S[V] . 
If we now perform the rescalings it is easy to see that 

((Vxr)2>^ = m^((Vxr)2>-^ and ^^(0^ = ^Hl)' • (G23) 

We define the rescaled renormalized couplings and g as for the bare couplings (G9)-(G10) 

= m^g^ and b, = m^^% . (G24) 
Then the rescaled renormalized effective action defined by 

Sren[V]=m''S^^^[V] (G25) 

is given by 

^.en[Z] = m-aJ^W-^^lSm (G26) 

with the rescaled one-loop counterterm Ai5[V] given by 

AiS[V] = m'AiS[V] (G27) 
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and using (G23) we write AiiS[]/] as 

AiS[V] = 
with the rescaled counterterms 



= "^Ci and Cj = m 



D-e, 



(G28) 



(G29) 



Now we use the exphcit perturbative results (4.85)-(4.87) for the counterterms Ci and C2 
and Eq. (Gl) for the variational mass m and obtain for the counterterms C^^ and 



Co 

do 



1 1+^ 



2Sd T[D/i2-D)]^ 
{2 - Dy r[2D / {2 - D)] 



Co 

do 



1 1+^ 



2D 

where we remind that 

Sd = 2 7r^/Vr[D/2] , Co/do = -2''''r[i2 - D)/2]/r[iD - 2)/2] 
and in the limit d ^ 00, e fixed, is of order 0{1) since 

Ci = -TT 2^-' e-(^-')^^ [1 + 0{l/d)] , 7e Euler's constant . 
while Cj is exponentially small. 



(G30) 



(G31) 



(G32) 



3. Final results 

With these notations, the final results for the large orders have the same form, with the 
unrescaled quantities replaced by the rescales ones. 

In the bulk of the paper, when we use these normalisations, we rely on (Gil) and (G21) 
and omit the underlinings _±_ for all the quantities and the fields such as x, r, V, g, S etc. 



